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body of kinetic evidence which specilies N,0, 


Wo HNO, svstem sugyvested that the 
Apparent 


equilibrium data for this system would be more logically expressed in terms of 0, 


temperature inicle ym ndence is achieved by this manner of correlation 
ionization of N,O, in solutions are in accord with the 


nature of the reaction would seem to precluck 


between these nitrogen oxides and water 


Résume 
dans le svstéme NO, No, 


les donnees ce Péquilibre Cette methode ce 


apparente ce la temperature 


solution, sont en accord avec les données cinetiques et cle Peéquilibre 


tne étude des nombreux résultats cinétiques indiquant que N est le re 


Les resultats de la littérature 


Literature data proving the 
| 
kinetic and equilibrium data rhe ion 


the existence of i purely yus-phiase reaction 


etil ipal 


HNO, i sugvere Lutilisation ke pour exprimer 


orrelation permet dobtentr um independance 
prouvant Pionisation de re 


Le cauractére de la 


réaction semblerait exclure existence d'une réaction purement en phase gazeuse entre ces oxvdes 


et Peau. 


Pazote 


Zusammenfassung 
Reaktanten inn svstem NO). N,0, 


Gleichgewichtsdaten dieses Systems mit mehr Berechtigung 1 
sind. Man erhdlt bei dieser Zuordnung cine scheinbare Temperaturunabhangigkeit 


in Lésungen sind mit den kinetischen und Gleichgewichts 


werte tiber die Lonisation von N O, 


werten im 


Iv would be dillicult to cite a process so common 
to the industrial chemical scene and vet posing 
so severe a challenge to our compre hension of its 
fundamental mechanisms as that of the absorp 
tion of nitrogen oxides in water to produce 
nitric acid. The reactions involved are recognized 
to be the vapour phase oxidation of nitric oxide 
by oxygen followed by absorption and reaction 
of the nitrogen dioxide in the aqueous phase to 


produce nitric acid and gaseous nitric oxide. 


The nitric oxide oxidation 


Aus ciner kinetischen Betrachtung 


bereinstimmung. Der lonen-Charakter der Reaktion scheint dic 


reinen Gasph wenreaktion zwischen diesen Stickoxvden und Wasser auszuschliessen 


welche NO, ils den mass gebenden 


kennzeichnet folyt dass di 


HINO, 
n Ausdriicken von anzugeben 


L.iteratur- 


Annahme einer 


2NO + 0, —>2 NO, (1) 


has been the subject of rigorous study for many 


initiated this work, 


decades. BopENSTEIN 


revealing the third-order kinetics and the well- 
known negative temperature coetlicient of this 


reaction, Thirty vears later. Tri acy and 
Danxrets [19] postulated that the mechanism of 


Reaction (1) Is 
NO 
NO, NO 


> NO, (2) 
-2 NO, (3) 
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where Reaction (3) is rate-controlling. 
the negative tempe rature cde yx ndenev and over- 
kinetics of the 
apparently rationalized. 

With the 


dioxide and 


all third-order reaction are 


respect to reaction of nitrogen 


with water aqueous 
nitric acid, there has been evidenced an equally 
intense interest. Chemical equilibrium has been 
al, [2], Burpick 


SHERWOOD 


‘ 


determined by Ane. et 
Freep [4], 
Denpicn and Prince [9], and Ersurrein [10). 


CHAMBERS and 


The reaction has usually been written as 


2 NO, H,O —— HNO, HNO, (4) 
and 
3 HNO, ~— HNO, H,O + 2 NO (5) 
The over-all reaction becomes 
3NO, + HNO, NO 
At equilibrium, therefore 
(7) 
P’xo 
where K, is the thermodynamic equilibrium con- 


stant which is conveniently expressed as the pro- 


Fic. 1. Nitric acid - NO — NO, equilibrium. 


Thus, 


solutions of 


CARBERRY 


duct of A, 
In the 


ments, the above equilibrium is commonly ex- 


"HNO. / 


activity 


PNo and AK, 


absence of reliabk measure- 


pressed as 


log ie A, vs. weight per cent. HNOg. 

: 
The data of numerous workers are shown in 


this form in Fig. 1. 


ABSORPTION KINETICS 


The kinetics of absorption have been subjected 
to investigation principally by Cuampers and 
SuERwoop [7], Densicn and Prince [9], Caupi 
6], Pereers [16], and more recently 
by Wenper and Picrorp [20], 
failed to into 
First, while the 


and Denricn 
Interpretations 
have conn avreement two 
work of Chambers and 


diffusion 


as rate-controlling, Dexpian and Prince believed 


respects. 


Sherwood supports gas-film step 


that chemical reaction in the liquid phase controls 


the over-all absorption. The more recent work 


of and appears to reconcile 
the issue in that it was shown that the gas rate 
did affect the absorption insofar as local eddying 
and mixing in the reactant liquid phase is a 
function of gas rate. 

The second point of disagreement pertains to 
the 


water vapour and nitrogen dioxide, 


existence of a gas-phase reaction between 
There are 
compelling arguments both for and against the 
CHAMBERS 
[7], Perers [16] and McHanry 


[15] submit evidence in favour of 


J 


gas-phase production of nitric acid. 
and SHuerwoop 
a gas-phase 
reaction occurring simultancously with the liquid- 


WENDEI 


issue at some length and concludes that the data 


phase process. [21] has discussed the 


supporting the gas-phase reaction are indirect 
and that inferences drawn from the existence of 
nitric acid mists, gas-phase exotherms, and the 
generation of NO in an alkali system are tentative 
since other mechanisms can be postulated to 
account for these observations. 

On the other hand, [18], Kuz’winyKnu 
[13], and the author [5] found that no reaction 
occurs between NO, (N,O,) and water vapour 
in experiments conducted in the absence of a 


liquid phase. It may be significant that support 
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of a gas-phase reaction is found in studies where 
a condensed phase exists, while data taken in 
completely gaseous systems do not indicate 
reactivity of NO, and water. 

Not unrelated to this issue is the problem of 
specifying which one or more of the nitrogen 
oxides N.O,. N,O3) are the reactant 
species in the absorption process. On the basis 
of the general equilibrium correlation (equation 7) 
one would infer that NO, is the reactant in the 
absorption step. However, a vast body of 
chemical kinetic evidence dictates otherwise. 
Denpicu [9] found a negative temperature 
dependence of rate constants evaluated in terms 
of NO,. However, positive temperature depen- 
dency was displayed when N,O, was considered 
as the reacting agent. and Denpicu 
[6] and Wenpet [21] have shown that the rate of 
nitrogen oxides absorption in water increases 
linearly with the partial pressure of N,O, in 
the gas. Since 

2 NO, = N,O, (8) 


DeNBIGH concluded that two elementary 


reactions are involved in the absorption process : 


N,O, + H,O => HNO, HNO, 


2 HNO, = NO + NO, H,O (10) 


The rate data were correlated by Densicu with 


the equation 


d (N,O,) 


= k NO ) 
dt 


C (N,O,)'* (NO)? ] (11) 
while the experimental equilibrium data were 
expressed in terms of NO, as shown in Fig. 1. 
These kinetic data certainly suggest that NO, 
is the key reactant in the NO, — H,O system. 
The reactive nature of N,O, relative to its 
monomer NQ, is indicated in systems other 
than those involving water. In a study of the 
kinetics of the reactions of nitrogen dioxide with 
various alcohols, Farnum et al, [11] found that 
third-order kinetics and a negative temperature 
coeflicient characterized the initial reaction rate. 
When, however, their rate equation was written 
in terms of N,O,, second-order kinetics and a 
normal temperature dependency resulted, They 
concluded that N,O,, not NO,, is the reactant 
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species with alcohols. Subsequent work [8] on 
the same system demonstrated that this reaction 
is heterogeneous in nature involving the reaction 
of NO, and adsorbed alcohol on polar surfaces. 

With this imposing body of evidence favouring 
N,O, as the chief reactant candidate in various 
systems, it appears that equilibrium data for the 
nitrogen oxides — water system would be more 
logically expressed in terms of NO, rather than 
NO,. Given the dimerization equilibrium constant 
for Reaction (3) 


(12) 


then the nitrogen oxides — acid equilibrium ex- 
pre ssed in terms of N,0,, NO, HO, and HNO, 


is 
, PNo 
A 
15 
PNO K kK, 
K, — K,K, 


Fic. 2 Nitric acid -NO — N,O, equilibrium. 
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where 


K, = K,/K;"* (13) 


K, P*®xo 


Therefore, the equilibrium data of Fig. 1 may b 
ding the 


NSTEIN [3] 


stated in terms of NO, by simply div 
publish d values of kK, by Bop 
measured A, for Reaction (8) over a temperature 
range of 9-113°C. Emploving B 
data, all the data in Fig. 1 have bee 


accordingly and the resulting correlativ 


NSTEIN'S 
translated 
m appears 
in Fig. 2. 

Evidently, within the limits of the experimental 
errors involved in the original equilibrium data, 
the new correlation reveals apparent tempera 
ture independence in the range of 10-55 per cent 
by wt. HNO,, The 
dilute ( 


per cent) regions of acid 


spread of the data in the 
10 per cent) and concentrated 55 
strength is beyond 
comment since the original data in those regions 
are not experimental but extrapolated values*. 
The correlation presented in Fig. 2 would 
suggest that one of the chief factors « ntributing 
to the reduction in equilibrium concentration of 
HNO, with increasing temperature is indeed the 
N,O, in equilibrium wit! NO, at 


levels. Certainly a 


paucity of 
the higher 


dynamic equilibrium model stated ix 


temperature 
terms of 
N,f VP is in accord with the kinetic evidence cited 
above. 

Given the NO,-N,O, equilibrium 
constant as a function of temperature, it follows 
that the correlation in Fig. 2 may be utilized to 
predict acid equilibrium in temperature regions 


pre sently devoid of experi ntal data. 


REACTION MECHANISM 


While the kinetic studies and Implic ations of 
the correlation in Fig. 2 signify the key role of 


N,O, in the reaction under discussion, further 


*Subsequent to the initial preparation of this paper 
it has been brought to the author's attention that 
R. L. Menecus in an unpublished report, in 1943 
employed an N,0, correlation of acid equilibrium as a 
basis for extrapolating values into the range of acid 


concentrations greater than 60 per cent by wt 


evidence of even greater import is to be found in 
studies devoted to NO, behaviour in solutions. 
Gray and Yorr [12], in their review of 1955 on 
the reactivity and structure of nitrogen dioxide, 
cite and discuss a wealth of evidence for the 


ionization of N,O, according to the reaction 
N,O, = NO* 4 NO, (14) 


This ionization occurs even in pure liquid N,O,, 
and solvents of high dielectric constant (e.g.. 
Reaction (14), 
17] attribute the conductance 
HNO, -H,O-NO, (N,O,) system 
to the presence of NO* and NO,-. Ve ry recently, 
Miter and Warson [14] found that 
centrations up to on molal N04. the concentra 
tion of undissociated N,O, is small.” Further, the 
concentration of NO, does not exceed 0-02 molal 
and 


The reaction seq uence involved In the 


water) obviously enhance 
RoBERTSON et al. 


vaues of the 


“in con 


im coneentrated solutions of nitric 
acid 14). 
absorption and reaction of NO, with water may 


be writen as 


2 (N,O,),... (15a) 

(N30 4) (Ne NO" NO, 
(15b) 

2 HNO, = H,O + NO + NO, (15c) 


Following GRAY and Yorri [12], the rate of 


reaction of NO, with water can be written 


d(N,O 
t {ky (N,O,) (H,0 
dt 
(EENO,) | (16) 
Sines 
(HNQO,) N,O,)! *(NO)! (17) 
then 
d(N,O 
3 4k, (H,O) 
dt 
k, GINO 
(N,O,) — « 4 (18) 
ky (H,O)'? 
k [(N,O,) — (19) 


The above expression is identical with that 


found experimentally by Densicnu (equation 11), 
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Densicn and Prince [9] first demonstrated 
accord between their experimental correlation 
(equation 11) and the theoretical rate expression 
derived from Reactions (9) and (10). Equation 
(19) oby iously results from a 
and (10), 
alter the 
theoretical rate expression. 


only Reactions (9) 
of the 


It was the idk ntity 


ionization does not final form 


of the expr rimental and theoretical expressions 
which led Densicn and Prince [9] to the 
that 
diffusion, 


than 


abse ion of 


reaction, rather 
controlled the 
into water and nitric acid, 

It is of interest to note that although DENBIGH 


and Princ correlated their e uilibrium data 
| 


conclusion chemical 
gas-tilm 


nitrogen dioxidck 


in terms of NO,,. as shown in Fig. 1, an examina 
tion of the constant C in their rate expression 
(equation 11) reveals that 


(RT) *¢ 


In other virtual 


inde pe ndence of the pres nt equilibrium correla 


words, the temperature 
tion (Fig. 2) was implied in Denrnicu’s experi- 
mental rate expression, 

\ possible explanation for the unusual equill- 
brium correlation is suggested by examination of 
equation (18) in the lisht oftheionic mechanism. At 
equilibrium expressed in terms of partial pressure 
oft th N,O, and NO, it is that the 


‘equilibrium constant,” contains solubility 


evident 


equilibrium constants for the gases as well as an 
ionization constant and the appropriate activity 
Thus. 


dependencies of the various physical and chemical 


coctlicrents, compensating temperature 
equilibria would seem to account for insignificant 


temperature effect upon over-all equilibrium 
expressed in terms of N,O,. 

Clearly Fig. 2. does not prove N,O, 
to be the chief Primarily, the kinetic 
and ionization data provide the chief support 
for the reactant candidacy of N,O,. It should be 
noted, however, that while the kinetic data cited 
specify N,O, rather than NO, as reactant, in the 


absence of knowledge of the NO, ionization, the 


jt r 


reactant. 


locale (vas and or liquid phase) of the reaction 
remains open to de bate. 


As the ionic nature of the NO, water reaction 


is evident, one is compelled to make the inevitable 
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consideration of 


evidence of 
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water system 


that 


are quite unlikely, a purely gas-phase reaction 


suggestion since gas-phase ionic reactions 
between NO, (N,O,) and water is implausible. 
This preclude 
between gaseous NO, (N,O,) and water droplets 


contention does not reaction 


or a water fog. Such an event would be truly 
heterogeneous and in accord with the above 
mechanistic picture. It cannot be overem- 


phasized that, with the exce ption of McHaney’s 
data 


a vapour-phase reaction were 


work [15], the reliabk of those who infer 
the 


obtained in 


existence of 


svstems containing a condensed 


aqueous phase, whil experimenters who have 


worked with gaseous systems find no evidence of 
a gas-phase generation of acid from NO, (N,O,) 
and 


water. KUZ MINYKH [14] performed his 


experiments under both conditions and _ his 


results verify this generalization. Obviously, 
only a few droplets of water need exist in the gas 
phase to constitute acid and NO veneration with 
the possible conse quent formation of a fog or mist. 

The role of NO in absorption kinetics requires 
comment at this point. That its presence in a 
nitrogen oxides gas enhances the rate of absorp- 
tion in water, acid, and alkali has been attributed 


to the re action 


NO NO, N,0, (20) 
The species N,f ds is the anhydride of nitrous 
acid as N,O, is the anhydride of nitrie acid. 


Whi thermodynamics indicates that only small 
amounts of NO, exist in equilibrium with NO 
and NO, at 


attainment of equilibrium is quite rapid, 


process temperatures, the rate of 
Reac- 
tion of N,O, with alkali yields primarily the 
nitrite salt, while in the absence of NO the alkali 
system produces nearly equimolar quantities of 
nitrite and nitrate salts. The reaction in alkali is 


expressed by the reactions 
N,f )s 
2 HNO, 


20H 


-2 HNO, 
2 NO,- 


(21) 


(22) 


-2 H,0 


The formation of nitrous acid from its anhydride, 
N,O . by Reaction (21) is known to oceur quite 
rapidly in the gas as well as liquid phase [12]. 
lonization in the case of the N,0, H,O system 
is apparently the cons quence of OH™- ion attack 
upon the HNO, formed via the hydration reaction. 
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The direct evidence for the rapid gas-phase icknowledgement— The author is indebted to Dr. T. H. 
(21) Chilton for his most valuable advice during the preparation 


formation of nitrous acid by Reaction 
of this manuscript. 
would suggest that perhaps the mist or fog es 
encountered in some nitrogen peroxide absorption . 
NOTATION 

studies is nitrous rather than nitric acid. 

“uNO activity of nitric acid in solution 

CONCLUSIONS 44.9 = activity of water in solution 

constant in equations (11) and (19) 
The results of numeorus kinetic studies of the 
rate constant in equations (11) and (19) 
nitrogen oxides-water system clearly indicate , equilibrium constant for Reactions (9) 
that NO, is the chief reactant species in this and (10) 
system. This suggestion prompted a correlation : equilibrium constant for reaction (6) 
of the nitrogen oxides — aqueous acid equilibrium PROP’ NO, 
data in terms of N,O,, rather than NO,, with @" HNO, “HO 
the result that appare nt tempe rature ind pen- equilibrium constant for Reaction (8) 


dence of the equilibrium data is realized. Phy sical 


15 
A, 


measurements revealing the ionization of dissolved : 
partial pressure of NO, 


N,0, into NO* and NO,~ ions in aqueous media 


partial pressure of N,O 
not only support the kinetic and chemical i 


partial pressure of NO 

equilibrium data but provide potent support to 
it law constant 

the argument that a purely gas-phase reaction T bsolute temperature 

between NO, and water Vapour Is an unlikely denotes concentration of enclosed species 


event. : constant in equation (17) 
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Abstract 


Graphical method for design of cooler condensers by MizusHiIna and Koroo \9 


is extended to general gas vapour systems including air water having high \ apour concentrations. 


Moreover this is applicable regardless of whether the mixture continues saturated or superheated, 


and also simpler than other methods developed previously though it may be less rigorous. 


Two examples of air-water system and one example of air—benzene system are solved by 


this method. It is shown that there are good agreement between the results of this method and 


those of more rigorous ones, i.e. CounpuRN and Hougen [6 method or Bras [1] modification. 


Résumé Les auteurs ¢tendent au svstéme général gaz-vapeur comprenant air-eau, a forte 


concentration en vapeur, la méthode graphique de Mizusui~a et Koroo (9) pour le calcul d'un 


condenseur Cette méthode est applicable indifféremment si k melange se maintient saturé 


ou surchauffe, et elle est plus simple que d'autres méthodes décrites précédemment bien qu'elle 


soil moms 


Deux exemples du systéme air-eau, et un exen ple du systéme air-benzéne, sont résolus 


par cette method Les auteurs montrent quil vy a concordance entre les résultats de cette 


methode et ceux de méthodes plus rigoureuses par exemple la méthode de CoLpurNn et Hovces 


(} ou la modification de Bras [1 


Zusammenfassung —Die graphische Methode fiir den Entwurf von Kiibler Kondensatoren 


nach Mizusmina und Kovroo [9] wird auf allvemeine Gas Dampf-Systeme einschliesslich Luft- 


Wasser mit hohen Dampfkonzentrationen ausgedehnt. Ferner ist diese Methode sowohl fiir 


gesittigte wie tiberhitzte Systeme anwendbar und ist auch einfacher als friiher entwickelte 


Verfahren, wenn sie auch weniger exakt ist 


Zwei Beispiele des Luft-Wasser-Systems und ein Be ispiel des Luft-Benzol-Systems werden 


nach dieser Methode behandelt. Ks ergibt sich eine cute (ber instimmung mit den Ergebnissen 


exakter Methoden, z.B. nach Cotsurn and Houcen [6] oder auch nach der modifizierten Methode 


von Bras 


INTRODUCTION 


general method suitable even in case of super- 


IN a cooler condenser, vapour condenses at the heating, the tediousness of the stepwise method 


cooling surface and the te ripe rature of ras IS still Its disadvantage. COLBt RN [5] in a recent 


vapour mixture decreases. The simultaneous P@per emphasized the necessity of a simplified 
transfer of heat and mass makes the design ™ethod of the calculation and proposed a method 


problem complicated. Since Cotsurnx and based on the terminal conditions only. Its result 


HouceEn [6 published their paper on this subject is, however, sometimes far from the rigorous 


first, various methods wer proposed, Colburn value as shown by Carrns [4]. 


and Hougen’s stepwise calculation is the most For simpler calculation which can be applied 


rigorous but limited to saturated vas Vapour regardless of whether the mixture is saturated 


nuxtures. Though this disadvantage was or not, graphical methods may be most hopeful, 


removed by Bas [1] who proposed a new and have been developed by several authors. 


Hitachi Ltd., 
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T. Mizusmrea, 


MickLey’s method [8] can be applied only 


water vapour mixture having 


to the case of air 
low vapour content. 
Just MickLey’s paper was published, 


Mizusutna and Koroo [9] proposed a graphical 


be re 


method which is useful for general systems of 


gas—vapour mixture but with low vapour concen- 
trations. 
graphical 


NELSON (7 


calculation method which operates s iwcessfully 


( RIBB and developed 


for a mixture even with high vapour 


concentration. However, it Is teclious t construct 


and use two charts, Le. temperatur enthalpy 


and temperature humidity charts. Thev assumed 


that the ratio of heat 


transfer coethcent 


nt to mass 
heat 
the 


transfer cocth 


is constant becaust the 


capa itv of the gas mixture is constant over 


This 


particular svstem used by 


range of humidity sumption ts 


the 


authors but may 


whole 
thre 


thei 


valid for 


not tor rs, Moy 


assumption on mass transfer coectl nt beme 


constant over thre whol cooling Is quite 


doubtful when the mass velocity of gas mixture 


im the 
sequence ol Vapour condensation. 
Recently, Bras [2, 8 


procedures for the determination 


apparatus changes considerably mn con 


developed graphical 
ot tie 
which ti nad error 


i ver it ts 


should by 


vas interface conditions, by 


computations were elimin ited. 


still tedious that values of propertr 
obtained by point to point caleculat 

The graphical method which will be deve loped 
in this paper is simpler and applicable to any 


gas vapour mixture having vapour 
concentration though 1t less 


Moreover, this method can handle tin 


svstem ol 
may bn rigorous, 
situation 
saturated. Essen 


Miz 


is 


even when the mixture ts not 


tially this is an extension of SHINA and 


Koroo's graphical method [9 


to the case of low vapour concentration, 


limited 
to that 


of high vapour concentration, 


OF A GRAPHICAI 


WATER 


THE 
METHOD 
Basi equations 


ratio of heat transfer coeflicient to nm 


DEVELOPMENT 
AIR 


FOR SYSTEM 


the 


ss transter 


counter current cook r condense r. 


coethicient ts 


HAsuIMorTo 


and M. NAKAJIMA 


h 


Pam | 


relation between k, and k,’ is 


(ps (1. HjdA 


M 
(P 
equation 
k,’ M 
p 
| P 
M, \(P —p,)(P 


Substituting H p)| (M My») into the 


above 


From equations (1) and (2) one obtains 


h, Se\} 


p \(pP pol 
where 


Pum | 4 
cM, \(P Py) (P 


pi ) 


Se \? 
P — Pr) 


(4) 


Even in the case of high vapour content, the 


p,) is ordinarily 0-8 ~ 1, and 
Pry? taken 
Therefore, the 


Of Paw P 
AY 


value 


may as Ot for 


the value of 


air-water svstem. value of « 


approximates to unity for aimr-water system. 
This approximation its conservative. 


Hk ret 


re gardl ss of Vapour concentration, 


This is an extension of Lewis’ relation. 
Therefore, the 
equation is applicable to 
the 


vapour. 


following conventional enthalpy 
a cooler condenser of 


amount of water- 


humid air containing any 


= (6) 


Cc 


Gy di, i,jdA 


where, the enthalphy of humid air is defined as 


*In the case of gas system [10], the power to Schmidt 


amd Prandtl modulus is } rather than 4. 


a, 
k, 
Pi 
(2) 
VOL 
ae, (3) 9 
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The surface area is obtained by where 


(dG,'\" 
c, di JB ( | ") (18) 
A Gy, (3) 
’ Assuming the values of viscosity and Prandtl 


number of the stagnant gas ; ame as tl f 
Heat balance gives the inclination of operating 


the mixture respectively, and comparing equa- 


line. 4 Le tions (16), (17) and (18), the following equations 
(9) are obtained. 
dt, Gr h hy, 
‘ 
The inclination of path line of air condition is C, Bey oe 


where 


(10) (20) 


Substituting equations (19) into equations (8) 


and (11), the following equations are obtained. 


(11) 
4 ai 


h, 
t; by 4 Gp 
hy i 


Consideration of the change of values of h, and « 


in the condenser. 8 di 
If there is a large amount of the change of the hp (; i.) 
values of h, and ¢, from the inlet to the outlet a 
of a condenser owing to the condensation of i. F he Cp 
(22) 
the vapour, it is desirable to take them into t, —t, hy 


account. 


For gas systems, heat transfer coeflicient may a 
be expressed as follows. | 
h, je, @ (Pr) (12)* 
where 
16} 
_{dG'\" 
( | (13) 


The relations between the quantities on dry 


basis and those on total basis are 


G = G,' (1 + H,) (14) 


i.) (15) 
Introducing equations (14) and (15) into equations 
(12) and (13), one obtains HUMIDITY 
iG..’\ Fic. 1. Relations between 8 and humidity and humid 
h, (1 i1,) (Pr)} (16) heat and humidity. 


Detining hy as fictitious heat transfer coefficient Design procedure 


for the flow of stagnant gas only gives The values of 8 are plotted against the 
“ae? : humidity H, on Fig. 1, » 0-4 is for the flow 

hy ip Gy (Pr) (17) 4 
' across tube banks and » 0-2 is for the flow 


*Refer to the footnote on page 196. inside tubes. c¢, vs. i, are also plotted on Fig. 1. 


dt t, tj 
9 
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TEMPERATURE T - 
Fic. 2. Enthalpy temperature chart 
Fig. 2 is the enthalpy-temperature chart with 
a group of constant humidity lines. 
If it is that the of inlet 


gas, the enthalpy of outlet gas, the te mperature 


assumed condition 


water and gas flow 


of water available and the 
rates are known, the procedure for the graphical 
method is as follows. 

(1) Caleulate the heat 
gas film h. defined by equation (12) at the gas 
the 


c,/h 


transfer coetlicient of 


inlet. Knowing the values of c, and 6 at 

inlet gas condition from Fig. 1, 

by equation (19). 
Also, calculate the 


he. made up of the conductances of th 


combined conductance, 
conde nsate 
film, metal pipe, scale and the cooling water film. 


water 


In most Vapour 


practical cases of air 
mixture, hi, can be taken as constant. 

(2) Plot the operating line 1Z on the enthalpy 
temperature chart. 

(3) Subdivide the operating line into several 
equal increments (4.B.C...Z). 

(4) Plot point A’ (t,,, 

(5) Construct a tie line of slope Bho Cy/hp 
A (ty, i) and locate A’ 
The value of 


through the point 
(t,,, on the saturation curve. 
8 can be obtained from Fig. 1 for the humidity 
(6) Join 
A (ti, and intersect 


B parallel to the abscissa at 


point A’(t.,, with 
a line through point 


point BR’. BR 


represents the condition (f,, of the second 


point 


point. 


Hasumoro and M. 


NAKAJIMA 


(7) Repeat steps (5) and (6) for the remaining 
increments. 8 should be obtained for the humidity 
ABC .Z is the path line 


of the gas condition. 


of each point. 


(8) If the gas becomes saturated at any point, 
the 


the saturation curve. 


gas condition changes hereafter following 


(9) Plot £/(i, i;) vs. and obtain 


‘21 
(i, i,) graphically. 


(10) Calculate the value of A by 
(21). 


equation 


The calculation results and the comparison of them 
with the results of the other methods 


(a) Example of Cotmurn and Hovcer 


(saturated mixture). Using the procedure men- 


tioned above, the example in Colburn and 


Hougen’s r was solved. 


¢, O41 P.c.u./Tb °C 
G 10,800 |b /hr ft? 
Pr 0-76 
19.600 
j 0-0063 
All of these values are taken from Colburn and 
Hougen’s paper. 
Substituting these 
the heat 


calculated. 


values into equation (12), 


transfer coellicients of gas film, h, is 


h, (0-0063) (0-41) (10,800) /(0-76)! 


31-9 P.c.u. ft® hr 


Humidity H, is calculated from vapour pressure. 


(O-S35) (18) 
3-14 Ib Ib 
(0-165) (29) 


Hence, the value of ¢, is obtained by Fig. 1 as 
1-69 °C 


Since the gas flows across tube banks, the value 


of 8 is obtained by n 0-4 curve in Fig. 1. 


p 1-765 


| 
A 
| 
| 
D4 ao 
L 
| 
_ 9 
1958/ 
195 


From equation (19) 


1-69 
(1-765) (81-9) 


CB Cs 


Bh, 


0-03 
hy 


Colburn and Hougen gave the value of the 
combined conductance as 

hy = 310 P.c.u./hr °C 
Since the gas condition at the inlet is ¢,, 95 ( 


and 


at this point is calculated as follows. 


3-14, the enthalpy of the gas mixturé 


(1-69) (95) (595) (3-14) 
2,030 P.c.u./Ib 


Since the mixture at the outlet is saturated at 


tag = 40°C, 


0-049 Ib Ib 
0-262 P.c.u. Ib 
10 P.c.u. Ib 


The cooling water temperatures are shown in 
Colburn and Hougen’s paper as follows, 
At the gas inlet 


th, 60°C 


At the gas outlet 
te = 25°C 


Thus, the operating line can be plotted on the 


enthalpy-temperature chart. 


Since the gas is saturated with vapour from 


Table 1. Factors 
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the beginning in this example, the path line of 


determined from the graphical method on CoLwurN and Hovucen example 


the gas condition coincide with the saturation 


curve, 


the mentioned above, 


results tabulated in Table 1 


Following procedure 


the are obtained. 


And the surface area of the condenser is calculated 


by equation (21) 


igi 


Gin Cp di, 


hp i 


Gy» 2-800 lb hr 


B 0-03 
‘BR 


di 


> 


A (2-800) (0-038) (8-295) 696 


Hougen’s result is 695 ft?. 


both 


The Colburn and 
The 


apparent. 


good agreement between results is 


In Fig. 3, the te mperatures of gas, condensate 


surface and water are plotted against the surface 
The 
results, and the broken 


area of the condenser. solid 


lines are for 


CoLBURN and HouGEnN 


lines are for the results of this method. This 
again shows the comparatively good agreement 
between both methods. Enthalpy difference 
between the two points at higher temperature 


transferred 
the 
gas mixture owing to 


the amount of heat 


condensate 


is larger than 


the 


same temperature as the 


where is assumed to be at 


the fact that the enthalpy is based upon 0°C 


Point i, iy ty B 
(No.) (°C) (Ib/Ib)  (P.e.u. /Ib) 


l 05° 314 2030 
2 93-2 2-20 1421 1-2 1-593 
90-0 1-40 1-419 
0-80 528 33-6 1-265 
5 5 270 1- 
30 206 1- 
20 i”) 
8 510) 0-10 75 25-6 1-039 
16-5 0-07 55 1427 
10 wo O-O49 ) 250 


16-45 1480 3-21 696 
14-83 S32 S80 2-70 550 
13-21 80-5 380 2-67 £36 
11-79 65-6 148 380 3°33 
10-65 j 5- 


10-35 


10-02 10-25 

9-68 30-8 25 20-80 
O56 28-6 32-5 31-60 53 
9-49 270 20-5 19-5 52-30 0 


10 1 
hp (°C) (P.c.u. /Ib) (P.c.u./Ib) B/G, - 


CR t i i i 


6-98 


: 
= 
VOL. 
iga 
9 
— 
‘ 
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0-00385 
11-0 hr ft? °C 
0-11 Ib Ib 


0-272 Ib °C 


In this case the gas flows in tubes, therefore 
n 0-2. 


1-021 
0-0242 


113 P.c.u. /hr ft? 
105-5 P.c.u. /ib 


23 


At the gas outlet 


Fic. 3. Gas temperature, ¢,, interface temperature, { 
and cooling water temperature, ¢,; vs. surface area of a 


condenser. The calculation results are shown in Table 2. 


The value of integration is 8-31, and 

water. Therefore, th temperature decreasing 
rate at higher temperature of this calculation is Gy 7,050 Ib ‘hr 
lower than that of Cotpurn and Hovucen. The surface area is. therefore 

However, water temperature distribution lines 
of both Thit thods coincide wh other « losely. A (7,050) (0-0242) (3-31) 565 Tt? 
This means that the calculated amount of heat The result of Bras is 561 ft?. 
transferred across a certain range of the surface The difference is only 0-7 per cent. 
area is the same for both methods. 
A Grapuicat Meruop ror Varour-Gas 

Mixture Oruer THAN 


SYSTEM 


(b) Example of Bras [1] (unsaturated miature). 
At the gas inlet. 


‘ 0-245 P.cu. Ib Basic equations 


G 11,200 Ib hr ft? Equations (3) and (4) can be applied to this 
O76 case too. However. (Sc is not equal to 


Re 14,600 unity in this case. However, it may be assumed 


Table Factors detern from the graphical method on Bras ‘ ram ple 


Point il, 5 i 


0-100 
O-O85 


oe 


oom 


tS ts te te 


q 
BURN - MOUGE J 
NAKAJIMA 
/ 
Cr 
ty 
22-8 P.c.u. ‘Ib VOI 
ty 195°C 9 
‘ 1958/ 
| 
l 136-5 105-5 23-0 125 2-78 Los 1-78 
2 117-7 03-5 22-5 1024 2-77 240 
3 93-0 774 21-8 277 37-1 200 
$ 73-0 62-40 21-1 1017 2 33-3 20-0 33-0 3-09 
5 62-3 5340 20-8 2 30-9 2540 28-0) 3462 
7 12-8 B40 20-0 1-o1o 2 5-7 mo 6-32 
274) 109-7 147 2 3-5 wo 10-07 
20-4 22-8 19-5 1-006 2 2-2 160 is 14-78 
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that « is constant in a condenser. Therefore, Design procedure 


equation (5) should be written as Fig. 4 is an enthalpy-temperature chart of a 


h./k.’ ne (5’) Vapour-gas mixture other than air—water system 


with the saturation curve for modified enthalpy. 


Hence, 


Gy, di. (7, i, id { (6) 


where is modified enthalpy (9) and defined J ya, 


as 


i, 


Therefore, 


The 


equation (9), = 


inclination of operating line is given by f 


And the inclination of path line of eas condition 
ic. 4. Enthalpy and modified enthalpy 


Is temperature chart. 


di, | i } 
dt, 


(10>) 


And the 


inclination of tie line is { 


(11) 


lhe diff renee bn tween Fe and Uy is given from Fic. 5. Difference between modified enthalpy 


equations (7) and (7°). and enthalpy. 


i i (x —I)e,t = (a —1)(c, + c, H)t (23) 


Fig. 5 is a diagram of the relation between 


And this value is a function of ¢ and H. 


enthalpy and modified enthalpy. A group of 


inclined straight lines of parameter H is plotted 


Consideration of the change of values o h, and e 
£ 


by equation (23) and shows the relation between 
(7,.’ i,) and ¢ at constant H. It is convenient 


in the condenser 


Equations (12)(20) can be applied to this to take the ordinates and abscissae of Figs. 4 


case too. and 5 of the same scales. 


Hence, It is assumed that the condition of inlet gas, 


the enthalpy of outlet gas, the temperature of 


cooling water available and the water and gas 


flow rates are known as in the above section. 
Then, the operating line AZ and the point 


re presenting the condition of inlet gas, 


G 
hy (7 ) A tay) can be plotted on Fig. 4. 
= Plotting point A,’ (t,,, H,,) on Fig. 5, one can 
obtain — from the length of the vertical 
t; —t,; hy line, dy Ay. Plot points a and a’ on Fig. 4 so 
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} 4 tp 
a4 
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| A 4, tu / 
ie 
an ° 
h. (ty a; ) | 
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| 2 
9 H- | 
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SQ 
[ 
- 
di 
A= 2G, | 
J hy ti) 


T. N., 


that 4a A’ a’ 
represent ty) ard (i,,', t.,) respectively. 


A,’ Ay. The points a and a’ 


The values of cp/h,, 8 and h, can be obtained 
In the 


gas organic vapour mixture, however, the thermal 


just like the preceding section. case of a 
resistance of the condensate film predominates in 


hy and varies considerably in a condenser. 
Therefore, the value of h, should be calculated 
at And the 


can be calculated by equation 


Peau (P Pi) 
Pr and Sc. Thus, one 


each point. also value of xz 


assuming 
l. and 


using average 


calculates the inclination 


of the tie line a Bhocy hy, and constructs 
the tie line from a, and locates A” (i.’. t.) on 
the saturation curve for 7,’ ws. ¢.. 

Join point a’ with A Since the inclination 
of a’ is (i,, t equation (10°) 


)s 
of the 
point 4° is equal to (1 x) times th 
of a’ A Now, one « 


of gas condition from A 


path line at 
of that 
the path lin 
through B.. B 
ponds to the next point Bon the operating line. 


indicates that the inclination 
Lille 
‘an construct 


corres 


The ste ps ned above shou cl re peated 
until one reaches the final point 


Then. plot 3a i ) VS, and calculate the 


value of the integral in equation (21 Finally 


one can calculate the surface area 1 by equation 


(21°). 


The 


experiment 


com partson of re sults of Cah ilation and 


The procedure mentioned above is apphed to 
Run No. 14D of and Rorson’s experi- 
ment [1] 
At the gas inlet 


with air—benzern system 


0-698 


At the gas outlet 


O-5S4 
In averayve 


730 
Se 0-463 


0-730 


\ alue Ss of 


Hasnimoro 


Factors di termined from the graphical method on and Ropson rample 


3. 


Tabli 
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M. 


Point 


hn 


} 


ft* he 


{ 
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: 
| 
nS 
at 
Se * 
tz 
| ~ 
= a 
| 
= 
Hence 
202 


Now, one can make an enthaply and modified 


enthalpy-temperature chart and a diagram of A 
the relation between enthalpy and modified C 


enthalpy. 


The condition at the gas inlet and outlet are 


Cs 
2,754 P.c.u. lb d 
32-9°C 
ii Gr 
tie 29-7°( 
The inclination of operating line is 
Gp 
LC, 60 H 
680 
Gan O-l4 29 h 
Now, the operating line can bi fixed on enthalpy Mn 
temperature chart. 
At the gas inlet, the heat transfer coc flicient ho 
h, = 10-1 P.e.u. /ft® hr °C 
And 
H, 20-41b Ib 
{1 20-4)" 
(P24 20-4) (0-38) 
6°36 Ib °C 


Therefor 


Cr 6-36 L 
0-342 
hy (10-1) (1-845) M 


The value of hy is obtained at each point using n 
a conventional Nusselt equation for condensate P 
film conductance. p 
Point to point calculation results are shown in 
PBM 
lable 3. 
The values of integration in equation (21’) is t 
3-08, and 
x 
Gp 8 
(0-80) (0-342) (O-14 29) 1-11 
hy Ao 
The surface area the refore 


me P 

A = (1-11) (3-08) — 3-42 ft? ; 

Se 

Compared to 3-09ft® as calculated by the Re 


rigorous method of and Houcey, and 
3-3 ft® 


as found experimentally by Roperrson 


and Smrru, the agreement is good. B 
& 
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NOTATION 
area of transfer surface, ft® 
constant in equations (13) and (18) 
specific heat, P.c.u./Ib 
humid heat, P.c.u./Ib of dry gas °€ 
diameter of tube, ft 
flow rate of stagnant gas, lb /hr 
mass velocity of gas vapour mixture, lb/hr ft? 
mass velocity of stagnant gas, lb/hr ft? 


humidity, lb/lb of dry gas 


coefficient of heat transfer, P.c.u. /ft? hr 


fictitious coefficient of heat transfer when only 


stagnant gas flows, P.c.u./ft® hr “( 


combined conductances other than the gas film 
P.c.u. ft® hr 


enthalpy of gas vapour mixture, P.c.u./Ib of dry 


modified enthalpy of mixture, 


u 


gas 


vapour 
Ib of dry gas 


mass transfer coefficient, lb mole /ft? hr atm 


mass transfer coefficient, lb/ft? hr of dry 


gas) 


flow rate of cooling water, lb/hr 
molecular weight, lb /mole 
number of power 

total pressure in condenser, atm 
partial vapour pressure, atm 


average partial pressure of stagnant gas in gas 


film, atm 

temperature, 

constant defined by equation (4) 

variable defined by equation (20) 

latent heat of vaporization at 0°C, P.c.u. ‘Ib 


Viscosity lb/hr ft 


Prandtl] number 
Schmidt number 


Reynolds number 


Subscripts : 


stagnant gas 

gas—vapour mixture or gas film 
interface 

cooling water 


vapour 
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Abstract. -T'o obtain information on gas liquid flow undet operational processing conditions, 


im investigation was made of the flow of air water and air-—oil mixtures in horizontal smooth 


pipes with inner diameters ranging from 24 mm to 140 mm, and rough pipes with an inner diameter 
of 50 mm 


can be predicted with reasonable accuracy for all pipe diameters within the range investigated. 


i result of this, diagrams can be given from which the flow patterns for these mixtures 


For the pressure drop the Lockhart-Martinelli correlation appeared to be valid for plug, 


slug and froth flow at atmospheric pressure only ; for wave and mist-annular flow it was inad- 


equate under any conditions Separate new correlations are given for these cases those for 


plug, slug and froth tlow include the influence of gas density. 


The results of holdup measurements have been correlated « mpirk illy with the slip velocity 


between the phases. A special holdup meter has been deve loped enabling accurate measurements 


to be mack 


Résumé —0)n a étudié Pécoulement de meélanves 4 deux phases (air-cau et air -huile) en tubes 


rivontaux dun diametre interieur de 244 14) mm, et en tubes rugueux d'un diameétre 


lisses h 


intérieur de 50 mm, pour mieux comprendre le comportement des courants 4 deux phases dans les 


conditions de pratique 


Les résultats ont été utilisés pour tracer des diagrammes permettant d’établir avee un 


precision raisonnable les configurations de écoulement de ces mélanges dans des tubes d'un 


diametre intéricur quelconque, compris dans les limites de nos mesures 


Les valeurs de Lockhart-Martinelli pour la perte de charge ont paru étre valables pour les zones 


de circulation sous forme de tampons, de bouchons et de mousse a pression atmosphérique 


seulement ; pour les zones de circulation ondulatoire et annulaire avec pulvérisation, elles n’étaient 


applic ible sous aucune condition Pour ces cas de nouvelles valeurs sont données : celles qui se 


rapportent aux zones de circulation sous forme de tampons, de bouchons et de mousse, tiennent 


compte de Vinthuence de la densité du gaz 


Une relation entre les résultats des mesures du hold-up et la vitesse de glissement entre les 
| 


phases a ete ctablie dune facon empirique On a mis au point un jauge special qui permet ce 


mesurer le hold-up avec précision 


Zusammenfassung tm Kenntnis tiber die ZweiphasenstrOmung unter Betriebsbedingunyen 


zu erhalten, wurde die Stromung von Luft-Wasser- und Luft-6l-Mischungen in waagerechten 


glatten Rohren mit lichten Durchmessern zwischen 24 und 140 mm und rauhen Rohren mit 


Durchmessern ve untersucht 


lichten 


Als cm re bnis dieser Arbeit werden ben wus denen sich dic Stromunys 


form fiir diese Mischungen aussagen lisst und zwar mit hinreichender Genauigkeit fiir alk 


Rohrdurchmesser im untersuchten Bereich 


Fir den Druckabfall scheint die Beziehung von Lockhart — Martinelli fiir Kolben-, Tropf 


und SchaumstrOmung our bei atmosphirischem Druck giiltig zu sein Wellen- und 


Propfchen-Ring-Stromung war sie in allen Fallen nicht brauchbar Fiir diese StrOmungsarten 


werden neve Bezichungen mitgeteilt, Die Gleichungen fiir Kolben-, Tropf- und Schaumstrémung 


enthalten den Kinthuss der Gasdicht« 


Die Krvebnisse von Messungen des Fliissigkeitsgchaltes konnten empirisch mit = der 


Schlupfgeschwindigkeit zwischen den Phasen in Beziehung gesetzt werden besonderes 


Messyeritt fiir den Flissigkeitsgehalt wurde entwickelt und ergab ausreichende Messyenauigkeiten, 
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lL INTRODUCTION 


A numper of authors have studied gas-liquid 


flow in pipes of diameters small 
with the diameters used in technology. Quite 
recently some experimental data were given by 
Reip et al. [16] for air-water flow in 4 in. and 
6in. pipes, but only for a narrow range of air 


m OMmparison 


and water velocities. Baker [3] gives some field 
results for 4in. to 10in. pipes, the validity of 
which is restricted to the sper ial cases he studied. 

An investigation of gas liquid flow through 
horizontal large-diameter pipes was therefor 
made, to extend the region of pipe diameters 
covered. 

The influence of the following parameters was 


established : 


1. Pipe diameter. Smooth pipes of inner 
diameter D = 24, 50, 91 and 140 mm were 
used. The smaller pipe diameters were 
included not only for check measurements, 
but also for a series of measurements 
covering a wider range of gas and liquid 
mass velocities than reported in literature. 

2 Liquid properties, Water. spindle oil, 
and gas oil were used as liquids. 

3. Gas density over a narrow range. Air 


pressure ranged between 1 and 3 atm. 


° 


$. Pipe roughness. Various degrees of rough- 
ness wer investigated in pipes of 
D 50 mm. 


The flow pattern, the pressure drop, and the 
liquid holdup were studied, 


2. EXPERIMENTAL SETUP AND 
MEASURING PROCEDURE 

Figure 1 gives a flow scheme of the experi 
mental setup. The 25m horizontal test section 
was of transpare nt perspex pipes. 

The pressure drop was measured only at 
locations beyond an inflow length of 60 D. For 
the 24 and 50 mm pipes 3 successive measuring 
sections each 4m in length were used. and for 
the 91 and 140 mm pipes one measuring section 
of 8m. The measurements in rough pipes 
were made in 8 sections of 2m each, preceded by 
17 m of smooth perspex pipes ol the same dia- 
meter. Pressure taps were located in special 
flanges and the pressure was transmitted from the 
underside of these flanges to normal U-tube 
manometers through liquid-filled lines. Care was 
taken to keep air out of these lines. 

Differs nt entrance devices for the air and liquid 


were tested; the infiuence of change in inlet 


Flow scheme of 25 m test setup. 
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device was found to be negligible beyond 60 D, 
except in one case mentioned below. 

Four centrifugal pumps and a pressure vessel 
made a wide range of liquid flow rates possible 
(0-02 to 820 m®/h). High flow rates were measured 
with an orifice meter, the low rates by weighing. 
The liquid used was recirculated after de-acration 
in a channel 25 m long. 

Compressed air was used from a 6-5 atm. supply 
with a maximum capacity of 1800 kg/hr. The 
air flow rate was measured with orifice meters of 
different sizes. 

Liquid and air temperatures were measured in 
the separate supply lines and at the end of the 
test section. The oil temperature was regulated 
with a steam coil in the de-acration channel. 

Properties of the circulating liquid (see Table 1) 
such as density, viscosity, and surface tension 
were regularly checked. For the properties of 
the air the data for dry air could be taken sinc: 
its morsture content was nm vligibly low, 
liquid holdup a capacitiv: 


Three 


dises (spacing 9 mm) were mounted in the pipe 


For measuring the 


method was developed. stainless steel 


parallel to the direction of flow (see Fig. 


They form the electrodes of a capacitor, the 
capacitance of which depends on the liquid holdup 
in the pipe. The discs (same diameter as the pipe) 
are circular in shape except at the point wher 


they are mounted in the pipe wall. A linear 


Fic. 2. Holdup meter. 
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relationship between holdup and capacitance is 
thus obtained. This makes the method almost 
inde pendent of the configuration of the gas liquid 
mixture and makes it possible to determine mean 
values for the holdup in a fluctuating flow. Such 
holdup meters were used in the perspex pipes of 
D 50, 91 and 140 mm. 

The holdup meter was calibrated under different 

conditions : 

1. By filling a closed horizontal section of 4m 
length to a known liquid level, the holdup 
meter being fixed at its mid-point. This 
represents the case of stratified and wave 
flow flow patte rs as described below). 
\ linear rv lationship between capacitance 
and holdup was found. 

2. By 


position, 


placing a 3m section in a_ vertical 
again with the holdup meter at 
through 


its mid-point. Gas-—liquid flow 


this pipe may show different dispersions 
repre senting « ither plug, slug or froth flow. 
arrangement the holdup 


In a_ vertical 


can easily be determined by suddenly 
closing the air and oil supply. (Only in the 
case of a fine dispersion did the calibration 
curve differ somewhat from that found by 


method 1). 


The capacitive method was checked by means 
of a series of « xperiments in the 50 mm L.D. pipe 
with a radio-active tracer. For this purpose an 
oil-solublk radio-active tracer was homog ously 
mixed with the oil. In this case the radiation 
intensity of a part of the pipe depends linearly 
on the liquid holdup. The results showed that the 
differences in holdup between both methods were 
random, giving a standard deviation of 0-03. 
Only for froth flow was a significant difference 
found and this was probably due to a finer 
dispersion than obtained in the vertical calibra- 
tion. The capacitive method was corrected 
accordingly in that region. 

For each combination of pipe diameter and 
liquid a series of about 12 runs was done, each 
run being performed under eight or more different 
measuring conditions. The ranges of the super- 
ficial mass velocities of liquid and air are given in 


Table 2. Each series started with a run of one- 
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phase liquid flow in order to check the measur 
ing procedure. The friction factors thus obtained 
were in agreement within 5 per cent with the 
factors for smooth pipes given by Moopy [14}. 
For the rough pipes the roughness was found 
from the one-phase liquid friction factors by 
comparing them with the factors given by 


Moopy [14]. The roughness values thus obtaimed 
for the different cases were 


0-0030 and 0-0012, 

All the pressure drop data were corrected for a 
compre ssibility effect. When the pressure drop 
becomes appreciable th gas expansion vecelerates 
both the gas and liquid phase. An approximats 
impulse correction AP. based on mean velocities 
was used to obtain the friction pressure drop 
between sections 1 and 2 from the measured 


pressure drop. It was found from: 


AP. 


(1) 


The correction was small for the tests in the 91 
and 140 mm pipes, but for the 24 and 50mm 
pipes it amounted to 15 per cent of the total 


pressure drop in cases of high ve locity 


Table 1. Liquid properties at ope rating temperature 
of 28°C 


| 


Viscosity Density Surface tension 
Gas oil 2-9 815 27 


Table 2. Range of superficial mass velocities 


Pipe dia. D (m) (kg m*s) m, (ke s) 
05-165 1-5 6000 
O-140 | 0-2-55 > 5000 
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3. DISCUSSION O} RESULTS 
3.1 Flow patterns 


Visual observation shows that different flow 
patterns may occur with gas liqquid flow in hort- 
zontal pipes. In accordance with results obtained 
with small diameter pipes (Bercenin and 
[+], Atves [1], [9]) we observed 


the following flow patterns*) in our tests; 


a. Stratified flow. The liquid flows in the 
lower part of the pipe and the air over it 
with a smooth interface between the two 


phase 


b. Wave flow. Similar to stratified flow, 
except for a wavy interface, due to a 
velocity difference between the two phases, 
In wave mist flow the difference in 
velocity is so great, that a certain amount 


of atomization takes place. 


Plug flow. The gas moves in bubbles or 


plugs along th upperside of the pipe. 


d. Slug flow. Splashes or slugs of liquid 
occasionally pass through the pipe with a 
higher velocity than the bulk of the liquid. 
Pressure tluctuations are typical for this 


type of flow. 


‘ Mist-annular flow. The liquid is partly 
atomized in the gas phase and _ partly 
flowing in an annular film along the pipe 
wall. 

f. Froth flow. The gas is lisp rsed in fine 


bubbles through the liquid phase. 


Charts have been proposed (ISBIN et al. [7], 
Kosrertn [9] and Wurre and Huntineron [18]) 
to show the ranges of flow variables for the 
different flow patterns. They all use the gas and 
liquid velocities or combinations of these as 
variables and the different charts can be trans- 
formed into each other. The chart proposed by 
KosTertn [9] ts used for representing our results 
(see Fig. 3). 


*A film entitled “ Simultaneous flow of gas and liquid ” 
and showing these flow patterns for air-gas-oil flow in 
91 mm LD. pipes is in preparation and will be obtainabk 
on loan from the local Shell organization. 
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Here the relevant variables are the mixture 


velocity v,, and the gas percentage ¢,, with: 


= — and ¢, i < 100°), 
: 


Q, = gas flow rate by volume, 


@, = liquid flow rate by volume. 


Qur results with air spindle-oil and air-vas-oil 
flow in pipes of different diameters showed that 
the influence of liquid \ iscosity and pipe diameter 
on the location of the lines of transition between 
the various flow pattern regions is small. These 
transition lines are in fact transition zones with a 
breadth comparable to the shift of the transition 
lines due to change in pipe diameter or oil. Thus 
one chart (Fig. 3) can be given presenting all our 


Deviations from this 


results for air—oil flow. 


Fic. 4 (a). Flow regions in 91 mm LD. pipe. 


Fie. 4(b). Flow regions for high eg’s in 91 mm LD. pipe. 
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generalized chart may occur in cases where the 
gas density is different from that used in our 
experiments (1-2 — 3 kg/m!) especially near the 
transition lines. The same holds for cases of high 
liquid viscosities (> 50 

Our results for air—-water flow showed nearly 
the same transition lines (se¢ Fig. 4(a)) except for 
wave flow, which occupies a much larger region. 

Although the influences of pipe diameter and 
liquid properties on the various transition lines 
are small and have been neglected for the general 
diagram of Fig. 3 they can be discussed by a more 
detail dstudy of our results. Figures 4(a), 4(b), 5(a) 


and 5(b) rive some typical cases, 


a. The transition from stratified flow to plug 
and slug flow. This occurs at superficial 
liquid 


velocities v, Of about 0-2 m/sec 


ince pendent of superficial air velocity. For 


Stratified 


Fic. 5 (a). Flow regions for air—gas-oil flow. 


Fic. 5 (b). Flow regions for high cg’s for air—gas-oil flow, 
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the smaller diameters the transition occurs 
at somewhat lower velocities, because 
instabilities in the liquid surface can easily 
reach the wall at the top of the pipe 
(see Fig. 5(a), 24 mm pipe). 

The transition from wave mist flow to 
slug flow. For the largest pipe diameter 
(140 mm) this occurs at an air percentage 
of between 94-97 per cent. For the smaller 
diameters the waves reach the wall at the 
top of the pipe more easily ; therefore here 
again the transition occurs at higher air 
percentages. There is, moreover, the in- 
fluence of surface tension. For air-water 
flow the wave flow region is greater than 


for air—oil flow (Fig. 4(a)). 


The transition from wave tlow to wave 
mist flow. Atomization of the waves begins 


in all cases for v,, ~ 11 m, sec. 


The transition from wave mist flow to 
mist-annular flow. The liquid droplets 
formed by atomization of the liquid in 
the lower part of the pipe are deposited 
on the pipe wall and form a film of liquid. 
For ¢ « > BO m, sec, this film and the liquid 
layer in the lower part of the pip form one 
fully developed liquid) annulus. Thus 
generally the liquid annulus co-exists with 


a mist in the gas phase. We observed a 
tendency for the surface to curve upwards 
in the neighbourhood of the pipewall at 
Uy 15 msec. However, at these veloci- 
ties we always found atomization, even at 
the lowest oil flow rates O-O02 msec). 
This contradicts results of Krasyakova 


[11] for air-water flow in 30 mm pipes. 


The transition from mist-annular flow to 
slug flow. This occurs at air percentages of 
about 98 per cent. Again for the small 
pipe diameters the transition occurs at 


somewhat higher air percentages. 


The transition from plug flow to slug flow. 
This zone ts very wick no ral conclu- 


sion can, therefore, be given. 


The transition from plug flow and slug 
flow to froth flow. This occurs at 


v,, ~ Bm/sec. Here the pipe diameter 
shows the effect of shifting the transition 
to higher velocities for the greater pipe 
diameters. 


3.2. Pressure drop. 

The correlation generally used for two-phase 
pressure drop is the one given by LOCKHART and 
Martinetot [13], based on their results with 
air-liquid flow in 0-06-1 in. pipes at atmospheric 
pressure. They use the parameters >, and x, 


given by: 


AP,, AP, 
AP, AP, 


p? (2) 
pressure drop over pipe length 
AL for two-phase flow. 
pressure drop over pipe length 
AL for the gas-phase flowing 
alone through the pipe. 
pressure drop over pipe length 
AL for the liqquid phase flowing 
alone through the pipe. 


Our results with smooth pipes showed that for 
plug, slug, and froth flow under atmospheric 
conditions their correlation also holds for the 
larger pipe diameters. The standard deviation 
between our pressure drop data and their correla- 
tion was 13 per cent, while maximum deviations 
of + 30 per cent occurred. Our results showed, 
however, that in other cases the deviations were 
rather appreciable, for all pipe diameters. Closer 
examination of the relevant data indicated that 
the Lockhart-Martinelli correlation does not hold : 


for plug, slug and froth flow if gas densities 
differ from that of air at atmospheric 
pressure. 

b. for stratified flow and wave flow. 


c. for mist -annular flow. 


This probably explains the discrepancies as found 
by various authors when using this correlation 
(Berceun and Gaz [4], Cuknowern and 
Martin [5], Baker [3] and Rerp et al. [16)). 

For the results in the 50mm rough pipes 
comparison with the Lockhart-Martinelli correla- 
tion shows discrepancies when for AP, and AP, 
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in their parameter x the values for rough pipes 
are used. By using, however, in x the values for 
smooth pipes at the same conditions and in ®, 
the value for rough pipes reasonable agreement 
is again obtained for plug, slug, and froth flow 
at atmospheric pressure. This is also confirmed 
by CutsHo_m and Larrp [6] in a recent paper on 
They tind, 


larger values of 


air-water flow in rough 1 in. 
that at 


« Dan increase in pipe roughness has less influence 


pipes. 
however, in addition 
with two-phase pressure drop than with one 
phase pressure drop. This did not appear from 
our measurements. 

We developed new correlations for pressure 
drop for those cases where the existing ones 


appeared to be inadequate, as follows. 


a. Plug, slug and froth flow ; 


de nsily 


influence of air 


For plug, slug, and froth flow it was found 
that though the Lockhart-Martinelli correlation 
holds at atmospheric pressure, the values it gives 
for the pressure drop at higher air densities were 
too high. Consequently we tried by using dimen- 
sionless groups to develop a new correlation for 
these flow patterns, which would, in addition, 
reflect the influence of gas density. 

The two phase pressure drop is considered to be 
a function of the following quantities : pipe length 
AL, pipe diameter D, superticial liquid mass 
velocity m,, superticial 


gas 


mass velocity mn. 


liquid density p,, gas density p,. liquid dvnamic 
Py Y Pg | 


viscosity 7 ras dynamic viscosity »,, surface 


tension o, and acceleration of vravity g. 


Since in the Lockhart-Martinelli correlation the 


influence of surface tension and gravity is 
neglected without loss of accuracy, it is felt that 
these quantities can be ignored for the flow 
patterns under consideration within the range of 
variables investigated. The following six dimen- 


sionless groups can then be formed : 


AP,, AL m, ”) 
where m, = m, + m,is the mass velocity of total 

gas and liquid flow. 


Since the pressure drop is directly proportional 


Gas-—liquid flow in horizontal pipes 


to the pipe length (with the other variables 


constant) we have the following relationship : 


t g 
Here : AP,, 
1 AL m? (4) 
2D 
and 
m,D 
(5) 


analogous to the friction factor and Reynolds 
Tre choice of these 
the 


sideration of the flow of an evaporating medium 


number in one-phase flow. 


groups based on m, stemmed from con- 


in a pipe where the mass flow of liquid and gas 
separately is continually changing, but where the 
total mass flow is constant. 

The intluence of the gas density could be found 
by analysing our data for the 24 and 50 mm pipes 
taken over 3 separate sections of 4m _ length 
each. Due to the pressure drop a difference in 
P, p, existed between the 3 sections for each test 
condition the other dimensionless groups remain- 
ing constant. We assumed a power law to exist 


between A 


tp and Py/ giving : 


x 
P 


in (6) 
(Pg/ Pt) 
tpn A,, for the fixed value 


Here (A,.).. is the value of A 
1-38 x 10°. Figure 6 gives some 


ol 


p 
p D 


Fic. 6. Influence of air density on pressure drop with air 


gas-oil flow in 24 mm pipe. 
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re sults showing the 


influence of density on 
Here AP pressure drop over pipe length AL 


\ Here a friction factor A, defined below and 

which would occur if the whol 
independent of the was density, ts used to bring 

zy mass of gas and liquid wer 

t he eurves together m one heure, nt 

flowing in the liquid phase through 
soa function of other ivy was 

the pip 
found to be a function of m, m, a It is 
easily seen that in the case of all yuid flow Furthermore s given by (7). 
(wm. 0), x equals zero, in the case ol all was 
flow (mm, m 1), x equals unity For inter 
mediate values of m, m, our result ive for 


the « Trapani al correlation 


The ol the other ould bv 4 ey 4 6 
analysed separately from that of 
because hia. 7 New pressure drop correlation for plug, slug and 


J Re, 


The standard deviation for our 492 test data 


Our data showed that curves could be plotted 


of versus Re and m, (see I 7 thre obtained in 24 50, 91 and 140mm smooth 
ira net of the viscosity ratio ould be pipes with air-water, air-gas oil and air-spindle- 
neglected in the range of viscosities vestigated tl flow compared with this new correlation 
from 60 to 1500), These cur can be Was 6 percent form, n O-OOT and 14 per cent 

lor O-O0T m, 0-05; in the second range 


pre nted by 


most data were for slug tlow type where the 


\ pressure fluctuates strongly and an accurate 

230 | ‘| determination of th pressure drop is impossible. 

‘ The influence of gas density was investigated 

for O05 and Re SOOO, over a range of Il 3-3ke the factor 
we 0-O0138 (p, py) then varies from 1-3. In Fig. 6 
this influence on A, Ay, is shown for some test 


Here. ly oS the friction factor wi! ich would data at various values of mom, it is seen that 


apply if the whole mass of gas and hund wer extrapolation to much higher densities is not 
flowing in the hquid phase. Including the im permissible, since the curves would intersect the 
fluence of density, this leads to the general / At l. Further research at high vas 


densities would bn required, 


correlation : 


For 


0-00138 


(10) 


= 
x l tor VOI 
1958/ 
froth fthow 
A \ 230 | | m, 
my, 
| x 
or AP. 230 *) 1o-00138 | Re, 
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The data obtained in the 50mm rough pipes xX 1. The pressure drop shows a sharp decrease 


with air-gas-oil flow were analysed in the same there. For 1 X < 2 the conditions just cover 
way as those of the smooth pipes. Pipe rough the transition from wave to slug flow. By using 
nesses studied were « D 0-030 and 0-019 different inflow devices it was possible be obtain 
{cast tron pipe) and « D 0-0030 and 0-0612  cither slug or wave flow. th other conditions 
(ground cast tron pipe). The same correlation being the same. In this way the overlap in the 


formula (10) as developed for smooth pipes holds two parts of this curve was found: this clearly 


with the same restrictions for all values of « D. showing the influence of flow pattern on pressure 
if for Ay the value is taken corresponding to thy drop. 
€, D of the rough pipe. The results for the pressure drop with wave 


llow were analysed s« parately. It was found that 


b. Stratified flow and wave flox AP 
the dimensionless vroup , could be corre- 
For stratified tlow and wave tlow the pressure I 1 AL m? ‘ 
drop was lower than predicted by the Lockhart » 
Martinelli correlation, as can be seen for wave lated with m m, (see Fig. 9) 


flow in Fig. 8. For stratified flow only afew data giving 


are available, indicating the same trend as for \p aie 


ip 
9 wave flow. For various superticial liquid mass = ( a (11 
| ) 
AL m, 
velocities different curves are found when th. ; 
2 D op 
Lockhart-Martinelli parameters are used, This ré 


Fic. 8 Data for air—gas-oil flow in 140 mm pipe correlated 


shows that these parameters are inadequate for 


correlating the pressure drop in the region of wave 


flow. This behaviour coincides with the transition 


to a different flow pattern as can also be seen in hic. 9% Pressure drop correlation for wave-flow results 


Fig. 8. The curve for the constant supertie ial in 14) mm tube with air—gas-oil flow. 
liquid velocity vy = 031m see consists of 2 
parts. For x > 2 the air mass velocity is small 


and slug flow exists, the results being in agreement Here C was found to depend slightly on pipe 


with the Lockhart-Martinelli correlation. With diameter and liquid viscosity, see Table 3. The 


increasing air velocity the flow pattern changes air density could not be varied during these 


to wave tlow this being the part corresponding to tests and amounted to 1-2 kg /m*. 
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Table 3. Influence of pipe diameter and liquid 


viscosity on C 


Diameter lir -gas-oil fir spindle-oil 


m moe 


mool m mor 


mart 


In this case again a marked influence of pips 
roughness on the two phase pressure drop was 


found. The experiments with air -gas-oil flow 
in the rough 50 mm pipes showed that the same 
correlation formula (11) could be used, with the 
constant C, depending on the pipe roughness as 
given in Table 4. 

Table 


Influence of pipe roughness on C, for 


aur gas oul flox in SO mm pipes 


Pipe roughness « D 


moor 


mor 


mosz 


Mist-annular flow 


For this flow pattern measurements were made 
the D 24. 50 


and 91 mm for 


the 


again 


only in 
and in 50 mm for 
air Here the Lockhart- 
Martinelli parameters were found to be inadequate 
Fig. 10. This 


contradicts Levy's [12] theoretical prediction. 


air—gas-oil pipes 


spindl -oil tlow. 


for the correlation of our data, se« 


10. Data for mist-annular flow in 50 mm_ pipes 


correlated with Lockhart-Martinelli parameters. 
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For mist-annular tlow the gas is the continuous 


phase. Therefore a new dimensionless group 


was based on the gas flow alone: 
AP, 
1 AL 
2D op, 


(12) 


The results are thus correlated with a Re vnolds 
Fig. 11. 


number based on the gas flow alone. see 


Experimental pressure drop curves for mist 


annular thow. 


It was found that for liquid flow rates which are 
the 
independent of this rate. 
oil flow Ar, then be 


alone 


not too low, friction factor Avy becomes 


kor the case of air-gas 
could correlated with m, 
giving 


A, a ‘ (13) 


Pp 
holds 


influence of 


This result for 

No the 
was found over the experi- 
1-2-3 kg 


For lower superficial liquid mass velocities A,, 


a 0-12 kg* m ? see 4. 
30 m, 200 ko m®? sec. 
air density on Avy 
of air densities of 


mental range 


decreases with 
For rough 50 mm pipes with air-gas-oil flow it was 


found that in the case of not too great roughness 


0-0012 and 0-0030, or ¢ of 0-06 and 0-15 mm} 


m kg/m se 
- air / go 4mm 
——— Air/gas-ai SO mm pipe 
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the pressure drop was nearly equal to that for 
mist-annular flow in smooth pipes. For great 


€ 
roughness | 0-030 and O-O19 or « 15 and 


0-95 mm | the pressure drop was about 10 per 


cent higher than for air flowing alone. 


3.3 Holdup 


In most of our experiments with air-oil flow 
the holdup was measured. This was done to 
obtain a more accurate compressibility correc- 
tion for the pressure drop and more information 
about the gas-liquid flow. In the literature data 
for the holdup are reported by Lockuarr and 
MARTINELLI [13], Jounson and Anovu Sane [8] 
Auves [1] and Koziov [10]. These results wer 
all obtained with the quick closing valve 
technique, and they mutually show differences 
when correlated with the Lockhart-Martinelli 
parameter x as is done by most authors. Our 
results for the 50, 91 and 140 mm pipes obtained 
with the capacitance method also showed devia- 
tions from this correlation at various liquid 
mass velocities, indicating that x, which is only 
determined by the ratio Ugg inadequate 


and that both v have to be accounted 


s 


; and vy 
for. 

For vertical two-phase flow it is found by 
et al. [17], Moore and [15] and by 
Baitey et al. [2] that the holdup is mainly a 
function of the slip velocity (v') between both 


phases, where : 


“sg 5 (14) 


The use of this type of correlation for our hori- 
zontal test results was found to be successful, 


and gave for the gas holdup R, with air—oil flow: 


(15) 


Figure 12. gives 


Here A 


some test data. 


0-60 (m ‘sec) 
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Fic. 12. Gas holdup for air-spindle-oil flow in 140 mm pipe. 


No significant effect of pipe diameter or liquid 
viscosity was found in our velocity ranges. 
The gas density influence was not investigated. 
The standard deviation in R, for our 444 data 
compared with this correlation was 0-04 (absolute 
value). 

CONCLUSIONS 
1. The flow patterns occuring with gas-oil 
flow can be predicted from one diagram, 
the influence of pipe diameter and liquid 
viscosity being small. For gas—water flow 
the same diagram with an enlarged wave 


flow region cap be used. 


The Lockhart-Martinelli correlation has 3 

limitations ; it cannot be used: 

a. for plug, slug and froth flow, if gas den 
sities differ from that of air at atmos- 
pheric pressure, 

b. in the stratified and wave flow regions, 


c. in the mist-annular flow region. 


3. For plug, slug and froth flow a new experi- 


mental correlation which includes the 


influence of gas density and pipe roughness 


GIVES : 


0-84 x 


0-00138 
Pg 


AP,, |) 280 


with: 


F 
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rs =z. x0 
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correlation gives: 


For wave flow a new 


AP,, 
1 AL 


Pe 


for 


Here C depends mainly on pipe diameter 


and pipe roughness (Tables 3 and 4). 


For mist-annular flow AP,, is independent 
of the if 
than 30 kg m*sec; for air gas-oil 


liquid flow rate this is greater 


flow a 


new correlation then gives: 

1 AL m2 

dD 


0-12 (m_) 4. 
Pe 

The gas holdup correlates with the slip 


velocity between both phases : 


0-60 i 
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acceleration of gravity (m sec?) 
length of measuring section (m) 
superticial mass velocity (kg m® sec) 
two-phase drop over measuring 
section (N 


drop over measuring section if the 


pressure 
m*) 
pressure 
liquid were flowing alone m*) 

drop over measuring section which 
if the of 


liquid were flowing in the liquid phase (N /m*) 


pressure 
and 


would occur whole 


drop over measuring section if the 


over measuring section due to 


pre ssure 


vas were flowing alone mim) 


drop 


pore ssure 


the expansion of the iN 


flow rate by volume (m° sec) 


holdup 

Revnolds number 
Revnolds number detined by equation (5) 
velocity (m 


mixture 
superticial velocity sec) 

slip velocity (m ‘sec) 
Lockhart-Martinelli: parameter 
empirical exponent in equation (6) 
equivalent sand grain roughness (m) 


dvnamic viseosity (Ns m*) 


friction factor for two-phase flow defined by 


equation (4) 


\ for flow at 


tp two-phase 


friction factor 


Choo ros 


friction factor for mist-annular flow defined 


by equation (12) 

friction factor which would apply if the whole 
mass of gas and liquid were flowing in the 
liquid phase 

(m? seu) 


kinematic viscosity 
fluid density (kg/m 
surface tension (N 
Lockhart-Martinelli parameter 


l 


) 
mi) 
Subscripts liquid 


t 


pas 


total of gas and liquid 
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Mass-transfer from single solid spheres—IlI 


Transfer in free convection 


H. Garner and R. B. Kees 


Chemical Engineering Department, Birmingham University 
(Received 13 January 1958) 


Abstract \ theoretical analvsis of the problem of free convection shows that for Schmidt 


numbers exceeding LOO 


(N Ra)" laminar convection 
x turbulent convection 
od \p 
where Ra the mass-transfer Ravieigh number, defined by N Ra . 
Preliminary experimental observations, emploving the dissolution of compressed spheres of 

adipic and benzoic acids between gin. and jin. in diameter, confirm the above relationship 
for laminar convection. There ts a small difference between the mean curve through the results 
and the limiting value of Ny at low Reynolds numbers found previously by the authors (Part I*) 


It seems that turbulence sets in over the lower hemisphere at a critical Rayleigh number of 


3.5 10*, where the characteristic dimension is based upon the length measured along the 


surface from the upper pol 


Résumé Une analvse theortque du probleme de convection libre montre «jue pour des indices che 
Schmidt superieurs a 100 

N gn (Np) ‘ convection haniinaire 

Nig, (Np)! convection turbulent: 

gd" Ap 
ou Np, represente Pindice de transfert massique de Rayleigh detini par N Ra 
Des observations experimentales preliminaires utilisant la dissolution de sphéres comprimces 

des acides benzoiques et adiy ws dun diamétre compris entre gin. et }in., confirmérent la 
relation ci-dessus. Il vy a une faible différence entre la courbe moyenne passant par les résultats 
et la valeur limite de Ng pour les nombres de Reynolds faibles, trouvee précédemment par les 
auteurs (Part I* Il semble que la turbulence commence sur Themisphére inférieure pour 
un nombre critique de Rayleigh de 3-5 10, of la dimension caractéristique s‘établit dapres 


la longueur mesuree le long ce surface a partir du pole superieur 


Zusammenfassung —\us ciner theoretischen Untersuchung tiber die freie Konvektion ergeben 


sich folgende Proportionalitaten fiir Schmidtzahlen tiber 100 


Vg, ~ (CN — turbulente Konvektion 
Ap 
wortn die Rayleighzahl der Stoffiibertragung bedeutet, definiert durch N’,, 


Frihere Versuche itiber dic \ullosung von komprimierten’ Kugeln aus Adipin — und 
Benzoesiure von 0.5 bis 19 mm Durchmesser haben die obigen Bezichungen im laminaren Bereich 


bestatigt. Zwischen der mittleren, durch die Messpunkte gelegten Kurve und dem Grenzwert 


von Ne bei kleinen Reyvnoldszahlen, wie er friiher von den Autoren gefunden wurde (Teil 1*) 
besteht cin geringer Unterschied ks scheint, dass an der unteren Halbkugel bei einer 


kritischen Ravleighzahl von 3.5 1° Turbulenz cinsetzt, wenn als charakteristische Liinge die 


Bogenlinge vom oberen Pol cingesetzt wird 


*Chem. Engng. Sci. 1958 9 119. 
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INTRODUCTION 


In Part I of this paper it was shown that forced 


convection could be assumed to predominate, 
when the Rey nolds number had reached a certain 
critical value defined by 


> 0-4 (N,,)' No, 


In the case of } in. diameter benzoic acid spheres 


dissolving in water at 30-0°C. free-convective 


effects were still discernabk up to Reynolds 
number of 750 and at Reynolds numbers below 
10 the variation of the Sherwood number with 
Reynolds number is small. Thus a study of free 
convection should hy lp toevaluate more accurately 
mass-transfer at low flows. 

Apart from the recent work of Marurers 
et al, [8], there seem to be no data on mass transfer 
from spheres in free convection, (Problems of 
transfer from very small spherules in still media 
are those of molecular diffusion). Thus a preli- 
mipary investigation was undertaken to find out 


the factors of significance. 


Turory 


\ dimensional analysis of the problem vields : 
Ad f Ap | 


p 


(1) 


where the dimensionless groups are respectively 


the Sherwood number. the modified Grashof 
number and the Schmidt number. It is con- 
venient to rewrite this relationship in the form : 
Nig Ny. ) Nia N, (2) 


Sec’ 


since for practical purposes the Sherwood number 
is dependent only upon the modified Rayleigh 
number when the Schmidt number exceeds 100. 

It was Scumipr and Beckmann [11] who first 
suggested that the problem could be solved by 
using boundary-layer theory, in which all 
velocity and concentration changes are confined 
to a thin zone near the surface. There are two 
principal differences between this laver and the 
one associated with forced convection, Firstly 
force F, 


Secondly, the radial con- 


extraneous 


we have to consider an 
where F 
centration and velocity profiles will each show a 


Ap p- 
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Transfer in free convection 


maximum within the layer and the parameters 
at the boundaries will be zero. By comparing 


buoyancy forces with viscous forces over a slab, 


14 
On | | 
~ 


L 


we find 


where a» is the distance from the leading edge 
and L the total length. 
greater than unity and if 5,, is to be of an order 
less than L then 


Since vo L is never 


N 104, 


Ra 
From equations (1), (2) and (9) in (Part I) 
we can get the boundary -layer equations for free 


convection : 


(ru) d (rv) 
0, 
ou Ou 
au ou Ap 
i g sine > (3) 
oY p 
Ar oc 


where ¢ is the angle between the normal to the 
surface and the vertical. 
y 0 to y h (h 


By integrating from 
Oo), we get 


h h 
Iidf, du 
ru*dy £ sin y dy 
r da oY y=0 
0 0 
and 
h 
l d dub 
r uy dy Y | (4) 
r da y=0 
0 


in which (c 
pair of 


— ¢,). The equivalent 


equations for heat-transfer has been 
solved by Merk and Prins [9] by substituting 


the polynomials in »=y/5 for u and 


—T,)/(T, —T7,). To a first approxi- 
mation these are ; 
= (I and u = u, — 


A further approximation is given. Their mathe- 
matical argument is not given here. The method 
of solution is similar except that @ and Np, 


are replaced everywhere by ¢% and Ny. Their 
final result may be written as : 

4 y 

Ns 0-559 Np, * N sc > 1. (5) 
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Graphs showing the variations of local transfer 
rates and variation of maximum tangential 
velocity in the boundary-layer are given in 


Figs. 1 and 2. 


Local heat transfer 


Maximum veloeitv bounda 


The method breaks down towar 
pole, as the boundary-laver thi 
in this region and is infinite at Iso Furth 
more it is assumed that there is dilference 
between the momentum and mass-transfer tiln 
Theoretically this is only true wher 
at a Schmidt number of unity 
solutions obtained by Merk 
for a flat plate agree quite well 
more rigorous calculations [12 
of heat-transfer, the theoretical expression 
predicts Nusselt numbers about 7 er cent 
lower than those obtained 
Jaxon and Linke [7 

If the Ravleigh number is large enou 
transition to turbulence boundary 
takes place (10° N pe 


of Grirrrrus and Davies [5 


surements 

protil 
My 

Differentiation of tl 


in this case. cx pression 


vives an intinite velocity gradient at the surface 
and to give a tinite value to the skin-friction, 


a laminar sub-layver can be postulated. Zisnen’s 


data [14] show that this sub-layer extends up to 


° 12, whe re y* a dimensionless radial 


distance, based upon the friction velocity u* ; 


Between this laminar layer and the fully turbulent 
core is a quasi-turbulent “ buffer” zone and for 
calculation purposes an effective laminar film 
can be assumed up to y* 15. Using this 
concept Bayiey [1] has obtained an expression 
for free-convective heat-transfer in two-dimen 
sional flow, which can be written in the form for 


mass-transter 


Nin (Nw) 110; 10° \, (6) 


It is seen that this relationship iS independent 
of the characteristic dimension of the body. 
equation 6) should be applicable im the 
case of the dissolution of spheres in turbulent 


tree convectior 


EXPERIMEN’ 


which was described in Part 1 

of * Analar rack acipu and 

four sizes from gin. 

ere compre ssed in the 
hand-press, described in Part L As before, 
photo raph plates wer exposed at suitabl 
time intervals during the dissolution period. 
After development these were subsequently 
projected in a diascope and superimposed to 
obtain the diminution of the sphere. From the 
radial diminution local mass-transfer coetlicients 
and overall transfer coetlicients can be estimated. 
The following data are required pellet 
density, the solubility and diffusivity of the 
solute in the solvent and the density and viscosity 
of the saturated solutions. Pellet densities were 
found by Archimedes’ principle. Solubilities and 
diffusivities were calculated from information 
in the literature [2, 3, 13). Viscosities of the 
saturated solutions were assumed to be the 


same as that ol pure wate! at the same tempera- 


ture. Densities of the saturated solutions were 


| 
™ 
* 
T 
y* and u* 
4} 
> \ 
i 
\ 
4 4 
Fic. 1. rates 
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fi 
/ 
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4 
4 
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measured in a 10 ml pyenometer, which enabled 
10.000, 


densities to be measured to 1 part in 


The mean results of three determinations at 


each temperature are given underneath. 


Saturation densities (g ml) 


Table 1. 


Distilled water Benzoic acid in 


idipic acid in 


wales 


mater 


O-09056 


25 


iff rene 


accuracy of tha 


The average 


density between a saturated solution of adipic 


acid and water is 2 per cent and that between 


a saturated solution of benzoic and water is 


The latter inaccuracy is not as 


25 per cent, 


serious as it would appear on first sight, for 


iS proportional to 
This 


a factor 


the mass-transfer coeflicient 
the fourth root of the 


reduces the significance of the error by 


density difference. 


of about four. 


5 
hia. 3. Free convection overall transfer rates. 
Adipic acid upper pole support 
acid lower ole support 
Benzoic acid upper pole support 
Benzoie acid lower pole support 
= 
= 
8, 
i 
Fic. 4. Local mass-transfer rates in laminar free conves tion 


Nara 
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Transfer in free convection 


The results are presented graphically in Figs. 3, 
fand 5. The overall transfer results are applied 


to the region defined by 


Ng, 10° < 588 and 798 Ng < 2,200, 


Lawal mass-transfer rates in turbulent free 
convection 


Naw 2-32 x 10°. 


Dis 


l.xperiments were performed with the support 
both at the upper and lower pole of the sphere. 
From Fig. 3 there is seen to be no significant 
difference between the effects of the two methods 


Unlike 


convection there is no back-flow re gion, where a 


of support. forced convection, in free 


support can hinder or promote separation of 
the forward-flow. The effect of the support in 
free-convective flow is to add vorticity into the 
effect Is 


vorticity 1s 


boundary layer and the resultant 


independent of the point where 
added, 
Theory suggests that there is a linear relation- 


Ny, and Np, 


assumption, the results for the overall transfer 


ship between Based on this 


rate were correlated by the principle of least 


squares. The correlation coefficient of the line 


so obtained was 0.59. As this represents a low 


degree of correlation, individual points were 


re-examined, It was noticed that all spheres 


with a modified Rayleigh number exceeding 
149 10° showed a pronounced minimum in the 


local mass-transfer rate as some position other 
than the rear pole. Accordingly it was decided 
to neglect these in the general correlation. The 
0-68 and the 


is 


new correlation coellicient was 


regression line so obtained 


Ns, = 0-58, (Np,)'/* + 23. (7) 


15 
20 | 2 
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The this This 


may be compared to the t-value for a probability 


i-value of correlation is 3-20. 
of 1 per cent that this had occurred from chance ; 
the value so obtained is 2-92. We thus concluds 
that the 


noticed 


will be 
fall 


remainder It is 


correlation significant It 
that 


higher 


there are four which 
than the 


correlation 


pomnts 
somewhat 
these which render the coctherent 
so low. This could not be ascribed to any physical 
cause and the points remain anomalous, « spec ially 
since these are consistent with the limiting 
value of the 
Paper 1. 
all the 


the same order as the maximum error 


Sherwood numb r reported in 


results is 14-2 per cent 
ittributabk 
to the CX} rimental technique ( 16 per cent 


To reduce these unavoidable errors. greater 
magnification of the processed plate ss required 
needed, in which the saturation 
different 

In connexion with the 


high 


somewhat 


and solvents are 


densities are much from the solvent 


density. latter condition. 
solubility generally 


Mi re 


however, solutes of 
tend to 


accurate work is contemplated, 


dissolve irregularly, 


that the 


of th 


It is seen from equation (7 
of the 
number agrees well with Merk 
That the 


ri 


mental value coetlicient Ravleigh 


and Prins’ [9 


theoretical value. condition : 


Lim Ne 
0 


Ra 


is not satistied is not surprising, since the lowest 


Rayleigh number used in this work is 3-62 10°, 


Data are given by Maruers and ef al. [8] on the 


sublimation of benzene- and naphtl alene-coated 


sphe res at reduced pressures, when thermal 


The therein 


gradients are insignificant. results 


Fic. 6. Isothermal mass-transfer in free convection 
Maruers ef 


© Present work 


However, the mean probable error of 


which is of 


Garner and R, 


covered a much lower Rayleigh number range 
than those of the present work and thus form a 
useful comparison (Fig. 6). The two sets of data 
scem to be consistant. 


Phe 


seatter 


results for the local mass-transfer rate 


more than those for the overall rate. 


latter 


ereat 


since the transfer 


rate at the 


approximates to the 


circle (¢d 00°) and superim- 


position errors here are at a minimum. To 


overcome this difficulty local mass-transfer rates 
(for Ny, < 1-49 
of the product 


10°) were expressed in terms 
Ne, No 


should be 


which according to 


theoretical reasoning almost inde- 


igh 


Based on this assumption the 


pendent of number, (see Theory). 


mean probable 
rates at successive angles, 


errors of transfer 


60° apart, beginning at the upper pole, are : 
15.9 yp r cent. 16-7 ye r cent, 17-3 yp» r cent and 
These 


overall transfer rate 


24-58 per cent. of the same order 


as that of the 


errors ar 
and the above 
assumption is considered justifiable. 

The 


upper pol 


maximum mass-transfer rate occurs at the 
and progressively decreases around 
110 
until it is a minimum at the rear pole. Over the 
140 ) 


agrees 


the surface with a point of intflexion at é 


majority of the sphere surface (up to ¢ 
though not the 
well with that predicted from laminar boundary- 


the shape, magnitude, 


laver theory. There is a very rapid increase in 


the boundary layer thickness over the rear 


hemisphere until it is infinitely thick at the rear 


This thickening gives rise to the formation 
tail, 
rendered visible by the density striac. 


pole. 
sometimes 
Photo 
grahs of the dissolution of sugar-coated spheres 


illustrate the 


of characteristic which is 


in Fig. 7 higher transfer rate over 
the upper hemisphere and show that the boundary 
laver is infinitely thick at the rear pole, 

This 
the surface can be explained by considering the 
This 


concentration, 


distribution of mass-transfer rate over 


nature of the free-convective flow. flow 


is set up due to a difference in 
which gives rise to a density difference between 
the sphere and the surrounding fluid. It is thus 
expected that there is a general inflow, roughly 
the From 
considerations of continuity the heavier fluid 
Solute 


radially towards upper surface. 


must fall downwards around the surface. 
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-free thuid will then stream upwards to take the 
place of the inflowing fluid, giving a general 
circulation as shown in Fig. 11 (Part I), which 
Hence the 


boundary-layer is finite at the upper pole and 


has no effect at infinity upstream. 


progressively thickens as the fluid diverges from 
the sphere surface until it is infinite at the rear 
pole. Since the driving force for mass-transfer 
remains the same, this means that the concentra- 
tion gradients become less, whence the mass- 
transfer rate decreases. Boundary-layer theory 
predicts that the tangential velocities increase 
with increasing 4. If this condition holds when 
boundary-layer theory fails, then the increasing 
rate at which material can be convected from 
the surface would tend to counterbalance the 
thickening of the boundary-layer. This would 
explain the point of inflexion. 

Table 2 summarizes the data for those points, 
which were neglected in the general correlation. 
Rayleigh 

length, 


The third column (N,,). gives the 


number based upon a_ characteristic 
measured from the upper pole to the ring of 


minimum transfer. 


Overall transfer rates at high Rayleigh 


numbers 


Table 2. 


10 d, (Nya), » 108 N sp, 
149, 154 83 
201 146) 64 
232 102 206 
316 118} 351 80 
118 76 
550 100 304 131 


If (Ny). is assumed to be 
Rayleigh number, the mean probable error of 
these values is 13-9 per cent, which is of the 
sane order as that for the overall transfer results. 
constant 


Thus either (N,,). isa (approxi- 


mately 3.5 108) or else the variation thereof 


is masked by the scatter. This value is the same 
as that which accompanies the onset of turbulence 
in the free-convective flow around a horizontal 
Eckuert and Soruncen [4] have 


cylinder [6]. 
studied the flow 


around a_ heated, vertical 


Mass-transfer from single solid spheres 


independent of 


Il. Transfer in free convection 


aluminium plate in air with an interferometer 


and found that at small thermal potentials, 


transition to turbulence took place at 


4 10%, 


correlate heat-transfer from spheres by a linear 


Scumipr [10] was able to 


plot of Ny, against N,,'* over the 
10° Ny, < 10, Thus in the 


turbulence may be present over the rear hemis- 


range 
present case 
phere. It would be expected that turbulence 
would set in in this manner, as maximum tangen- 
tial velocities occur here. 

A typical distribution of local mass-transfer 
rates ts given in Fig. 5. The pronounced minimum 
eV idently corresponds to the onset of turbulence. 
The eddies within the boundary-layer enable 
material to be convected more rapidly, whence 
the increase in mass-transfer rate below this 
position. 

CONCLUSIONS 


Mass-transfer in laminar free convection can 
be correlated by a linear relationship between 
gd3 Ap 


and *, where N ita Maximum 


transfer takes place over the upper hemisphere. 
At a Ravleigh 


(where the characteristic dimension is given by 


critical number of 3-5 108 


¢éd 2), it appears that turbulence sets in. This 
causes an Increase in the mass-transfer rate over 


the lower hemisphere. 


NOTATION 


concentration 
d = sphere diameter L 
Y — diffusivity 
g — gravitational acceleration LT? 
mass- transfer coeflicient LT 
distance along the surface of the 
sphere from the upper pole L 
Ne, — Grashof number 
Nyy Nusselt number 
Nig = Rayleigh number 
= critical Rayleigh number (defined in 
Text) 
Nye = Reynolds number 
Ns, Schmidt number 
Sherwood number 
r radial distance L 
ro distance from axis of rotation L 
uit velocities in and y directions LT 
T — temperature 0 


u* friction velocity (defined in text) 
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distances mutually at right-angles dimensionless concent ration 


dimensionless distance, normal to surface (defined in text) 
(defined in text) angle, from upper pole 
boundary-laver thickness Subseripts 
angle (defined in text) at surface 
at of be 3 
dimensionless distance, normal to surface ut oundary-laver 
Operators 
(defined in text) 


partial differential 
shear viscosity L d total differential 


kinematic Viscosity 127 ! \ Anite-differenc 
fluid density ML” Superseribed dashes indicate the dimensionless groups 


skin-friction intensity VWLiT? are modified for mass-transfer 
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On the action of surface active agents on mass transfer in liquid=liquid 
extraction 


G. and S. G. Tersesen 


The Technical University of Norway, Chemical Engineering Laboratory, Trondheim, Norway 


(Received 28 


Ipril 1958) 


Abstract—The effect of the surface active agents sodium hexadecanyl sulphate and sodium 


undecany! sulphate on the rate of mass transfer, droplet oscillation and terminal velocity has 


been studied during extraction of o-nitrophenol and iodine from aqueous solutions by single 


drops of carbon tetrachloride. The effect is caused by the pure surface active agents and does 


not depend on the presence of impurities. The hydrodynamic character of the effect has been 


contirmed \ comparison with solid sph res shows that the surface active agents can make the 


rates of mass transfer and the terminal velocities equal to those for solid bodies, internal circula- 


tion, oscillation and the zig-zag path being completely eliminated. Evidence is presented 


indicating that the high rates of mass transfer obtained in the absence of surface active agents, 


about three times that for solid spheres, are not caused by the droplet oscillation, and that, 


at the most, they can be accounted for only partly by internal circulation. It is suggested that 


the explanation must be sought in the micro structure of the flow pattern. 


Résumé —Les auteurs étudient au cours de Vextraction de le nitrophenol et de Viode a partir 


des solutions aqueuses par gouttes individuelles de tétrachlorure de Carbonne leffet des agents 


tensio actifs hexadecanyle sulfate de Na et undecanyvle sulfate de Na sur la vitesse de transfert 


massique, Poscillation des gouttes et les vitesses a leéquilibre. Leffet est produit par les agents 


tensio actifs purs et ne depend pas de la présence dimpuretés, Le caractére hydrodynamique 


a été contirmée., Une comparaison avec les spheres solides montre que les agents tensio actifs 


peuvent donner des vitesses de transfert massique et des vitesses & léquilibre égales a celles 


obtenues avec des corps solides, la circulation interne Voscillation et la direction en zig-zag étant 


compléetement éliminées. Les auteurs prouvent que les vitesses de transfert de masse ¢levés 


obtenues en Pabsence d'agents tensio actifs trois fois plus plus environ qu’avee les sphéres solides, 


he sont pas occastonnees par loscillation des gouttelettes, et, au plus, elles peuvent etre interpreteées 


en partie seulement par la circulation interne, Ils suggérent que l'¢ xplication peut étre cherchee 


dans la micro structure du modéle d’écoulement. 


Zusammenfassung —Die Wirkunyg der obertliichenaktiven Stoffe Natriumhexadekanyl sulfat 


und Natriumundekanyvisulfat auf die Grosse der Stoffiibertragunyg, Tropfenschwingung und 


KEndgeschwindigkeit wurde bei der Extraktion von o-Nitrophenol und Jod aus wiissrigen 


Laisungen durch einzelne Tropfen von Tetrachlorkohlenstoff untersucht. Die Wirkung wurde 


durch die reinen obertlichenaktiven Stoffe verursacht und hing nicht ab von der Anwesenheit von 


Verunreinigungen. Der hyvdrodynamische Charakter de Wirkung wurde bestitigt. Ein Vergleich 


mit festen Kugeln zeigt, dass die oberfliichenaktiven Stoffe die Gross der Stoffiibertragung und 


der Endgeschwindigkeiten gleich denen fiir feste KOrper machen kOnnen, wobei inner Zirkulation, 


Schwingung und Zick-Zack-Bahn vollstiinding ausyeschaltet wurden. Es wird nachvewiesen, 


dass die hohen Betriige der Stoffiibertragung ohne obertliichenaktive Stoffe, die etwa dreimal 


so gross sind als fiir feste Kugeln, nicht durch die Tropfenschwingung verursacht werden und 


dass man sie in den meisten Fiillen nur teilweise der inneren Zirkulation zuschreiben kann. Es 


wird empfohlen, die Eerklirung in der Mikrostruktur des StrOmungsbildes zu suchen. 


* Present address: N. V. De Bataafsche Petroleum Maatschappij, The Hague, Netherlands. 
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INTRODUCTION 
Previous papers in this series [12], [19], [1] and 
[20] have dealt with the effect of the surface 
active agent sodium oleyl-p-anisidine sulphonate 
on the rate of extraction of iodine or o nitrophenol 
solutions with single drops of 


from aqueous 


carbon tetrachloride. The results were given in 
terms of an interfacial resistance coeflicient which 
was defined simply as the difference between the 
resistance to mass transfer in the presence and 


This addi 


tional resistance was shown to be closely related 


absence of the surface active agent. 
to the concentration of surface active molecules 
adsorbed at the interface, and as a first approxi 
mation these two quantities were assumed to be 
proportional, Ky idence was also pres nted indicat 
ing that the additional resistance caused by the 
surface active agent is in effect the result of 
some modification to the normal mechanism of 
mass transfer, and does not involve specific 
physical or chemical forces acting between the 
surface active molecules and the diffusing solute. 
It was suggested that the general nature of this 
could be explained by 


additional resistance 


assuming that the adsorbed surface active 
molecules somehow are capable of modifying the 
hydrodynamic conditions in the viemity of the 
interface. The term hydrodynamic its in this 
context to be interpreted in its widest sense as 
involving any movement of aggregates of solvent 
molecules, however small. 

In the present investigation other surface active 
agents have been used in an attempt to correlate 
interfacial resistance with the character of the 


adsorbed layer of surface active molecules. 


Lewis [17] has already from results obtained 
with his transfer cell suggested that only rigid 
films are capable of reducing the rate of mass 
Under the 


conditions prevailing in his transfer cell, Lewis 


transfer. particular experimental 


considers transfer of turbulence through the 


interface to be the predominant mechanism of 
mass transfer, and that rigid layers of surface 
active molecules interfere with this. 

The second object of the present work is to 
clarify the effect of impurities on the action of 
agents. This particularly 


surface active 


important, since the results of CULLEN and 


Davipson [2], appertaining to the effect of 
surface active agents on the absorption of carbon 
dioxide in water, point to the conclusion that 
pure surface active agents cannot alone produce 
interfacial resistance. The results of LinpLanp 
and ‘Trrsesen [19] and and 
Texsesen [1] with purified sodium oleyl-p- 
anisidine sulphonate seem, however, to contradict 
this idea. It must be pointed out, however, that, 
whereas and Davipson consider their 
effect to be of the barrier type, the latter authors 
interpret theirs in purely hydrodynamic terms. 
Although subjected to repeated recrystallizations, 
the surface active agent used by the latter 
authors cannot be claimed to have a really high 
degree of purity. On the other hand, it seems 
unlikely that 


for the effect of an agent 


impurities could be responsible 
which ts active in 
To settle 


the point conclusively it) was nevertheless con- 


concentrations as low as 107 mole 1. 


sidered necessary to synthesize surface active 
agents with less complex molecules, using routes 
which ensure the greatest possible freedom from 
homologues and other impurities. Results obtained 
during the present work with highly pure sodium 
hexadecanyl sulphate (C,, H,, Na*) show 
that impurities do not contribute noticeably to 
the effect of strong surface active agents on the 
type of resistance to mass transfer encountered 


with our system. 


There remains the important question of the 
nature of the hydrodynamic effect responsible 
for the retardation of mass transfer, and more 
particularly, of the part played by the internal 
circulation in the falling drops. Garner and 
SKELLAND [9] have shown that surface active 
agents can stop this internal circulation and 
thereby reduce both the rate of fall of the drops 
and the rate of mass transfer. Garner and Hate 
[6], however, have reported that in extracting 
diethvlamine from toluene by drops of water in 
the presence of Teepol, their technique of visual 
observation failed to reveal any internal circula- 
tion even in the absence of the surface active 
agent. These findings led them to conclude that 
mass transfer was a_ true 
Tersesen [19] 
have pointed out that with the system water- 


the retardation of 


barrier-effect. LiInpLAND and 


VOI 
| 
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carbon tetrachloride with iodine or o nitrophenol 
as solute the partition coeflicients are so greatly 
in favour of the solvent phase that the resistance 
to mass transfer is almost exclusively located in 
the continuous aqueous phase outside the drops. 
These authors did. hows ver, stress that, although 
internal circulation therefore cannot directly 
assist in the mass transfer it might well do so 
indirectly by altering the hydrodynamic con- 
ditions outside the drops. It has been made the 
third object of the present work to provide 
further evidence on this question, The effect 
of the surface active agents on the rate of fall 
of the drops will be taken as an indication of their 
effect on the internal circulation, and will be 
compared with the effect on the rate of mass 
transfer. If the rate of fall and the rate of mass 
transfer both decrease in a similar manner with 
mereasing concentration of surface active agents 
it is reasonable to assume that the surface active 
agents in both cases act through a retardation of 
the internal circulation. Lf, on the other hand. the Vv 
vary in fundamentally different ways, it must b 
assumed that retardation of internal circulation 
and of mass transfer are different. and largely 
manifestations of the adsorbed 
Already 


and Tersesen [19] have pomted out that the 


independent 
surface active molecules, LINDLAND 
sharp primary retardation in the rate of fall of 
the drops appear to reach its maximum value 
at a concentration of surface active agent about 
one-third to one-half of that necessary to develop 
fully the additional resistance. It was not at the 
time considered safe to draw any conclusions 
from this observation, but the results of the 
present investigation contlirm and extend this 
finding and make it unlikely that surface active 
agents reduce the rate of mass transfer mainly 
through their action on internal circulation. 
This conclusion is, of course, limited to cases 
where the normal resistance to mass transfer is 
located almost exclusively in the continuous 
phase outside the drops. 

GARNER [5] has also pointed to the effect of 
surface active agents on the oscillation of the 
drops as another possible mechanism for their 
ability to slow down the rate of mass transfer. 
It will be shown that with the present system 


the intensity of oscillation increases with increas- 


ing concentration of surface active agents, 


reaches a maximum and then disappears com- 
plete ly, the whole cycle being completed at con- 
centrations much lower than those necessary to 
vive the full reduction in the rate of mass transfer. 


Again, no correlation appears to exist. 


EXPERIMENTAL 


The apparatus and experimental procedure have been 
described in detail elsewhere [19]. The glass rod in the 
adsorption section of the column was rotated at 750 r p.m., 
and the effective length of the extraction section proper 
was 143em. The stop-cock on the collecting burette at the 
bottom of the column was fitted with a Teflon plug which 
did not require lubrication. The drops were formed at 
the comparatively slow rate of 1-25 per see in accordance 
with the observation of Garner and SKELLAND [8] that 
slow rates of formation promote constant drop size. As 
before 19), capillaries of thin wall glass tubing were 
found most satisfactory 

The greatest care was taken to avoid contamination 
which, even in minute quantities, interfered with the 
results Wherever practical considerations made it 
essential to use flexible connexions the joints were so 
designed that the rubber o1 plast was always separated 
from the fluid by a cushion of air as described previously 


(19). An improved version of the drop-forming apparatus 


FROM CONSTANT 
PRESSURE RESERVOIR 


| PV.G. TUBING 


PV.C. TUBING 


AQUEOUS LEVEL IN COLUMN 


CAPILLARY 


GLASS SCREENING TUBE 


Fic, 1, Improved arrangement for droplet formation. 
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shown in Fig. 1 was used. This allowed easy change of 
capillary without danger of contamination 

The 
cleaned 


followed by 


were thorough, 


column and auxiliary equipment 


sulphuric acids 
But 
pres ed uw 
the high 


surtace 


mixture of chromic and 


with distilled 


with a 


rinsing water even with 


then precautions it ocensionalls possible for 


rates 


this 


considerable periods to reproduce transfer 


obtained without addition of active agent 


being by far the most susceptible condition 
For the bulk of the 


maintained at 25 ¢ 


experiments the ter was 


thermostatically 


perature 
by circulating 
controlled water through the column jackets 

When the 


adsorption and extraction sections of the colu 


drops passed the restriction between the 
nn they were 
slowed down and subsequently accelerated in the extraction 


The drop velocities determined in the extraction 


do 


were 


section 
terminal velocities 


section net therefore represent 


Thes« determined in a column long with 


internal diameter 75cm by timing the drops over a 


distance of 180 em after they had fallen 200 c1 \ccording 


the column diameter used is not quite 


wall effect 


to [16 


suflicient to eliminate the 


MATERIALS 


Synthesis of sodium undecanyl sulphate 1 Na 


distillation 
\ccording to 


Undecenoic acid derived from castor oil by 


in vacuo was used as a starting material 


Jones and PyMawn [14] this process gives an acid practically 


The distilled 


five 


free from homologue crude acid was 


in vacuo, re-crystallized times from petroleum ether 


(40-60°C), once from ethanol and once from acetone 


The ervstallizations were carried out at In © with om 
liquer 


tinal cistilla 


part of acid to one part of solvent. and the 
After 


melting interval was 24-1 


was removed by hvdrauli 
the 


about 80 per cent melting within 24-5 245 © as compared 


pre ssure 


tion in vacuo 240 C* with 


with the m.p. of 245 given by Jones and PywMawn [14 
The 
Kean ef al 
re-crystallizations from 
melted the 
per cent within 28-6. 28-8°C as compared with Ralston and 
a4 
effected 
aluminium NYSTROD 
After repeated washing with alkali and acid, followed 
the elted im the 


with 80 per cent within 15-0165 ¢ 


acid was hvdrogenated according to 


to 


undecenot 


15 undecanoic acid which after three 


petroleum ether is) deseribed 


above, over interval with 


melting point of 28-15 ¢ 
Redu« thon to alcohol 


according to 


was with lithium 


hvdrick Brown 


distillation in vacuo undecanol 
This melting point appears to be exceptionally difficult 
Beilstein from 11 to 194 


VERKAD! 


by a 


interval 16 


vives values ranging 


et al, [29 


and 


while obtained 16-5 4 


*The melting intervals were determined by the capillary 
the 
The temperature of the water flowing passed the capillary 
the 
still 


results. 


method material being observed through a lens 


was never raised faster than min Around 


the 


melting point the rate was reduced to 0-1 ¢ 


slower rates of heating made no difference to 


and S. G 


TeRIFSEN 


The undecanol was sulphated with chlorosulphonic acid 


t}, and the resulting sodium undecany! sulphate extracted 


in glacial acetic at by the procedure of Drrcen ef al 


for eight hours with ether in a Soxhlet apparatus 


It had the form of skinny ervstalline plates 


twee 


Synthesis of sodium heradecanyl sulphate 4 16 Na 
is prepared by the method of PLeseK 
of the 


a Grignard reaction with the ethyl ester 


Palmitic acid w 
2 Ihe 


was condensed by 


chlorice undecanoic acid obtained above 


>-hvdroxy hexade 


carboxvlic acid was made from purified adipic acid with 


of cyclopentanone carboxylic acid to 


canoie acid (mp 705-797 ©) cyclopentanone 
7 152-46 C as compared with the m.p. 
The 
und treated with sodium and 
The 


wis 


melting interval 151 
of Verkape ef al. [29) of 


distilled in vacuo 


adipic acid was 


esterified 
Organic Synthesis (30 
thre 
reduced to palmitic acid with hydrazine hydrate according 
After five 
tetrachlorict 


dry toluene as described in 


hvdroxv-acid obtained from mam reaction 


to [25 recrvystallizations from ethanol, 
neetone, as 


melted 


and one from 
the 
as compared with the 
The 
for 


one from carbon 


neid acid 


> 


described for undecanon ie 
the 
highest m of ¢ 


aleohol 


and 


over interval 62-55 62-7 


28 


given by TiIMMERMANS 


reduction to was carried out as described 


ervstallizations from 
9-15-40-3 
The 


des ribed 


undecanoic acid after two re 


the interval was whereas 


acetome melting 


TIMMERMANS [28] gives mp sulphate was 


finally prepared and purified as for sodium 


undecany! sulphate, and its appearance was similar 


Carbon tetrachloride 


For the majority of experiments the carbon tetrachloride 


was purified by distillation as described previously [19], 


but in some cases it was treated with chlorine and chlorine 


dioxide according to Dickinson and LeerMakers [3]. 
This elaborate purification did not affect the results, and 


was consequently discontinued 


Other materials 


The o-nitrophenol was a erystalline Merck product, with 
The Pro 


Both substances were dissolved in distilled water 


a melting point of 45°C, iodine was of Analvsi 


quality 
containing 0-002 mole lof citric acid giving a pH of 3-15 


Hops oF ANALYSIS 


Solutions of iodine and o-nitrophenol in water 
and in carbon tetrachloride were analysed photo- 
metrically as before [19] and [1], using a Beckman 
DU Quartz Spectrophotometer. 

Dis< tT SSION 


RESULTS AND 


\ few exploratory experiments were carried 


out with v-valeric and lauric acids as surface 


active agents, but the sodium salts of hexadecany! 
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falling droplets of carbon tetrachloride, and the 


sulphate and undecanyl sulphate were used for 


the bulk of the work. In all cases the solvents temperature was maintained at 25°C, 


were water and carbon tetrachloride with the 


former as the continuous phase. The diffusing n- Valerie acid 


solute was mainly o-nitrophenol, but iodine was n-Valeric acid had no effect on the rate of mass 


used for a few experiments. The transfer was transfer in the concentration range investigated, 


always from the continuous aqueous phase to viz. 7-5 10> to 23 10°? ¢ 100 ml aqueous 


Table 1. Extraction of o-nitrophenol in presence of highly pure, synthetic sodium hevadecanyl sulphate 
and 0-002 M citric 


id 


7 


(mole 1) (mm) (em sec) m Cra (sec com) (sec em) 


R 


Average droplet volume: 18:1 0-2 mm* 


Average droplet volume: 191 0-8 

0 1M 1F 21-97 0 
19-45 
19-9 


wo 


145 


19-2 


18-05 


18-95 
18-6 


*Average of 13 observations. For Rpg: yn 


tAverage of 8 observations. For Rpg: yn = 06, 


EO 
VOL, 0-29 22-1 1-955 18-6 2-8 
In 22-1 101 31 
52/50 O-58 Is-2 21-8 1-77 53-4 7 
OOS 22-4) 1-71 0-2 
OST 21-2 1-2 14-3 
21-2 1-53 63-3 17-5 
1-16 19-7 1-47 719 26-1 
116 17-0 19-7 140 75-4 20-6 
1-45 19-85 1-45 7246 26-8 
1-75 74-6 25-8 
1-75 Is-2 19-85 1-27 36-0 
2-62 Int 125 82-8 37-0) 
22 20-1 1-27 81-2 35-4 
2-01 1-12 92-0 47-1 
201 20-0 1-15 91-7 5-0 
678 157-1 109-0 
0-678 156-5 108-4 
12:1 19-1 O-650 165-5 117-4 
12-1 19-2 0-646 164-6 116-5 
15-1 21-4 162-0 114-8 
16-7 158-1 110-0 
0-665 162-2 114-1 
30-5 O79l 85-0 87-8 
6140 19-1 0-690 158-6 110-5 
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phase. The distribution coeflicient for n-valeric 


acid between carbon tetrachloride and water 


was found to be approximately 1.55, and prior 
to the extraction experiments the corresponding 
the 


phase to prevent any appreciable net transfer of 


amount of n-valeric was added to solvent 


surface active agent across the interface. 


Lauric acid 

Lauric acid is sparingly soluble in water, and 
was consequently added to the solvent phase. 
An impure, commercial sample reduced the rate 
the 


sulphonate 


of mass transfer to approximately same 


extent as sodium oleyl-p-anisidine 
[19], but after purification the effect disappeared 


completely. 


Sodium hexadecanyl sulphate 


The results obtained with the highly pure, 


synthetic product are given in Table 1, and those 
obtained with a sulphate prepared from commer- 
Both series have 
the 
droplets 
the 


preceeding papers in this series the action of the 


cial cetyl alcohol in Table 2. 


been plotted in Fig. 2 which also shows 
the fall of the 


caused by the surface active agent. As in 


reduction in rate of 


surface active agent on the rate of mass transfer 


Table 2. 


and S. G. 

is expressed as an additional resistance defined 

as the difference between the resistance with and 

without the agent present : 
R Rey — Ryo (1) 

The 


aqueous basis. Only averages for each concentra- 


resistances are calculated on a raflinate or 


tion of surface active agent are plotted to avoid 
overloading the tigure. 

It is evident from Fig. 2 that the action of the 
pure and impure hexadecany! 


The 


DAVIDSON [2] that the presence of impurities 


sulphates are 


identical. suggestion by CuLLEN and 


constitutes an essential condition for reduction 


in the rate of mass transfer is not therefore 


generally valid. However, impurities can have 
such an effect also in this case as shown by the 
results reported above for lauric acid. As pointed 
out earlier, the effect observed by the present 
authors is hydrodynamic in nature, whereas 
CULLEN and Davipson claim to have a true barrier 
effect, and, whereas we use liquid-liquid extrac- 
tion, their system is gas absorption. 

The fully drawn curve in the lower part of 


Fig. 2 represents equation (2 and is seen 


R 
R,,/112)*" 


(2 


Extraction of o-nitrophenol in presence of sodium hexvadecanyl sulphate prepared from commercial 


cetyl alcohol, and 0-002 M citric acid 


Ns 


(mole 1) (mm*) (em sec) 


Average 
21-8* 


droplet volume 


| Rp 
Cra (sec 


Rp; 


em) 


120-0 
119-2 
117-3 
110-5 


S78 Iz 


9035 
O17 


166 
743 
735 145 
718 157 
158 
Tol 158-! 
O713 161-4 
0-689 167-4 


— 
zs 


*As for Table 1. 
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56 Is 
5 17-05 
5-6 Is 
2 
16-8 18-1 
17-85 
32-5 13-1 
17-55 
17-45 
121-6 
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Additional mass transfer resistance R,; and reduction in rate 


of fall \ v with sodium hexadecany! sulphate and 0-002 M citric acid 


Droplet size 
O Droplet size 
@ Droplet size : 


to fit 
experimental error. This equation is very similar 
to that 


anisidine sulphonate [1] and given below on a 


the observations within the limits of the 


prev iously found for sodium olevl-p- 


molar basis (3) 


R 


Ry, 


Ns 


(3) 


but the 


resistance 


The slopes are practically identical, 
the 


R,, differ by about 10 per cent. This discrepancy 


limiting values for additional 


appears to be too large to be attributed to 
experimental error, but no explanation for it 
can be advanced at present. 

The plot of the reduction in the rate of fall 
of the droplets exhibits a complex pattern. 


It appears that additions of surface active agent 


| 


18-1 pure agent 
19-1 


18-2 mm?; 


pure agent 


impure agent 


less than about 0-5 10-7 mole | have no effect 
on the rate of fall although they appreciably 
With additions 


the 


reduce the rate of mass transfer. 


vreater than 1 10-7 mole 1, however, 


reduction in the rate of fall increases very sharply 


with increasing concentration up to a_ local 
maximum, whereupon it falls and passes through 
a minimum. Subsequently it rises slowly and 
becomes ultimately constant. The local maximum 
and minimum can be explained by the action of 
the surface agent on the oscillation of the droplets. 
It was observed that the intensity of oscillation 
was slight in the absence of surface active agent, 
but increased strongly along the rising branch 
of the curve. The maximum oscillation and the 
maximum reduction in the rate of fall occurred 


at about the same concentration of surface active 
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5 
| 0 
x 60 > > > > + + 
VOL. S$ 
9 


G Boy im S. z 
HRISTENS 
ERIESEN 


. i t f 
able ction / } igh ly ure, & it} el d. 
decant pric 


and 0-002 M citric acid 


ne 107 
Rp 


(mole 1) (see cm) 


1s 
18 
18°. 
Is 
is 
Is 
20 
18 
Is 
Is 
17 
Is 
18 
18 
18 
18 
18 
17-4 
18 

18 

Is 

Is 

17-8. 

18 

18 

18 


= te 


tS 
te 


2 
25 


ms 


Rr, 
(sec em) 
21-85 170 wo 
| 21-35 1-77 7-5 
VOL 
9 
1958/7 
Mo 124 So 
Ta 19-5 1-12 95-7 50-9 
19-65 1407 50-6 
730 ) 1-12 47-2 
91-3 10-5 lol 57-7 
oo 19-45 538 
110 1-6 517 
123 O-OR87 106-1 50-3 
123 o-os2 63-1 
131 1-45 lil 65-7 
134 0-955 111 O42 
134 19-5 114 67-5 
133 ool 113 
153 115 ony 
153 owe 111 
153 O-006 116 oo" 
10-5 115 ooo 
183 1-45 008 116 69-5 
10-45 120 73-2 
274 18-65 19-5 126 
207 18-5 10-4 | 82-3 
70 1-45 145-6 
1100 18:1 19-3 145-9 


the 


action of surface 


agent, viz. 1-2 10 


With 


decreased and dis- 


7 mole 1. further 


additions the oscillations 
appeared completely when the reduction in the 
Krom 
there on the droplets fell in a pe rfectly quiet way, 
A reduction in the rate of fall due to oscillation 
has previously been reported by GARNER and 
SKELLAND [9]. The fact that the curve for the 
resistance coellicient exhibits no corresponding 
maximum and 


rate of fall passed through its minimum. 


minimum shows that in this 
case the effect of oscillation on the rate of mass 
transfer is insignificant. 

The complex pattern of the curve representing 
the reduction in the rate of fall is clearly due to 
effects, i.e. 
agent on the 
oscillation and on the internal circulation. Uf 


the supe rposition of two separate 


the action of the surface active 


the effect on oscillation could be climinated the 
effect on the 


internal circulation would be 


active agents on mass transfer in liquid-liquid extraction 


expected to follow a steadily rising curve of a 
character not unlike that of the resistance curve, 
This similarity indicates that in this case the 
action of the surface active agent on internal 
circulation and on mass transfer is parallel, 
However, the results given below for sodium 
undecanyl sulphate show that this parallelism 
is fortuitous, and that the reduction in the rate 
of mass transfer is only partly, if at all, caused 


by suppression of the internal circulation. 


Sodium undecanyl sulphate 


The results obtained with highly pure, synthetic 
sodium undecany!| sulphate are given in Table 3 
and the ay erage resistances for each concentration 
plotted in Fig. 3. 
results of Table 


solution was made 0.006 N with sodium chloride. 


In this figure are also shown the 


t obtained when the aqueous 


It is clear that the reduction of the number of 


Re 
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Additional mass transfer resistance R pj and reduction in rate 


of fall Av with sodium undecany! sulphate and 0-002 M citric acid o 


@ 0-006 M sodium chloride also added. 
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Table 4. Extraction of o-nitrophenol in presence of specially pure, synthetic sodium undecanyl sulphate 


and 0.002 M citric acid 


0-006 M sodium chloride 


ns t 
(mole 1) (mm?) (em sec) 
24 Is-4 19-6 
241 19-3 
241 19-2 
24-1 19-2 
19-5 
36-5 
610 19-2 
oi Is-4 125 
243 17-9 19-25 
245 19-25 


Rp Rp; 
m» ¢ RA (sec em) (sec om) 
1-20 
12 
141 116-6 
105 
0-969 108-8 62-2 
108-3 61-7 
Oss 125-9 
151-1 106-1 
O-743 145-5 100-5 


carbon atoms in the chain from 16 to 11 has 
greatly weakened the surface activity and changed 
the character of the curves. The results can no 
longer be desembed by equations of the type 
applicable to the sodium salt of oleyl-p-anisidine 
sulphonate and henack canvl sulphate Two 
equations were necessary for the results obtained 
in the absence of sodium chloride, equation (4) 
for the lower branch and equation (5) for the 
upper, and these are shown as solid drawn lines 


in Fig. 4. 


Rp, 
36 x 10-° 
1 — Rp, 112 
112 
exp (4) 
Ry, 
N, = 92 x 10-* (5) 


Ry 


In the presence of sodium chloride equation (6) 
was found to fit. 


Ng = 21-5 x 10-* Rw (6) 
1 — Ry 
Equations (5) and (6) are straight Langmuir 
equations whereas equation (4) Is a modified 
form. The equations can be seen to describe 
the results within the limits of exp rimental 
error. Without addition of sodium chloride there 
appears to be a break in the curve. This might 
be due to a phase transition, but sinee repro- 
ducibility was distinctly poorer with this surface 
active agent, restraint must be shown in drawing 


conclusions. 


The curve for the reduction in the rate of 
droplet fall is similar to that obtained with 
hexadecanyl sulphate. The maximum and 
minimum are clearly visible, and variation in 
droplet oscillation in relation to these was also 
the same. The striking feature of the results is 
that the reduction in the rate of fall reaches its 
terminal value when 80 10-7 mole | of surface 
active agent has been added, whereas the addi- 
tional resistance to mass transfer at this concen- 
tration Is only half developed. It is reasonable to 
conclude that the internal circulation has ceased 
when there is no further reduction in the rate 
of fall, and that the action of the surface active 
agent on the internal circulation therefore, at the 
most, can be only partly responsible for the 
reduction in the rate of mass transfer. 


Transfer of iodine in presence of sodium undecanyl 


sulphat 


It is of interest to compare the action of the 
weak surface active agent sodium undecany! 
sulphate on the extraction of o-nitrophenol and 
iodine. Such a comparison has previously been 
reported for the strong surface active agent 
1}. The 


present results with the weak agent are given 


sodium oleyl-p-anisidine sulphonate | 
in Table 5. The ratio of the additional resistance 
caused by the surface active agent R », to the 
resistance in the absence of the agent Ry, ts 
plott d in Fig. 4 for the extraction of o-nitrophenol 
and iodine. The results obtained with iodine 


234 


VOI 
1958/ 
| 


On the action of surface active agents on mass transfer in liquid-liquid extraction 


Table 5. Extraction of iodine in presence of highly pure, synthetic sodium undecanyl sulphate and 
0-002 N citric acid 


» 7 Cp, 100 
Ns 10 Rp R p; 
(mole /1) (em sec) m Cra (sec com) (sec /em) 


Average droplet volume : 85 0-2 mm? 
21-0" 


19-6 
200 
20-0 
19-45 
19-6 
1o-6 
19-45 
10-45 
19-3 
19-45 
115 
1060 2 115 


zoe 


168 
168 


*Average of 5 observations. For Rpg Itand ayn 


$00 = 


to 30 $0 70 
Fig. 4. The ratio Ry; plotted against concentration 
of sodium undecanyl sulphate. hic. 5. Terminal velocities with sodium hexadecany! 
@ for iodine sulphate. 


© for o-nitrophenol. @ without citric acid O with 0-002 M citric acid. 


were less reproducible, but the figure shows 


clearly that the action of sodium undecany! 


sulphate is the same for both diffusing solutes. 


This confirms the previous results with the 


strong agent. 


Terminal velocities in larger column. The effect of 


citric acid and salts 


The reduction in the rate of fall observed with 


sodium hexadecany! sulphate is shown in Fig. 5, 
Termin: ‘locitie ‘any 

and with sodium undecanyl sulphate in Fig. 6. erminal velocities with sodium undecany! 

: sulphate. 
© without citric acid @ with 0-002 M citric acid 
with an equivalent diameter of 0-337 cm. D with 0-002 M sodium citrate. 


The droplet volume was in all cases 20 mm? 


|... 
2-16 17-5 
1-1 75-2 34-6 
205 29-6 
1-63 90-8 50-2 
1-63 49-2 
1-62 
140 90-6 
oo 137-8 97-2 
1-11 133-6 93-0 
VOL, 1-16 131-1 90-5 
1-18 128-5 87-0 
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The rate of fall in the absence of surface active 
at 25 C. With the 
tion column the velocity was 21-9 cm sec showing 
that the 
the entrance to the 


agent was cm. sec extrac 


reached at 
the 


terminal velocity was not 


extraction section of 
column. The value of 22-8 cm. sec is considerably 
higher than the 19-8 em. see reported by Hu and 
Kintner [13] at 304°C. These 
that when plotted against the droplet diameter 
the velocity reached a maximum of 20-2 
0.35 om. 


authors found 


cm sec 
The 


difference in temperature does not explain the 


at an equivalent diameter of 


discrepancy between their results and those of 


the prese nt authors. It is ye rtinent to point to the 


possibility that the terminal velocities might have 


been affected by surface active impurities pres nt 
in the technical grade organic liquids 

On addition of sufficiently large quantities of 
the 


appeare d to 


surface active agents terminal velocities m 


all cases approa h 20-l em. see 


Irrespective of whether citric acid o sodium 


citrate were present as well. 


All curves have maxima coinciding with 


maximum oscillation and minima where the 


oscillations have ceased comple tely Phe chara 


teristic shape of the curves ts part larly well 


established in Fig. 6 containing a large number 
The 


continuation of the 


of observations. oscillations are not just a 


oscillations i duced by 
the 


droplets. This follows from the obse rvation that 


internal currents during formation of the 
in some cases the droplets fall fairly quietly at 
the start before the adsorption laver has been 
established. and then suddenly begin to oscillate 
violently further down the column. 

In addition to internal circulation and oscilla- 
tion the movement of the droplets is characterized 
by a zig-zag path. \ similar behaviour has been 
observed by Garner and for 


Whereas these 


the path became mux h more irregular on addition 


gas bubbles. authors found that 


of vaseline, the present observations showed that 


agents reduce, and 


When the 


path becann 


additions of surface active 


even climinate this tendency. 


fall 
practically straight. 


has reached its minimum the 


There were some indications that a stationary 


adsorption state had not been established even 


rate of 


and S. G. 


after a fall of 200 em. It was observed in some 
cases that oscillation started within the measuring 
distance, and in other cases that it died away. 
The surface active agents became more effective 
on addition of sodium chloride, sodium citrate, 
The effect on the 
velocity can be seen from Figs. 5 and 6 which 
that the weaker 
strongly affected. 
curves with the 


essentially unaltered, and the terminal velocity 


and even citric acid, terminal 


surface active 
The 


agent Is 
the 


remain 


show 
most character of 


maxima 


appears to approach the same limiting value on 
addition of sufliciently large quantities of surface 


active agent. Particularly illustrating is Fig. 7, 


where the terminal velocity in the presence of 


2.2 10-7 mole 1 of sodium undecanyl sulphate 


concentration of 
this 


has been plotted against the 


added sodium chlonde. Even with very 


Terminal velocities with ng 22 « 107 mole 


Fie 7 
of sodium undecany! sulphate plus varying additions of 


sodium chloride. 


low concentration of surface active agent the 


limiting condition is reached when the solution 
is made 0.004 M with respect to sodium chloride. 


In this connexion it must be remembered, that 


all experiments with o-nitrophe nol as diffusing 


solute, and also the present experiments with 


iodine. have been carried out with solutions 


containing 0-002 mole | of citric acid, 
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Comparison with solid spheres. The nature of the 


action of surface ac tive agents 
The 


diameter of 


terminal velocity of a sphere with a 
0.337 cm 


the method recommended by 


was calculated to be 
23-8 cm/sec by 
[16]. 

figure with the average terminal velocity of 20-1 


It is of interest to compare this 


obtained on addition of comparatively larg 


quantities of surface active agents. The results 
of Metuvs and Trrsesen [21] show that the 
shape of the droplets generally approximate to 
oblate spheroids. With sodium oleyl-p-anisidine 
sulphonate the ratio of the short to the long axis 
0-82 at a concentration of 1-2 « 10 

The 
of the terminal velocities for two bodies of the 
the 


was ba 


mole per litre and 0-78 at 8-2 10 ratio 


same material moving in same medium is 


given by equation (7) 


Yea 

where Ay 

coetlioent. 

of the 


is the projected area and C the drag 
For an oblate spheroid and a sphere 
this 


volume reduces to 


spheroid 
(8) 
Us (sphe roid b ( sphere) 


same 


i (sphere) 


As a first approximation the drag coetlicients 
for spheres and spheroids at the same Rey nolds’ 
By inter- 


polation in the table of drag coetlicients given by 


number are assumed to be the same. 


LapPLe the velocity ratios are calculated for 
the two additions of surface active agent to be 
0-89 and 0-87 respectively, whereas the ratio 
between the observed terminal velocity and that 
calculated for an equivalent sphere is 0.84. This 
difference between the calculated and observed 
ratio of less than 5 per cent can probably be 
explained by the fact that the drag coeftlicients 
for oblate 


than for spheres, but, unfortunately, measured 


spheroids undoubtedly are larger 


drag coeflicients for spheroids are not available 


in this range of Reynolds’ numbers. However, 


the conclusion is justified that the movement of 


a droplet in the presence of sutfliciently large 


amounts of surface active agents 


differ from that of 


appears to 


an equivalent sphere solely 


All other effects such as 


internal circulation, oscillation and the zig-zag 


by its spheroid shape. 


path have disappeared. 

Sremes [26], using the data of Sruke [27] on 
the rise of air bubbles in water containing various 
amounts of capronic acid, has shown that the 
limiting rate of rise on addition of sufliciently 
large quantities of the acid, about 0-01 mole 1, ts 
approximately the same as for solid spheres, 
provided the Reynolds’ number is below about 
250. Above the deformation of the 
bubbles the 


coetlicient to increase above that for solid spheres. 


this value 


was too great and caused drag 


\ similar comparison with solid spheres can 
Data for 


mass transfer to solid spheres in liquids do not 


be made with respect to mass transfer. 


seem to be available*, but Suerwoop and Picrorp 
[25] have correlated data for solid spheres and 
droplets falling in air by 


Np, X 


Re 


plotting Ny, against 
The curve recommended by these 
Fig. 8 together with the 
corresponding curve recommended by McAdams 
The 
liqguid droplets in the present and previous work 
Fig. 8. It 


authors is shown in 


for heat transfer. results obtained with 


in this series are also shown in can 
be seen that the limiting mass transfer coetlicients 
obtained with comparatively large additions of 
the 


mended by Suerwoop and Picrorp whereas the 


surface active agents fall on curve recom- 


results without addition of surface active agents 
The line through 


the points obtained without addition of surface 


are about three times as high. 


active agents is drawn parallel to the Suzrwoop 
and Picrorp line for solid spheres and in a 
distance of 2-8 on the logarithmic scale. 

In Fig. 8 are also shown the results of Guyer 
and Prisrer [10] appertaining to absorption of 
carbon dioxide from gas bubbles rising through 
water. For bubble diameters larger than about 
03cm the points fall on the same line as the 
present extraction results without surface active 
agents present, but for smaller bubbles they fall 
considerably below and appear to approach the 
line for solid spheres. 

The results have also been compared with those 
of Hrertses, et al. [11]. These authors plotted 


Ny, as a function of NN», Ns.°* for the 


*See note added in proof (p. 240). 
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HEAT TRANSFER 
Mc ADAMS 


MASS TRANSFER 
SHERWOOD AND PIGFORD 


5 
Naot Nsch OF Nee a Np, = 


Fic. 8. Comparison with solid spheres 


without addition of surface active agents [19], 


and [20 
wtive 


with comparatively large additions of surface 


agents [19 1] and [20 
results of Guyer and Ptister with bubbles of carbon 


dioxide [10). 


continuous phase for droplets of water im tsebu 


tanol and of isobutanol in water = together 


with the 
referred to above. 


observations of Guyer and Prisrer 


The pres nt results fall among 
theirs, but on a much steeper curve, more in 
line with those for carbon dioxide. The poimts 


of Guyer and Prisrer and of the present iuthors 


appear to correlate better by the plot proposed 
by Surrwoop and Picronp, but further work ts 
necessary before generally reliable correlations 
can be arrived at. 

active 


The nature of the action of the surface 


agents can now be understood. It consists 


simply in eliminating all effects associated with 
The 


make the liquid droplets equivalent to 


the liquid boundary. adsorbed molecules 


solid 


spheres, except for the change in shape which 


and S. G. 


little effect on mass transfer. 


These results rule out the possibility that the 


appears to have 
adsorbed surface active molecules form a barrier 
which interacts with the diffusing solute molecules 
and thus retards their passage through the inter- 
face. The previously reported conclusion that 
the effect is hydrodynamic in character has thus 
been contirmed., 

The problem is no longer to explain why the 
rate of mass transfer is so slow in the presence of 
surface active agents, but to explain why it is 
absence of these. 


so fast in the The present 


results, and particularly those with sodium 


undecanyl sulphate, show that droplet oscilla- 
that 


internal circulation is at most capable of pro iding 


tion has only an insignificant’ effect, and 


part of the explanation. This is in agreement 
with the observation by Garner and Hae that 
the surface active agent Teepol reduced the rate 
of transfer of diethyvlamine from toluene by drops 
of water, also when there was no internal circula- 
before the addition of the The 
their has been 
[21] as 
being caused by impurities present in the T'eepol, 


tion agent. 


maximum im resistance curve 


explained by and TrERJESEN 


and there is no reason to doubt that the same 
mechanism is responsible in their case as in ours; 
that is that their effect is also of a hydrody namic 
rather than of a barrier type. 

The situation can be summed up in the form 
of a paradox: The high rates of mass transfer 
with pure systems must be due to some feature 
of the hydrodynamics of the falling droplets 
unconnected with oscillation and the zig-zag 
path and at least partly unconnected with the 
internal circulation. If the whole explanation 
cannot be found in the macro picture of the flow 
pattern it must be sought in its micro picture. 

The concept of transfer of turbulence developed 
by Lewis [17] to explain the results with his 
transfer cell does not appear to be directly 
applicable to falling droplets. In his transfer 
cell the rate of mass transfer was independent 
of the molecular diffusivity, whereas the present 
with iodine and show 


results o-nitrophenol 


proportionality with the diffusivity to the power 


of about 0-4 also in the absence of surface active 


agents. The interfacial turbulence of Lewis and 
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Pratr [18] is not directly applicable either. 
These authors found that this effect occurred 
only when the transfer of the solute was accom- 
panied by evolution of heat*, whereas the present 
results show no difference between the exothermic 
transfer of o-nitrophenol and the endothermic 
transfer of iodine except that which is caused by 
the difference in their molecular diffusivities. 
The present results also throw some light on 
the more general question of establishment of 
equilibrium at phase boundaries. It is rather 
striking that the mass transfer coeflicients for 
droplets in presence of surface active agents in 
Fig. 8 fall so well on the curve recommended by 
Suerwoop and Picrorp for solid spheres and 
These 
The fact that so 


many different processes can be correlated by 


droplets falling in air. involve both 


evaporation and sublimation. 


purely hydrodynamic variables indicate that, if 
true interfacial resistance of the barrier type 
exist at all, its contribution to the total resistance 
to mass transfer is unlikely to be greater than 
about 20 per cent, which is the spread of the 
experimental observations used by SueRwoop and 
Picrorp in establishing their correlation. At least 
some of this spread is due to experimental 
error and possible imperfections in the method of 


correlations. 


CONCLUSIONS 
(1) The action of strongly surface active 
agents on mass transfer does not depend 
on the presence of impurities, but is 


caused by the pure compounds, 


(2) Whereas sodium oleyl-p-anisidine sul- 
phonate and sodium hexadecany! sulphate 
exert a strong action on both terminal 
velocity and rate of mass transfer, the 
action of sodium undecany! sulphate ts 


much weaker. 


(3) The limiting values of the additional 


resistance to mass transfer and _ the 


reduction of the terminal velocities 
appear to be approximately the same for 
all surface active agents when sufliciently 


large quantities are added. 


*See note added in proof (p. 240). 
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(4) 


(6) 


(3) 


(9) 


(10) 


In the case of strongly active agents the 
previously developed adsorption isotherm 
gives a satisfactory description of the 
mass transfer results, whereas equations 
of the Langmuir type are adequate for 
the weaker agents. 


The additional resistance to mass transfer 
caused by the relatively weak agent 
sodium undecany! sulphate, expressed as 
per cent of the resistance without addi- 
tion of this agent, is the same both with 
o-nitrophenol and iodine as diffusing 


solute. 


The surface active agents act strongly 
on the oscillation of the droplets, and 
this is reflected in the terminal velocities, 
but not in the rate of mass transfer, 
indicating that oscillation is of little 
or no importance for mass transfer. 


With sodium undecanyl sulphate the 
limiting terminal velocity is reached with 
much smaller concentrations than the 
limiting resistance to transfer. 
Taking the terminal velocity to be an 
internal 


indicator for the degree of 


circulation in the droplet, this shows 
that the 
must at least partly be due to other 


reduction in mass. transfer 


factors than elimination of internal 


circulation. 


On addition of sufficiently large quantities 
of surface active agents the droplets 
behave as solid spheres or spheroids 
with respect to both terminal velocity 
and rate of mass transfer, confirming that 
the effect is hydrodynamic in character. 
The rate of mass transfer without surface 
active agents present is about 2-8 times 
that for equivalent solid spheres or about 
equal to that observed by Guyer and 
Prister [10] for 
dioxide from gas bubbles with diameters 


transfer of carbon 


above 0-3 

Conclusions (6), (7) and (8) imply that 
the high rates of mass transfer observed 
in the absences of surface active agents 
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are not due to droplet oscillation, and at 
the most only partly to internal circula- 
tion. The liquid-liquid boundary must 
give rise to a flow pattern which affects 
greatly the rates of mass transfer without 
affecting the terminal velocities. 

(11) The phenomena of transfer of turbulenc 
as described by Lewis [17] do not appear 
to be directly re sponsible. 


*Since this paper was prepared the authors have received 
the paper by F. H. Garner and R. D. Suckling (4. 1. Ch. E 
Journal, 1958, 4. 114) on mass transfer from soluble solid 
spheres. Their results fall close to the curve recommended 
by Sherwood and Pigford and shown on Fig. 8, and thus 
give strong support to the conclusions of the present paper 

The work of K. Sigwart and H. Nassenstein (1D / 
Zeitschr. 1956, 98, 453) has shown that the interfacial 
turbulence, first observed by Lewis and Prarr [18 
occurs not only with exothermic extractions, but also with 
endothermic Interfacial turbulence and eruptions 
cannot therefore be rulled out and must be further con 


sidered us pros bole \planations of the hvdrod 


type 
for the high transfer coeflicients obtained in the absence 
of surface active agents 

Both these aspects have been discussed in some detail 
in a paper by S.G Terjesen to be published in Dechema- 
Monographien Vol. 32. 
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NOTATION 


A, = projected area of droplet on plan 
vertical to direction of movement em* 
a long axis of oblate spheroid em 
b = short axis of oblate spheroid em 
( drag coeflicient Dim. less 
Cr = solute concentration in solvent phase my 
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Sh 


equilibrium concentration in solvent phase 


solute concentration in aqueous phase 


mg 
average of initial and final concentra- 
tion in aqueous phase my em 
specific heat cal 
diffusivity sec 
diameter of equivalent sphere em 
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extraction coefficient on a time and 
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heat conductivity cal /em sec ( 
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resistance on a time and raflinate or 
aqueous sec om 
resistance without surface active agent sec ‘em 
limiting resistance on addition of 
relatively large amounts of surface 
active agents sec oom 


additional resistance caused by the 


surface active agent Rp . Rro sec om 
droplet volume moe 
velocity of fall em sec 
reduction in velocity of fall em (see 
terminal velocity em, sec 
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Viscosity vy sec cm 
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On the rate of mass transfer from a gas to a moving liquid film 


A. 


Department of Chemical Engineering, University of California, Berkeley, California 


(Received 20 August 1957) 


Tur need for calculating rates of mass and heat 
transfer to two-phase mixtures in motion arises disorganized and often unpredictable relative 


operations, but the complexity of the phenomena 


Abstract—A theoretical expression is derived for the rate of mass transfer of a substance from 
the gas phase to a liquid film, in which it is highly soluble, flowing down a vertical surface. It 
is assumed that the gas phase is in laminar motion. The problem is solved for both the special 
case of dilute gaseous mixtures, and the more general case of high rates of mass transfer. It is 
shown that when the Schmidt number is around unity, the Nusselt number is again proportional 
to the square root of a Reynolds number as in all laminar single phase flows, except that a new 
characteristic velocity has to be defined, and that the expression for the Nusselt number for 
mass transfer to a stationary vertical surface is also applicable for mass transfer to a falling 
liquid film, if this new characteristic velocity is used. The same argument holds of course for the 
analogous heat transfer problem, that is the transfer of heat from a surface through a moving 


thin liquid film into a laminar gas stream in motion, 


Lauteur deduit une expression theorique pour la vitesse du transfert massique d'un gaz 


Résumé 
aun film liquide 
de tilm le long d'une surface verticale. Ll suppose que la phase gazeuse est en écoulement laminaire. 


La solubilite du gaz dans le liquide est trés grande et ke licquricte coule sous forme 


Ce probleme est resolu d'une part pour des melanges de gaz dilués et d'autre part pour le cas plus 
général des grandes vitesses de transfert de masse. Quand le nombre de Schmidt est environ 
égal a I, Vindice de Nusselt est proportionnel a la racine carrée du nombre de Reynolds comme 
dans tout écoulement en phase laminaire Cependant une nouvelle vitesse caractéristique doit 
étre alors définie et Vexpression de Vindice de Nusselt pour le transfert de masse a une surface 
verticale stationnaire est applicable au transfert de masse a untilm liquide descendant. Le méme 
argument est valable pour le probleme analogue au transfert de chaleur dune surface a un 


courant de gaz laminaire par Pintermeédiaire d'un film liquide en mouvement 


Zusammenfassung -Fiir die Geschwindigkeit der Stoffiibertragung ciner Substanz aus der 
Gasphase in einen entlang einer senkrechten Wand strOmenden Flissigkeitstilm, in dem diese 
Substanz leicht loslich ist, wird ein theoretischer Ausdruck abgeleitet. Es wird dabei angenommen, 
dass die Gasphase laminar stromt. Das Problem wird sowohl fiir den speziellen Fall der verdiinnten 
Gasmischungen als auch fiir den allvemeinen Fall hoher Betrige der Stroffiibertragung gelést. 
Es zeigt sich, dass bei Werten der Schmidtzahl in der Nahe von 1 die Nusseltzahl proportional der 
Quadratwurzel einer Reyvnoldszah! ist wie bei allen laminaren EinphasenstrOmungen. Es wird 
hierbei allerdings cine neue charakteristische Geschwindigkeit detiniert und der Ausdruck fiir 
die Nusseltzahl der Stoffiibertragung an eine ruhende senkrechte Fliche ist auch anwenbar fiir 
die Stoffiibertragung an cinen fallenden Fliissigkeitstilm, wenn die neue charakteristische 
Geschwindigkeit verwendet wird. Die gleiche Losung gilt naturgemiiss auch fiir das analoge 
Wiarmetibergangsproblem, niimlich fir den Wirmetibergang von einer Fliche durch einen 


strOmenden diimnen Film an cinen laminaren Gasstrom, 


all too frequently in the chemical engineering — spatial distributions of the two phases, that is 


involved is such that in general a prior predictions operating conditions. 


in this field cannot be made with any degree of 


confidence. This is due primarily to the highly 


the phase diagrams, which exist under the usual 


Thus, it is well known that when a gas-liquid 
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mixture is made to flow through a tube, more or 


less random and non-reproducible flow con- 


figurations, loosely referred to as bubble flow, 
are observed. It 
that no 


mathematical model for predicting the behaviour 


slug flow, annular flow, ete., 


is not surprising therefore general 
of two-phase mixtures has emerged so far and 
that even the empirical correlations available in 
the literature are neither general nor very reliable. 

In an attempt to arrive at a better under- 
standing of heat and mass transfer in two-phase 
gas-liquid mixtures, a simple mathematical model 
is set up and solved in detail. As will be seen, 
the purpose of this study is not so much to cover 
as general a physical situation as possible, but 
to examine carefully an idealized and well 
defined model in the hope that the solution of a 
simple problem will yield results applicable to 
more general S\ stems and shed some light on the 
behaviour of more complicated two-phase flow 
phenomena. 

In this paper, a theoretical expression will be 
presented for predicting the rate of transfer of a 
given chemical substance from a moving gaseous 
phase to a flowing liquid film, in which it is 
soluble. 


highly This occurs for example in a 


and im. certain 


When the 


essentially stationary the rate of mass transfer 


wetted wall absorption tower 


other industrial operations. film is 
can in many instances be easily predicted, since 
the problems are of course analogous to the well 
worked out cases of mass and heat transfer from a 
solid surface to a moving fluid. To date however 
no attempt has been made to determine how the 
rate of mass transfer in such operations would 
be affected if the liquid surface were to move 
with a velocity at least comparable to that of 
the gas stream. 

The purpose of this paper, then, is to study 
Such 


a treatment is at present possible only when the 


theoretically the effect mentioned above. 


gas is in laminar motion and when the liquid film 
flows smoothly along a surface. Momentum and 
mass transfer between co-current laminar streams 
have been investigated in some detail by Lock 
[2] and by Porrer [3}. 


contains an extensive bibliography on the subject. 


The latter article also 


The problem treated here is however both simpler 


and in a sense more general. Thus, on one hand 
the liquid stream remains essentially uninfluenced 
by the motion of the gas phase and it will be 
assumed that, for very soluble substances, the 
resistance to mass transfer is confined entirely 
to the gas stream. On the other hand, the treat- 
ment has been extended to cover the case where 
the rate of the 


normal component of the velocity to the surface 


mass transfer is high, so that 


of the liquid stream cannot be neglected. 

Consider then a vertical flat plate. Gas and 
liquid flow co-currently, and the liquid is assumed 
to occupy a thin film over the surface of the plate. 
highly 
the 


A substance. soluble in the liquid, is 


which is in laminar 
inert. It 
of course that unless the film is exceedingly 
thin it 
assume 
im the 


commonly observed on the liquid surface can in 


transferred from 


vas 
motion and contains an is well known 


will not be laminar. But since we shall 


that the 


gas phase alone, the rippling which is 


resistance to mass transfer is 


this case be disregarded. 

This phenomenon is described by the standard 
the Navier- 
Stokes equations, the continuity equation, and 


set of hydrodynamic relations 


the diffusion equation. However, it is usually 
permissible to simplify the mathematical treat- 
Prandtl’s 


boundary layer equations [4] which are, 


ment and use instead well known 
when 


applied to the gas phase, 


The symbols have the following meaning : 


the distance parallel to the liquid surface. 
the normal distance from the liquid surface. 
gas phase velocity components in the z 
and y direction respectively. 
the concentration of the diffusing species 
in the gas phase. 

D the diffusion coefficient. 

vy = the kinematic viscosity of the gas. 
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yu du 
u + v— (1) 
dua oY oy* 
du oi 
0 (2) 
or 
oc oc 
u Ly - D (3) 
Mw dy dy? 


A. 


In addition, the solution to equations (1-3) must 
satisfy the boundary conditions, which for dilute 


mixtures are 


at y O u 1 0. ¢ c 


at y © and also at u (4) 


where u, = the velocity of the liquid surface, 


which for practical purposes is 


constant, 
Uv the velocity of the gas stream far 


from the surface. 


the concentration of the diffusing 


substance far from the surface. 


c the equilibrium concentration of the 
diffusing substance when in contact 
with the liquid (for the purposes of 


this discussion ( 


It is quite evident now that u and vr must be 


obtained first, from the solution to equations 


(1) and (2), and then substituted into equation 
(3) before the rate of mass transfer can be pre 
dicted. Actually, 


identical with the well known boundary 


equations (1) and (2) 


are 
laver 
equations for a flow past a flat plate (4), the only 
difference being that the velocity component u 
is not zero at the liquid surface (y 0) but is 
equal to a constant wu, 

Equations 1) and (2) were originally analvsed 
by Blasius (with u, 0) who noticed that a 
could reduce them 


similarity transformation 


to an ordinary differential equation which be 


then solved numerically. Surprisingly enough, 
the more general problem with u,F- 0 can also 


be solved by the Blasius transformation. Thus, if 


7 "and u U. f' (5) 
vi 


Equations (1) and (2) reduce to the single ordinary 


differential equation 


= © (4) 
with f=0@, at » = 0 


Furthermore, it is casy to show that if u and v 


are expressed in terms of the function {(»), then 


equation (3) can be transformed into 


Se .dé 


+ 0 (7) 
dy? 2 dy 
where @ “* and Se 
so that, at » 0 0 and at » co 


The solution to equation (7) is readily shown to be 


dy 


and therefore the rate of mass transfer at the 


liquid surface is 


< va 
ay 


Finally, the Nusselt number Nu becomes 


N 

dy 


which reduces to a particularly simple form when 


AY 1, in which case 


Nu [x2 


It is clear then that the Nusselt number for 


(10a) 


mass transfer can be calculated from equation 
(10) if the solution to « quation (6) can be obtained. 


This can be accomplishe dl 


(1) Exactly by a numerical procedure, or 


(2) Approximately, by the von Karman- 


Pohlhausen momentum integral method. 


The integral method will be discussed first. 
because it leads in this case to a simple and 
This 


apphed to 


intriguing result. procedure has been 


repeatedly boundary layer  fhuid 


2 
9 
alll 
| 


mechanical problems and is known to be in general 
quite accurate. The details can be obtained from 
any standard reference [4]. 

It is simpler in this case to start from equations 
(6) and (7)—although the usual procedure has 
been to go back to equations (1-3) — and attempt 
to satisfy them in the mean rather than solve 
them exactly. Thus, to replace the boundary 
layer equations one obtains the two momentum 


integral equations 
2/" (0) ( dy 0 (11) 
and 


20° (0) Se | {0 dy 0 (12) 


These equations do not have of course a Unique 
solution. It is customary to satisfy them by 
assuming a velocity and a concentration profile, 
each of which contains an adjustable shape 


factor. Thus. two reasonable and commonly 


used profiles are 


2 


for » 4, and f 1 for » é 
37 
0 — for 7 or (14) 
3 
-~' 7 - 7 
and 0 1 for » or 


which satisfy the required boundary conditions 
and also the two differential equations (6) and 
(7) at 7 0. 


The two shape-parameters, 6 and 


dy, can be determined from equations (11) and 
(12). It is found that 


and that AH, (A) 
V1 11692 
l 1-692 u 
— 
7-180 Se 

u Or 
where u, =", A rand 
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3A 3A? 3A‘ 
3A? 3A’ 
+ 4 4 ) for A l 
20 
(19) 
3A 3 3 3 , 


Therefore, if one defines as the characteristic 


velocity 


1-692 (17) 


one can readily show that the expression for the 
Nusselt number becomes 
0-323 Us 
Nu | Sc) (18) 
A(O, Se) 
where 
A(O, Se) 
(u,. Se) (19) 
A(u,, Se) 
It is interesting to note that the function 
@ (Uy, Sc), which can be evaluated from equation 
(15) and which is shown plotted in Fig. 1, is 
relatively insensitive to u, for most gases, that is 
for those with Schmidt 
So that the Nusselt number for such two-phase 


numbers close to 1. 


flows is again proportional to the square root of a 
Reynolds number, a relation which is well known 
to hold for single phase laminar flows, the only 
difference being however that the characteristic 
velocity of the system is a linear combination of 
the velocities of the gas and the liquuid, as given 
by equation (17). According to the momentum 
integral approximation, therefore, the expression 
for the Nusselt number for mass transfer for dilute 
mixtures, and for Se ~ 1, to a stationary vertical 
surface is also applicable for mass transfer to a 
falling liquid film if the characteristic velocity U, 
is used, The same argument holds of course for 
the analogous heat transfer problem, where the 
transfer is from a surface across a moving liquid 
film into a laminar gas stream and the resistance 
of the film is negligible. 

The momentum integral method is however 
only approximate and before the results derived 
from it can be accepted they should be compared, 


if possible, with those of the exact solution. This 
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A. Ackivos 


FIG | PLOT OF THE 


1, for then 
equation (10a) can be employed which requires 
the of f'' (0) of 


was obtained solution 


can be done most readily for Se 


only value as a function 


(0) 
of equation (6), and a comparison between the 


from a numerical 


exact and the approximate results is shown below. 


Table 1. Values of Nu Je for Se l as a 
ofl 


function of u, 


Exact Ippro. imate Error 
0-332 
O472 


O-4+42 
480 
0-682 


Moreover, as u, > © which occurs when 


bulk of the gas is essentially stationary, the exact 
expression is 


Nu O-445 J 


whereas the approximate method predicts 


Nu 0-420 


Ur 


u, 
\ 


FUNCTION 


It therefore that the Karman- 


Pohlhausen method, when applied to the present 


is secn von 
problem, gives results of acceptable accuracy, 
especially when it is recognized that the approxi- 
mate results depend somewhat, although to a 
slight extent it is true, on the profiles which 
were chosen to represent the velocity and the 


temperature distributions. 


Mass TRANSFER 


Hicu Raves 


Where the rate of mass transfer is high, them 


as was shown by Srauouae [5, 6], Eckert and 
ScuNEIDER [1] and others, the interaction between 
the equations of motion and the diffusion equa- 


tion cannot be neglected. Thus, even though 


equations (1-3) still hold, the boundary conditions 


given by equations (4) must be replaced by : 


the 


at y 


u c, and (p 


at y © and also a 


It is interesting to note that again the Blasius 
transformation will reduce the problem to the 
solution of two ordinary differential equations 
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which are now coupled on account of one of the . 
(a) for A <1 


above boundary conditions, Thus, if 


»=Yy J u =U, f'(y) and = @ (») A(; (1 — 


va 


s 


it is found that equations (6) and (7) remain 


unchanged, but that the boundary conditions 


become 


This system of equations will now be solved by 


the momentum integral method. The assumed 


20A 280 A3 


f'’-u, +(1 My) 
3 
7) ' for <8 (toa 12a? * 
1258 86283 
and f 1 for » 
=O7 =O7 7 7 b b, ( (26b) 
for 130A B0A2 105A} | 
7) O7 
7 J J Therefore, before the Nusselt number can be 


The two constants b, and b, are chosen so that caleulated, A, b, and 6, must be found from 
the assumed profiles will satisfy equations (6) equations (23) and (26) as functions of the three 


and (7) at the liquid surface, i.e. » 0. Thus parameters of this problem u,, 2, and Se. Unfor- 


tunately, the profiles given by equation (22 


b, cannot be used for n, > 2 3, because b, and b, 
2 My become imaginary, and the profiles chosen by 

am) SPALDING [6] must be employed. This extension 

n, ( tom will however be omitted here, as it adds very 

- little to the present development. 
b, 15 (15 4 bs) The problem has again a simple solution when 
(b, + 3/2 — 44 Se/n,) Se 1, for then it can be shown that A l 
where A = m also, and that b, = 6,. Thus, when Se = 1, 


and therefore, from the momentum integral equations (11) and (12) 


Nu; U. 2+b, /u, (3/8 — 6,/12) + (1 u,) (39 280 — b, 140 — 105) (24) 
A 3 (1 ny SecA + 2h, 2h,n, /ScA) 


A is obtained from 


8 + 2b, 2 
Sc AH, = (3 a) (I ad b, b, (25) 


3(1 — n,/Scd) + 2b,— Sed 8 12 140105 
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0 
2 3 4 
at » =1,0=1 (21) 30) 105 
(b) for A l 
Cs BA b, 
9g 
|| 


A. AcrRIvos 


2b, 


3 


0-323 | 1+ 


where 


39 b, 


140 105 
2b, 
| 1 + (] ny) 

3 


(28) 


According to equation (26) then, for high rates 
for Se ~ 1, the 
Nusselt number is again proportional to the 


of mass transfer as well, and 


square root of a Reynolds number, and the 
expression for the Nu number for mass transfer 
to a stationary vertical surface is also applicable 
for mass transfer to a falling liquid film, if the 
characteristic velocity, l 


given by equation 


(28), is used. 

Finally, a few calculations were performed to 
check the accuracy of the momentum integral 
method of solution as applied to this problem. 
It is known that this approximate approach 
On the other 
hand, when u, 1, a simple exact solution can 
Thus, for Se 


gives reliable results when u, = 0. 


be obtained. 


(a) From the exact solution: 


Nu 0-564 J 


[1 + 06366 n, + 04923 


(b) From the momentum integral method 


Nu — 0-530 [ee 


[1 + 0625 ny + 0-500 


The therefore to 


satisfactory. 


agreement is seen 


CONCLUSIONS 


An expression has been derived for the rate of 
mass transfer of a substance from the gas phase 
to a moving liquid film, in which it is highly 


soluble. It was assumed that the gas phase is in 


ir 


laminar motion. The problem was solved for 


both the special case of dilute gaseous mixtures 


Vas 


and the more general case of high rates of mass 
It was shown that when the Schmidt 
the Nusselt 


is again proportional to the square root of a 


transfer. 


number is around unity, number 


Reynolds number, and that the expression for the 
Nusselt number for mass transfer to a stationary 
vertical surface is also applicable for mass 
transfer to a falling liquid film, if the characteristic 
velocity U, given by equation (28) is used. 
The same argument holds of course for the analo- 
gous heat transfer problem. 

The mathematical equations were solved by the 
approximate von Karman-Pohlhausen momentum 
method. It 
comparing, wherever possible, with the results 
that this 


method is of very satisfactory accuracy for the 


integral was shown however, by 


of exact calculations approximate 


problem considered, 
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NOTATION 


velocity of the gas parallel to the plate and 
outside the boundary layer. 
and horizontal 


respectively, the vertical 


velocity components divided by U,.. 
velocity of the liquid surface u 

position co-ordinates along the vertical surface 


9 divided by 


and normal to it, respectively. 

concentration of the solute in the gas phase, 
in moles per unit volume. 

respectively, the concentration of the solute 
outside the boundary layer, and at the liquid 
surface (c, 0). 


dimensionless concentration. 


x 


vy D, the Schmidt 
diffusion coeflicient. 


number, where D is the 


Sand 5, — shape factors for the velocity and the tempera- 
ture profiles respectively. 

oT 

defined by equations (16) and (26). 

defined by equations (17) and (27). 

defined by equation (19). 
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Abstract 
control analysis are applied to the continuous flow stirred tank reactor 


In this paper some of the statistical techniques which have become part of modern 
Care is taken to present 
these methods without assuming the reader's prior knowledge, for, though they are widely 
used in other fields, they may as vet be unfamiliar to the chemical engineer. It is assumed that 
the system never departs very far from its steady state so that linearized equations are valid. 
then make 


parameters on the characteristics of the output. 


One may a statistical design of the reactor and determine the effect of various 


Résumé Dans cet article les auteurs appliquent au réacteur a cuve, agite a marche continue, 


les méthodes utilisées dans les analyses de contréle modernes. Ces meéthodes, largement utilisées 
dans certains domaines mais encore peu familiéres 4 lingénieur chimiste, sont présentées de telle 
sorte qu'elles ne nécessitent pas de connaissances particuliéres du lecteur. Les auteurs supposent 
que le systéme ne s'écarte jamais sensiblement de son état de régime, de sorte que des équations 
linéaires sont utilisables. Une étude statistique du réacteur est alors rendue pcssible, ainsi que 


la détermination de Vinfluence de divers paramétres sur les caractéristiques du produit obtenu. 


Arbeit statistische Methoden, die zur modernen 


Behandlung von Regelungsvorgingen vehoren, auf den kontinuierlich durchstrOmten Reaktor von 


Zusammenfassung In dicser werden cinige 


der Art des Riihrkessels angewandt 
Methoden dem 


Gebieten als weit verbreitet gelten miissen 


cliese 


seinem Beharrungszustand abweicht 
statistisehe 


auf die charakteristischen Grossen am 


INTRODUCTION 
Or the many interesting problems connected 
which 
the 
development of appropriate mathematical tools. 


with stirred tank reactors there is one 


has recently become tractable through 
This is the probl m of describing the behaviour of 
the reactor under the smal! random fluctuations 
which occur even under the strictest of controls. 
The values of input flow rate, temperature and 
but 


undergo fluctuations of a random nature because 


“actant concentration are seldom steady 


changes in ambient conditions, valve chatter 


and a multitude of other variations too complex 
The effect that this will have on 
the output may be studied by considering the 


to disentangle. 


input variables subject to small random fluctua- 


tions and using the techniques which have been 


Vertahrensinvenieur bisher 


so dass linearisiert« 


Vorkenntnisse werden beim Leser nicht vorausgesetzt, da 


wenig vertraut sind, obwohl sie auf anderen 


Es wird angenommen, dass das System nur wenig von 


Gleichungen gelten, Dann moge eine 


Berechnung des Reaktors durchgefiihrt und die Wirkung verschiedener Parameter 
\ustritt bestimmt werden. 


developed by electrical engineers and mathema- 
ticians for the study of linear systems and random 
functions. 

Since the inputs are known only statistically 
it is not possible to obtain more than a statistical 
But 


siderable value, for given thes tatistical properties 


description of the outputs. this is of con- 
of the input, including its standard deviation, one 
may calculate the standard deviation and other 
properties of the output. Such a calculation would 
show what limits could be tolerated in the input 
to achieve a product of given uniformity and 
show how the control parameters should be 
adjusted to accomplish this, 

As the mathematical techniques are compara- 
tively new and may be unfamiliar to the chemical 
engineer, a brief outline of the method of describ- 
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ing and using random functions is first given. 
It is shown how the random properties of the 
input may be obtained from strip chart data and 
how the properties of the output may be 
calculated. The method is illustrated by a 
numerical example. 


1 Tue STATISTICAL DESCRIPTION 
or Ranpom FuncTions 

The behaviour of a continuous flow. stirred 
tank reactor when certain of its input variables 
are subject to random fluctuations is to be 
studied. It is therefore necessary to specify the 
class of random functions that will be considered 
and the way in which they may be described. 
It may be stated at the outset that the class of 
fluctuations is that given by a stationary, ergodic 
Gaussian random process with zero mean. The 
meaning of this will become clear from the 
following summary of the chief properties of a 
random process. 

The conventional function a r (t) defined in 
the interval (a, b), means that for any value ¢, 


in the range a < t, < b it is possible to calculate 


the corresponding value of 7, namely 27 (¢,). If 


a(t) is a random function (or, more properly, 
a member of a random process) it is no longer 
possible to calculate the exact value a (t,) but 
only to make certain probability statements 
about it. Thus the probability that x (t,) is less 
than some prescribed constant, say 7, may be 
specified by a function F, 


F, (a, = Priz(t) <a; (1.1) 


If F, is differentiated partially with respect to 7, 


fy (ay. ty) day Pria, - 


r(t,) <a, dr,; (1.2) 
for very small dz,. f, is a probability density. 
It is as if a barrier is set up at ¢ t, with a shit 
from 2, to 7, + dx, in it; fy (7%, 4) dx, is then 
the probability that the curve « (f) should pass 
through this slit. For a slit of finite width, 
extending from 2, to 2,’ say, the probability of 


the curve a(t) passing through it ts 


Pria, < a(t) 
f, (é. t,) dé F, 4) (a4) (1-3) 


251 


Statistical analysis of a reactor linear theory 


Higher distribution functions can be set up 
giving the joint probability of two values a (t,) 
and 2 (t,). Thus 
Fy, (2, ty 3 te) = Prjaz(t,) <a, 


and x (tz) <a,} (1.4) 


and if 
OF, 
then 
So ty) da, da, = Pr{a, <2 (t,) <a, 


+ da, anda, (t,) <7, + dry} (1.6) 


for very small dz, and dz,. This corresponds to 
the setting up of two successive barriers, one at 
t, and the other at ¢, and saying that the pro- 
bability of any function passing through two 
slits between 2, and 7, + dr, and 2, and a, + dr, 
respectively is f, dz, dz, The n'" distribution 


function is defined by 


F,, (2, ty V9 ty 
Pr {ax <a, (tg) r(t,) <a,} (1.7) 
>” 
wy d n 43 t,) (1.8) 


A formal definition of a random process may 
now be given. It is a set of functions {2 (t)} 
and set of probability distributions 
1, 2, 3,...such that the 


relation (1.7) is true. The functions F 


are not 
independent, for, to be consistent, the following 
relations must clearly be satisfied 


F’, (00, 0, t,; o,f.) = I 

(24, t3 ©, te ; 00, t,) = F, 


An important class of random processes is such 
that a shift of the time origin does not change 
its description. Such a process is called stationary, 
and stationary processes may reasonably be 
expected when the overall conditions are essen- 
tially unchanged. For a stationary process only 
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Anis and N 
the differences of the ¢ values are significant ; 
thus 


and in particular 


F, (a, ty) 0) (1.11) 


is independent of ¢,. 
Any function of x, é(7), will have a mean or 
expected value 


(2) f, (a, da 


. 


If the 


is stationary /, is independent of ¢ and so is the 


which is generally a function of ¢. process 


mean of any function ¢(@). In particular the 


mean of itself 


IS independent of t for a stationary process. This 
constant may be subtracted from 2 (f) to give the 
class of perturbations with zero mean 

The Gaussian term in the specification of the 
class of perturbations given above is less easy to 
justify. A stationary Gaussian process with zero 
mean is such that all its probability distributions 


f, are normal distributions. Thus 


Sa Ze. ty) 


z,* 
€Xp 


207 (1 


etc.. where o is constant and p a function of 
t, — t, only, o ts the standard deviation of x (ft) 
i.e. the variance 


f, (x) da 


The 


important 


random has the 
that it 


linear 


Gaussian process very 
Caussian 


filter. A 


property remains 


after passing through a linear 


R. 


filter is the name given to a linear differential 
equation or equations for a function y(t) when 


x (t) occurs on the right hand side, that is, as the 


input. For example 


+ Ay(t) =a(t) 


is a linear filter and if the input a (¢) is Gaussian 
so is the output y(t). The analysis will be 
limited to the Gaussian process. 

The final term ergodic is of great importance. 
Speaking loosely it means that any particular 
random function from the set {a (t){ will even- 
tually behave like any other function from this 
this The 


d (x) f, (a, t) da may then be replaced by an 


set. ensemble average at a port 


. 


average over a long section of any one function 


r 


lim r)) dr 


T 


. 


lim (x (t)) dt, 


22 


the limit when the limits of integration 


infinite 


for in 
become it +r can be replaced by t. In 
particular for the functions we consider 
lim =a r(t) dt 0 
T->a 27 
T 
f 2 
oT [ (t)\* dt 
T 


lim (1.18) 


It is the equivalence of ensemble and time 
averages which permits the easy analysis of the 
effect of a linear filter, as in Section 3. 

A most important function for the deseription 
the 


For an ergodic process it is defined as 


of a random process 1s auto-correlation 


function. 
T 


lim [ (t + r)dt (1.19) 


7 


Pr, (7) 


It is an even function of + and if sharply peaked 
at the origin it implies that the value of x (ft) 
is not very strongly correlated with the value 


252 


q 
FP’, 0; te — t,) (1.10) 
| 
Sy = (4) 
r,? 
— 
| 
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but G,,(w) is generally a complex valued 


for slightly different time ¢ r; if it is a spread- ; : 
function and 


out function it means that two values are 


relatively strongly correlated even though at G,,, (w) G,, (@) G,,, (— w) (1.29) 


7] 


two more distant times. It may be thought of 


where the asterisk denotes the complex conjugate. 


as describing in a general way the quality of 


These bare facts on random functions will 


the random fluctuations. From (1.18) and 
(1.19) it is clear that 


sulfice for the discussion which follows: for a 


full treatment the reader is referred to LAanixc 


a” = (0) (1.20) and Barrtin [3). 


and this variance is the most convenient measure 


of the variability of +. 2. THE Eat ATIONS OO] \ REACTOR 
If there are two ergodic random processes (ft) with Ranpom Inpurs 

and y(t) the cross-correlation function The state of the reactor in which a single 
r reaction ts taking place may be described by 

1 f " two quantities, 7, the concentration of one 

Pay \7) 2T (4.31) reactant and the temperature of the contents. 
Referring to the list of svmbols for their meaning, 
may be defined. For this clearly the following equations, one a mass and the 


other a heat balance. ar readily derived. 


da 
Finally a function of great importance is the q (Xo VR (a, 7 2.1) 
spectral density, the Fourier transform of th. 
correlation function. AT 


gh { 


AH) VR (a, T) VU* (T, T 


where 


and inversech V volume of the reactor, 
t time, 
(w) dw. (1.24) q flow rate of reactants. 
i q flow rate of coolant, 
Le inlet concentration, 
but since ¢,,(7) is an even function G,,(«) is R = reaction rate, a function of x and T. 
real, and in fact h product of specific heat and density of 
e . reactant mixture, 
G,, (w) : COS wr dr (1.25 T, inlet temperature, 
T., = coolant temperature, 
and - AH heat of reaction. 
. U* rate of heat removal by cooling 
| G,,(w)coswrdr. (1.26 


The dependent variables z and 7' will be referred 


: to as the output variables, the independent 

G,, (w) Py (7) "dr (1.27) variable ¢ is the time and the quantities g, 
. T., Ye, Which are all functions of t are called the 

b input variables. It is these which will be regarded 

4, (7) 4 an (r) e'@’ dr (1.28) as subject to random fluctuations about their 


steady state values and their influence, severally 
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and N, 


or together, upon the mduced random fluctuations 
of x and T is to be studied. 
Equations (2.1) and (2.2) are rendered dimen- 


sionless by the new set of variables 


qt 


q | 


Th 
All) x, 
The 


and the 


» and « are quantities with mean value 1 


and x are Kas 


mean values of 


The equations become 


When the input variables are held constant at 
their mean values £ and » will approach steady 


state values and »,, which satisfy 


P (€,. »,) 


Subtracting (2.6) from (2.4) and (2.7) from (2.5) 
and 


z and y, namely 


putting &€ OF + gives the 


non linear 


equations for the perturbations 


dz 


y) — P(é,.»,)] 


dy 
dr 


(2.9) 


It is now assumed that the departures from the 


RK. AmMUNDSON 


steady state are small and a complete linearization 


is made. Introduce the following constants, 


(Pe), B= (Un), 


8 (Un),. € = (Ux), (2.10) 


the within the bracket denotes 
partial differentiation and the s without that the 
the 


Also let new random variables with zero mean 


where sullix 


derivative is evaluated at steady state. 


value be defined by 
1+A, & 


(2.11) 


Then, neglecting all squares and products of 


y, A, p and 


and the right sides of these equations are random 
functions of zero mean value ; they will be denoted 
by e(t) and f(t) respectively. 

The 


resolved into two parts, 


problem under consideration is thus 


The first is to study 
the effect of the random inputs ¢ and f on the 
solution of the equations; the second to study 
the composition of ¢ and f themselves from their 


component functions A, yp ete, 


Tur SoLuTion or tHe LINEARIZED 


EQUATIONS 


The equations to be considered are 


dr 


dy 


di . f(r) 


yy az 

Taking the Fourier Transform of these equa- 
tions and denoting the Fourier transform of a 
function by its letter the 


capital following 


solutions are obtained 


iw+1+8, 


E ’ F (3.8 
Alw) + (3.3) 


A(w) 


Fe 
Th 
3) 
Th q Ne 
x 
| 
+ a)2 yy ( E)A 4 (2.12 1958 
dt 
dy 
at 
OV, €A, (2.1%) 
dé 
— €) P (2.4) 
dr 
dy 
dr 
0 Ne P U (n,. * (2.7) 
3. 
| 
(8.1) 
+ (« — 1) + (1 — (2.8) 
KY [ +- 2, y) P(é,.»,) | 
+ (K + K) 
- 
— 
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and } (3.4) 
At w) Af 


where A (w) w* iw(1 x 8 
(1 + «)(1 + 6) 


or w* iwd, d, (3.5 


The inverse transforms of z will take the form 


W,, — badt 
7 


(t)/ r t)dt (3.6) 


The auto-correlation function 


1 
l 
(7) lim oT r) dt 3.7 


7 
7 


may therefore be calculated. It will involve the 


auto- and cross-correlation functions of e and 


Ve. Dees and Similar formula hold 


for ¢,,(7) and ¢,, (7). This will not be carried 
out in detail here as we rely for later calculations 
on the Fourier transforms of these relations. 
The derivation is given in detail in [3]. 

Equations (3.5) and (3.4) can be written in 


matrix form 


Z (w) 
y 


(3.8) 


Ho» (w) 


where H,,(«) is of course the transform of W(t). 
taking the 


transform of the resulting formula the following 


Using (3.6) in (3.7) and Fourier 


relation is obtained. 
2 he Hy, Gy 
Hy» Gy (3.9) 


(w) H,, 2G 


ee 


where G is the spectral density function, an 


asterisk denotes the complex conjugate and 
Re the real part. 


Similarly 
H,,?G 2 Re Hy Gy) 


H,,|* Gy 


(w) = 


(3.10) 
and 


linear theory 


Hs, Hey Ge + Hes Gy 
He, Hu, Gp, + Hig Ge 


je 


(q@) 


(3.11) 


From these formulae the correlation functions 


can be determined by the inversion formulae 


(1.24), (1.28) and in particular 


(w) dw (3.12) 


?.. (0 G.. 


The cross correlations between the input and 
the output can be found by the same technique. 


The relations for the spectral densities are 


Ga (w) = Bt, G, + He 
(3.138) 
G,,, (w) Hy, G,, Hos 
Gy, (w) = Hy, G, + Hy, Gy 


Before functions it is 


necessary to discuss the determination of spectral 


making use of these 
densities of e¢ and f, either empirically or as 
compounded from thy spectral densities of A, 
etc. 


t| Tne ANALYSIS AND COMPOSITION 


OF THE Random INPUTS 


On the assumption that the process ts stationary 
and ergodic the correlation functions of actual 
inputs may be determined from strip chart data. 
Thus, suppose that in a practical case x, (t) has 
been observed over a long period of time, how 
may this record be used to estimate the auto- 
correlation function’ This problem is treated 


in LANtNG and Bartin (Section 4.3) where a full 


derivation is given of the following results 
The function 
1 
V (7) (t) r) dt 
0 


where 7 is quite large is an estimate of the auto- 
correlation function. It is actually a random 
variable itself since it de pends on the particular 
function 2, (t) of the set |x, (t)} which ts the record 
being analysed. However its mean or expected 


value, 


| 
x 
x 
: 
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9 
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7 
E| V (r)] E (t) r)| dt 
ny 
r 
| dt 
7 


0 


is the true correlation function (the sullix 29 2% 
has been omitted for simplicity). 

The variance es," of this estimate of @ may he 
estimated for a Gaussian process and it ts shown 


in the above reference that 


this allows an estimate to be made of the length 
of the strip record which should be used to obtain 
a given accuracy. 

It mav also be shown that a given correlation 
function may be approximated by the sum of a 
series of functions of particularly simple type. 
The function exp(— ¢ |r|) is an even function, 
though its derivative is not continuous at the 
origin. It has, moreover, a particularly simple 
Fourier transform, namely (c* w*), If 
therefore a given correlation function can be 


represented approximately by 


d 
(7) A, exp (— he 


its spectral density is approximately 


G (w) = 4 (4.5) 
» 


For a complete discussion of the representation 
of auto- and cross-correlation functions in this 
way the reader is referred to chapter $ and 
appendix D of Lanxinc and Barriy’s book. It 
will be assumed here that the correlation func 
tions are approximated by the sum of functions 
of the type given in Table 2. The spectral 
densities will then be given in terms of the 


standard forms of this Table and be either 


2 P, iw 


7 a” 
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2 iw® + (aP, w?* 


2ZabQ),) iw + (a® + b*)(aP, 


w* 2(a* w* + (a* 


(4.7) 


these being the spectral densities of P, + 
and Py, + Qe ry + Py be + (see Table 2). 
The final enquiry of this Section must be as to 
how the functions ¢ and f depend on their com- 
ponent parts A, pv ete. From (2.12) and (2.13) 


é (1 (4.8) 
f = (q, A + €A, (4.9) 


Suppose first that the perturbations occur 
severally, e.g. A ja ! v, 0 and only A, is 
non-zero. Then clearly r 0 and f €A, 


so that 


(4.10) 


Similar formulae will hold if or ¢ the 
only active perturbations, but if ja 
0. A + 0 we have 


(%o Da, (4.11) 


If two or more perturbations are active at the 
same time their cross-correlations may be 
involved. The following case will show the 


method to be used in all cases. 


A(t) + p(t); 


p(t r) 


(4.13) 


Since the spectral densities, G, are linear 
transforms of the correlation functions, 4, the 
same relations also hold between the spectral 
densities G,, (@), Gy Gy (w) and their 
components G,,(@) ete. It follows from = the 
approximation technique outlined above that all 
the spectral densities may be approximated by 
expressions of the form (4.6) and (4.7). 
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5. Ture CorrReLvation Functions In general the coefficients in (3.9)-(3.11) may be 
OF THE OUTPUTS written 
It is now possible to show how the formulae Aw’ + Biw + C 
of Section 3 may be used systematically to w* +7 w(d2 2d,) + d/ 


determine the correlation functions of the 


Table 1 gives values of A, B and C for the other 


outputs. The quantities involved in the formulae hick 

(3.9)-(3.11) are of the form H,, (w)/*, Comparing 3.11) 

(3.8) with (3.3) it is easy to see that n each case one of these functions is multiplied 


by a spectral density, a. Gy or Gy, and these. 


following the result of the previous section are 


1+, 
H,, (w)!? 


d, w*) twd | presumed to be approximated by the sum of 
(1 } y)? w* 51 rational functions of w. These approximations 

w* (d* are substituted in formulae (3.9)-(3.11) and the 


Table 


Hy 0 0 y? x (1 x) (1 y) 
Hy, 0 0 x? H* Hy, 0 0 x) 

IH (w)|* 1 0 (i x)? Hye 0 y (1 + «) 


Table 7 Standard corre lation unctions and s vectral de nsities 
/ 


Correlation function Parity Spectral density 


ala* w*) 


Dio 


2 b?) iw 
r, (r) a>o odd (w) 


cos br. +r 7 D 


2 bia® w*) 


D ( w*) 


ig (7) = > 0 even by (w) 


2 2abiw 


(1) ce" sin br. a 0 odd (w) >< 
4 3 
D (w*) 


. 
| 
9 ( KB ( 
58/59 
| : 
Issoc. 
ele 2 
tr, (7) odd (w) 
D («*) 2 (a® b2) 
Py 
| 
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resulting formulae split into partial fractions. 
Four cases arise in doing this according as the 
roots of |A (w)/? 0 are (1) real or (II 


and the spectral ck nsities G ete.. a 


ee 


complex 
involve 
only vy and Yr, or (b) involve the functions 

Tables may be prepared which will 
this process of splitting into partial 
the 


., Inany event, the spectral 


simplify 


fractions and allow a direct calculation of 


coctlicients 
densities of the output are reduced by this means 
to sums of terms like (4.6) and (4.7) and by Table 2 
these may be immediately inverted to give the 
corre sponding correlation functions. 


Phe 


correlation 


variances are obtained from the auto 


functions immediately by setting 


An auto-correlation function ts even and 


of the 


so must always lhe 


Po, (7) + Pi (7) 


the 


corres 


that the 


function take on particular values, 


constants in 
The 


and one 


w he re the cle notes 
ponding variance Is thus 
obvious check on the correctness of numerical 
work is that this should be positive. 


6. AN EXAMPLI 


cond order irreve rsible 


taking 
ft 


\s an example a 
reaction, 24 
tank of 
If the other physical quantities are 
0-3 ft® sec, 1 tb ft®, 

But |b mole, 
ERT) 


lb mole 


products, place ith a 


stirred volume V may be 
considered. 
as follows: 4g mole 
T. 650K. AH 

h 60 Btu ft! K, 2) 
where k, 2°7 
R 1-087 Btu 


with 


20.000 
rk, xp 
Btu 

lb mole “K. Under these con 

a cooling coil of 500 ft® (heat 

coeflicient 100 Btu, br ft? 

with 520 kh. 

ft® sec, the reactor will run stably at T 

0-0843 Ib mol 


ditions and 
transter carrying 
0-138 
824-9 K 
ft? i.e. 


cooling water 
with a concentration 2, 
producing 989 |b moles/hr of product 

In the introduced in 


0-460, 


variables 


1-56, x, 


dimensionless 


1-05, 7), 


section 2 Wp 
UP 2-475, 


x 21-730, 8 


0-0843, and the parameters are 
10-091, — 0-428, 


0-428. 
‘ 1-052 With 
minator of all the functions i, Is 2 w,") 


17-309, w, 


10 ‘ these values the deno 


( w,*) where w 584. 
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Two series of results will be presented using this 
example. In the first the effect of “* white noise ” 
in each of the input variables will be considered 
to show (a) the type of auto-correlation function 
which results for the output and (b) how the 
deviation of sensitive 


standard the output ts 


to the different inputs. In the second the varia- 
tion of the auto-correlation function of the output 
with increasing spread of correlation of the 
perturbed inlet copcentration is shown. 

The 
to an idealized random process in which there ts 
a(t) at 


Since 


term “white noise’ used above refers 


no correlation between the values of 


different instants, Le. (+r) 0. 


dr is tinite im all cases, nitist be 


ae 


0 and it is proportional to the 
delta Naturally 
but it is 


iuifinite for + 

Dirac 
process cannot bn 
the limit of the 


tions are due to infinitely many 


so-called function. this 
realized in practice 
situation in which the perturba 
infinitesimal, 
that its 
spectral density ts a constant and so it is very 
the effect of different 


Accordingly, the following cases will be 


uncorrelated impulses. It has the virtus 


valuable for comparing 
inputs. 


considered : 


(1) Inlet alone is perturbed, 

(2) Inlet temperature alone is perturbed, 
G,, 

(3) Flow rate alone ts yn rturbed, 


concentration 


The results for unit perturbation of this kind in 


coolant temperature and flow rate may be 


those for inlet 


the 


functions by &* and & respectively. 


obtained immediately from 


temperature by multiplying correlation 
In all cases 


the correlation function can be written 


(6.1) 


le, 


the values of the constants C, being given in 
Table 8 together with the 
ao \ C for the autocorrelations. 


standard deviation 


Plots of d(+r) are shown in Figure 1. 
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correlat 


Output ternperature 


Dimensioniess time 


Fic. 1, Plots of auto- and cross-correlation functions for agitated reactor. 
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Table : 


Correlation function Perturbation 


\uto-correlation of Con 


concentration 


Te Thi}? 
Flow 


Auto-c« 


temperature 


rrelation of 


ross-correlation O-OSS 


of temperature and O-oM 


concentration 


O67 OO4S O-259 


0-168 


Owes 


Osi oO 


O-166 Oost 


Finally let the mecentrat perturbed 


a Gaussian random disturbance naving a 


ron tunction 


which the spectral density ts 


G 


whit tha 


the 


The case of 
Figure 2 


the 


nose 1s 


show resuits of these 


bearing of which will now be discussed 


7. Discussion 


The 


column of 


the 


last 


tlow 


deviations 
Table 3 that 


the greatest 


standard given 


show Variations in 


rate caus Variations in the output 
both 
the effect of pe rturbance of the inlet temperature 


the 


concentration and te rature. In 


and of inlet 
the 


rate Is 


concentration are very much 


same. but the effect of pe rturbations in flow 
60-70 
Fig. 1 

shown ; 
functions only half of 


the difference ol 


cent This Is 


the 


son per vreate! 


contirmed in where auto-correlation 


functions are Since these are even 


but 
Here 
again the curves corresponding to per urbations 


of inlet 


each need bn 


scales should bn noted, 


and concentration follow 


Those 

rate 
others, but fall off at much the 
rate 


te mipe rature 


one another closely. corresponding to 


the 
This 


bn judged by saying that when 7 l 


perturbations of flow well above 


same rate. 
may 


i.e., after one holding time of the reactor) d(+) has 


dropped to kk SS than bs per cent of its value at 
0), namely the 


and 


Variance, 


cross-correlation of output te rature 


1. As 


, meaning that 


concentration is also shown in 


would be expected this is negative 
increases of temperature are generally accom 


Also 


revion ol 


panied by decreases in 


the 


concentration 
the 


CTOSS-COrre lation Is 


correlation is for +r in 


This 


vreatest 


Ol. 


linn 
, 
7 
the 


fact that it 


negative 


is greatest when +r is small 


and (say as In 
the products 
This that 
to some perturbation comes about rather rapidly 
Thus 


0-05) because 


curve 2) 


implies that OOS) are 


vreatest. means deviation in y due 


than deviation im d(— 0-O1) will be 


the 
is still small although that in y 


smaller than 4 ( deviation in 


has eTown, 
and on the other hand 4 ( 


| 


0-1) is smaller than 
0-05) because the deviation in y has begun 
to decay although that in z is fully developed. 


Thus the 


that of some small multiple 
of the value of + for which 'b.,, (7) 


response of concentration lags behind 
temperature by 
is greatest. 


In the case of curve 1, corresponding to pertur 


bations of inlet concentration the lag is slightly 


greater. This is because a sudden slight increase 
the 


but 


in inlet concentration will 
the 


which 


mcrease reaction 
the fall 


rise in 


rate and hence temperature ; 


of concentration accompanies a 


0-168 
Fem 0-02) 0-337 
Flow o-109 
: 
VOL 
9 
or)! 1958/ 
b,, bee 0 6.2) 
q@* 
0 to 
leulations 
= z(t) y(t 
and 
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@ (WHITE NOISE) 


5 
z 
° 
= 
oO 
z 
2 


AUTOCORRELATION 


0.6 


DIMENSIONLESS TIME 


Fic. 2. Plots of auto-correlation function of output con- 
centration for correlated input concentration perturbations. 
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temperature is inhibited by the initial rise in 
concentration and so is greatest at slightly later 
time. Curve 3 corresponding to flow rate fluctua- 
tions shows that cross-correlation of output 
temperature and concentration is very distinctly 
smaller than in the previous cases. This may be 
explained as follows: an increase in flow rate 
reduces the holding time and so the extent of 


the reaction. But it also tends to increase the 


AMUNDSON 


NOTATION 


constants 

B 

(1 + a) (1 + 8) - 

random input =< (1 — £,)A 
Fourier transform of ¢ 

random input — 5) A 
ith probability density function 
Fourier transform of 
probability function 


spect censit 
temperature in that the cooling coils remove less pectral density of 2 
cross-spectral density of 7 and y 


heat and this increases the extent of the reaction. 


Fourier transform of W i transform function 


It appears that these two effects to a large degree product of specific heat and density of reactants 


balance one another and keep the pre xiuct of the dimensionless reaction rate 


concentration and temperatur deviations flow rate 

. . reaction rate function 

relatively small. 

time 

In Fig. 2 the output auto-correlation function 
temperature 

orrelated 


of concentration is shown for various heat transfer function, dimensionless 


inputs. The greater spread of the input auto- heat transfer function 


correlation (i.e., the smaller value of c) is reflected volume of reactor 

in that of the output auto-correlation function. estimate of auto-correlation function (§6) 
impulse response 
As a measure of this the values of 4(1 r. the 
concentration 


ratio of the auto-correlation function at + ! perturbation on dimensionless temperature 
to its value at r+ 0, are given below. Fourier transform of y 


perturbation on dimensionless concentration 


Fourier transform of = 


| (P,) 

(P,) 

0-055 | O18 0-55 | O79 

heat of reaction 


$(1)/46(0) input | 0 
(1) output 


0-38 0-61 


c's 
dimensionless temperature 


They show that the 
prolong the range of the correlation of the 


clearly process tends to dimensionless flow rate 


output ye rturbation on « 
marked 


when the input concentration is least correlated. 


perturbation on 


variables and this is naturally most 


perturbation on no 


dimensionless concentration 
autocovariance on z 
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Letters to the Editors 


Tue purpose of this discussion is twofold. Part 1 relates 
directly to the reference article, and presents a comparison 
of SKELLAND’s votator heat transfer correlations and data 
with theoretical predictions based on a simple unsteady 

state heat penetration model as it applies to scraped-wall 
heat transfer. Part 2 presents a similar model as it generally 
applies to a moving impeller, all or part of which may be 
made into effective heat transfer surface by coring and 
circulating a heat transfer fluid through it Although 
this part is not directly related to the article, its inclusion 
here seems appropriate since it does provide a way of 
mereasing the heat transfer capability of equipment 
stnilar in principle to the votator The resulting heat 
transfer gain may sometimes be substantial, e.g., in equip 
ment where it may be impractical to actually scrape the 


wall 


Parr 1 Scrarep-Fitm Hear TRansrer 


Ideally, i.e. under conditions which permit a simple 
and exact theoretical representation, material at a uniform 
bulk temperature continually moves down the rear surfaces 
of the scraper blades to the cylindrical heat transfer surface 
where it either heats or cools by molecular conduction 
until the following blade scrapes it up and thoroughly 
mixes it with the bulk fluid. Since the depth of conductive 
heat penetration per pass is small in the range of r.p.m’s 
where votators usually operatet, the small variation in 
peripheral fluid velocities within the thin heat transfer 
laver may be neglected. It is also assumed that the entrance 
and longitudinal flow effects are rendered unimportant 
by the intensity of the cross-sectional mixing Under 
these ideal conditions, the heat transfer mechanism is 
identical to molecular conduction into a semi-intinite solid, 
where the contact time is the time between successive 
blade passes. Using the surface temperature gradient 
equation in {1}, the average rate of heat influx for a given 
contact time may be readily calculated and expressed in 


terms of an effective scraped-film coefficient, h,} : 


12 
k pe 
h | (1) 


*Chem. Engng. Sci. 1958 7 166. 


Discussion of the paper ** Correlation of scraped film heat transfer in the votator '' (A. H. Skelland*) 


where & thermal conductivity 

density 

Cp heat capacity 

contact tine l N, 

\, number of 

" revolutions hour 

For the votator N, 2. and the final result is 
he 1454 pty : (2) 


Equation 2 represents the idealized situation as described 
» occur in 


above, and deviations would be expected t 
practice. At low Reynolds numbers, for example, the bulk 
mixing intensity would be low and temperature gradients 
Measured 


effective tilm coeflicients ought to be less than pred ted 


would penetrate deeper into the bulk fluid 


from equation (2) At high Reynolds numbers, on the 
other hand, the turbulence intensity would be high and 
eddy penetration of the theoretical heat transfer layer 
ought to give higher film coeflicients than predicted 
Nevertheless, because of the simplicity and theoretical 
aspects of equation (2), it would be interesting and perhaps 
informative to compare it with the empirical correlations 
and data of SKELLAND, whose paper also includes earlier 
data of [3 

Grouping terms in SKELLAND’s correlation (which is 
based on data mainly in the transition region) yields 
h, (D, (107,37) (453 57 47,17) (3) 

~ (yp 37) pc») 2», 17 

It is seen that &, p and ¢, appear approximately to the 
} power as in equation (2), but exponent n is quite different. 
Also, D, and « are relatively unimportant which is in 
agreement with equation (2) v and L, however, were found 
to be quite important. This may be attributed to entrance 
and axial flow effects which were neglected in the theoreti- 
cal treatment. 

If, however, the water data of Hounron (turbulent 
range) are used to establish the coefficients, grouping 


terms then yields 


tEven for water, a fluid of relatively high thermal diffusivity, the calculated maximum depth of molecular heat 


penetration is only about 0-007 in, at 500 r.p.m. in a 2-bladed votator. 


tJurson [5] used this same model in connexion with extruder heat transfer studies. 


coeflicients using a graphical integration of the temperature-distance profiles for various contact times. 


He calculated seraped-film 
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L 37 10 


4 


) 
(2) (k pe, 
L 


Functionally, equation (4) is in very good agreement with 
equation (2) except for the (D, L) term which may be 
effects. 


SKELLAND's and Houwron’s experimental 


attributed to entrance 

In Figure 1, 
seraped-film heat transfer coeflicients are compared with 
theoretical values based on the transfer 
Points on the 45° line w 


simplified heat 


model (equation 2) uld indicate 


perfect: agreement Average properties were used 


to caleulate the theoretical A's according to equation (2) 


h, (EXPERIMENTAL ) 


LEGEND 
LTONS WATER O«Ta (TURBULENT) 
»RELLANODS 
GiYCEROY 


( TRANSITI OW) 


1.6 


(THEORETICAL) 


Comparison of votator data (ordinate) with 


equation 2 (abscissa) 


HovuLron’s water data are in substantial agreement with 
theory and average about 15 percent higher than predicted 
thes 
values would be expected to be higher because of turbulent 
heat 


aspects of the film-penetration theory have been expounded 


Since data are in the turbulent region, experimental 


eddy penetration into the transfer laver. (Various 


by a number of writers [2, 6, 8, 9), just to mention a few.) 
As r.p.m's. increase, one might at first expect the diver 
gence between experimental and theoretical coeflicients to 
But at 
r.p.m’s. the theoretical heat transfer layer is also 
The 


divergence 


increase because of greater turbulence higher 


thinner 
and eddy penetration is more difficult two effects 


apparently compensate, and = the remains 
constant. 
In the 


seems to be 


there 
that the 


data of 
correlation The 


SKELLAND 


fact 


range (oil 
little 


experimental film coeflicients increase 


transition 
very 


with respect to 


theoretical values as r.p.m's. decrease is particularly 


difficult to rationalize. If it may be assumed that the data 
are substantially correct and that the votator operates com- 
pletely full at all times (ie., without boundary separation 
behind the 


I can draw is that the heat transfer mechanism, together 


blades), then about the only conclusion that 


with entrance and longitudinal flow effects, in the transition 
region is very complex indeed, and cannot even be approxi- 
heat transfer mechanism described above. 


mated by the 


Parr 2. Hear Transrer 


The votator will be used for illustration, although the 


following argument applies equally well for various heat 


transfer apparatus. As the impeller turns, fluid flow is 
forces 
heat 
This 


specific set of 


surface by drag and pressure 


that of 


into this moving film 


induced past its 


(see Figure 2). The problem is basically 


penetration from the surface 


problem has already been solved for a 


boundary conditions by Esier [3], and its mass transfer 


Pic. 2. Cross-sectional flow pattern in the votator. 


(Mutator assumed stationary with cylinder rotating.) 


analogue by Kramers and KreyGrer [7] In order that 
ipplicable to the problem at hand, it is 
necessary to assume that: (1) the 


is sufliciently intense so that the fluid is at a uniform bulk 


their solution be 


cross-sectional mixing 


temperature as it starts flowing along the impeller surface 


at point .1, (2) the boundary layer is laminar, and (3) 
the speed of rotation is sufliciently high so that the depth 
\to B 
a linear velocity 


film 


of heat penetration during the fluid travel from 


is small, thus permitting the use of 


gradient. Under these conditions, the coetlicient, 
h;. is given by 
13 
p Cy 
h, o-sos (5) 


where a path length in the direction of flow over the 


heat transfer surface. 
dv dy 


may be an effective or average velocity gradient 


velocity gradient at the surface, which 


if point values vary. 


€p, & as defined in Part 1, 
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In 


operating conditions we may write 


yD 


nD, 
a = dv dy, gn 


D, 


Equation (5) then becomes 


13 


y)nep,k* 
0-808 | (6) 
dD, 


13 


Cp ke 
Constant . | 
D, 


It is sometimes possible to estimate a and x with fair 
accuracy and obtain some idea of the relative heat transfer 
gain which may be obtained by coring the impeller. As 
in the case of scraped-tilm heat transfer, similar deviations 
(depending primarily on Reynolds number) would be 
expected in practice 


COMPARISON with DIMENSIONAL ANALYSIS 


Dimensional analysis, neglecting entrance and axial 


flow effects, would vield equations of the form 


m 


AD ),2 
1D, constant | (7) 
hk 


Kquations (2) and (6), respectively, may be written in 


similar form 


hn (" D? “) (2’) 
k 

AD 

k k 


where (nD? p) is the “mixing vessel” Reynolds number 
since the impeller and tube diameters are the same 

Equations (2°) and (6’) would have to be modified in order 
to take into account departures from the idealized flow 
models assumed in the derivations. A Reynolds number 
correction would account for differences in the character 
of the fluid motion near the heat transfer surfaces. The 
effect of Reynolds number, in turn, would be influenced 
by the thickness of the theoretical heat transfer laver 
This thickness depends partly on the same variables which 
determine the Reynolds number and partly on the thermal 
diffusivity of the fluid, and may be characterized by the 
Peclet number. These corrections, respectively, would be 
of the form 


nD? p 
constant | 
and, 
4 8 B 
D? pc, \' nD? p\" 
k 
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geometrically similar equipment under similar 


Equations (2°) and (6°) then became 


A.D, 


constant, . 
k 


Equations (2) and (6 are identical to equation ( 


except that the coeflicients m and n are expressed in terms 
of components which characterize the reference states 
(as assumed in the derivations) and deviations therefrom. 
Admittedly, this type of analysis may not be generally 
useful because of the interactions between x and 8. but 
it may, in some situations, allow a somewhat better 


appraisal of the heat transfer mechanism. 


Research Department G. A. LATINEN 
Monsanto Chemical ( ompany 
Springfield Vassachusetts 


NOTATION 


dD, internal diameter of votator tube ft 
N,, number of votator blades 
L length of votator tube tt 


h scraped film heat transfer 


BTU hr ft 


coellicient 


Her surface heat transfer 

coetlicient BTU hr ft 
" mutator shaft speed rev he 
0 1 contact time of a fluid 

element with evlinder heat transfer 

surface hr 
h thermal conductivity BTU hr ft* ft 

fluid density Ib 
Cy fluid heat capacity BTU Ib F 
t average axial flow velocity of fluid 

through votator ft hr 
absolute fluid viscosit, Ib ft hr 
a dv dy average surface velocity 

gradient along impeller heat transfer 

surface ft hr 
a path length in direction of flow past 

heat transfer surface ft 

constant dD, 
constant an 


exponents (see text) 


1/2 + By 1/2+ py 
AD 
constant, 
n 
1/3 + ae 1/3 2 
| 
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Notes on the diffusion-type model for longitudinal mixing in flow (Lievenspiel, Smith and Van der Laan) 


Tur value of Lievenseie. and recent paper [1] 
on diffusion models of mixing has been greatly enhanced 
by Van per Laan’s letter [2]. May I ask the courtesy of 
your columns to add two points which may be of practical 
interest. 

It is in principle impossible to inject a delta function 
of tracer into the flowing system, though in practice the 
time of injection may be very short indeed. However no 
great effort need be expended on devizing an injection 
time, is 


svstem if the concentration, as a function of 


measured at two points in the bed. In this case it is the 
difference of means which is nearly proportional to the 
distance between the points and the difference of variance 
effective 
Consider the equations (1) of Van der Laan’s letter but 
hamd side of 
» and measurement 
and z,, are points at which the concen- 
takes 


concen- 


which gives the diffusion coeflicient 


without the term on the right 
(ib) 


at that both 


injection 
Suppose, instead of injection at = 


function is measured and injection 


of the 


tration-time 


place upstream of the point =. The form 


tration function injected does not matter, it suflices to 
know that at = t_ the me asured concentration function 
is R, (o) say. It is now only necessary to solve Van der 
Laan’s equations (1b) (with right hand side zero) and (le) 
subject to the conditions (3b), (8d) and (3f) 

R, (3), R(z,0) = 0,2 > %. This 


obtained by the Laplace transformation and at = = 


solution is readily 


gives 


Ry 


4) Pe 


Ry (P) = R (ayy P) 


q+ Pec 


(1 + © 


q © 
q—q)e~ 


(1 rq” 


where { = z, — =, and the other symbols are as given in 
Van der Laan’s letter. If this expression is expanded in 
powers of p the coefficients give the moments as Van der 
Laan shows in equations (5) and (6) ; the conditions under 


which this may be done I have recently examined in more 


detail elsewhere [3]. We find that the expression in the 


9 
un) pe 


square brackets may be written | — up + § (o* 


where 


Now if Ry (p) 1 Ho P } (og? + po?) p®... so that 


Hy and o* are the mean and variance of the concentration 


time curve at the measuring point then 


R,, (p) 1 (Ho + + p 


em P (o,,° tm”) P* 


Thus the difference in means is ,,, Ho p and the 

difference in variances o,, 2 Oo" 
0 

Firstly they are independent 


Two remarks may 
be made concerning « and 
of the conditions at = — 0; all we require is that the 
injection of tracer should take place up stream of = 
taking the 


sufficiently far from the end z 


Secondly by second point of measurement 
z, the terms introduced by 
the conditions at this end may be made negligible. The 
terms are in any case so small that an iterative refinement 
of the approximate value Pe 2 o* can easily be made 

The second point concerns the calculation of ao" and Cn” 
from the observed concentration time curves. Firstly, it is 


clearly best to calculate as | (9 — R(b)d rather 


than as | # R(d) dd — p*, since the latter may lead to a 


0 

small difference of large quantities. Secondly, in practical 
evaluation of these integrals it will be necessary to stop 
the integration at some finite #, and it would be useful to 
estimate the error involved in this. The normal distribu- 


tion will provide a suitable, though rough, estimate ; 
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it may be shown that for very long columns the normal up to time @ so that 7 corresponds to @ — # and 


distribution is approached [3]. rt =(0 — #)/o' then + can be obtained from the following 
The normal distribution with zero mean and variance Table: 
o* is 


Table 1 


N (t) 


T 


and |f N(t)dt = exp N 2425 1-064 


30009 10005 
which clearly tends to o? quite rapidly as T+ «. If 


aT oT 


Thus 1°, accuracy is attained if r’ is approximately 3. 


R. Anis 


University of Minnesota 
Minneapolis 14. Minn. 


Suppose then that o* has been calculated by integration 
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Mass transfer to growing drops 


IN recent 


vears there have been some investigations of for extraction during the contraction of the droplet. 


the mechanism of diffusion into a drop growing at a \ consideration of the theory indicated that 70-6 per cent 
nozzle. COULSON and SKINNER [1] have extracted organic — of the extraction took place during drop growth, i.e. in half 
acids from water with benzene droplets, which were formed the “ drop time,” so the values of m and ¢ were altered 
at a constant rate and then pushed back into the nozzle accordingly. In view of the partition coefficients it was 
at the same rate before they could break away. The assumed that the transfer of the propionic acid was 
absorption of pure sulphur dioxide by growing water 


drops has been measured by Grooruuts and Kramers [2 
The results of these workers may be compared with an 


equation proposed by InKovic [3] and fully derived by lili st 
MacGittavry and Ripean [4]. The equation was used * 
to estimate the diffusion-controlled current in a dropping- 

mercury cathode The model for the derivation is a ; 

sphere expanding at a constant volume rate, with short * 7 

range diffusion of a solute towards its surface. where the 3 
concentration is zero. The mass of solute that has diffused ‘ 


in time ¢ after the appearance of the sphere is given by 


B57¢ Di VIA (1) 


Since the derivation assumes that the penetration depth 


is small compared with the sphere radius, the equation is 0-20} | | | 
valid for diffusion in either phase, so long as the resistance 

to transfer lies in one phase only. OO® aks 5.7 050 


\ representative selection of experimental results [1, 2] 2 


is plotted on Fig. 1, in comparison with equation (1). 


CouLson and Skinner's results were modified to allow Fic. 1, Comparison of equation (1) with published data. 
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phase and the recently found rates of gas absorption [6] and extraction 


controlled by diffusion in the benezene 
by forming drops are several times those predicted by 


transfer of the benzoic acid was controlled by diffusion 7 
equation (1). 


in the water phase. In the case of the sulphur dioxide 
Department of Che mical Engineering M. H. 1. Baten 


absorption [2], the driving-force concentration was taken 


as the solubility of the gas Tennis Court Road 


The points plotted in Fig. 1 show a mean deviation Cambrids 


from equation (1) of IS per cent. This may be due to 
heat effects (in the sulphur dioxide absorption), circulation, NOTATION 


drops. Nevertheless solute concentration in bulk of continuous phase, 


or the non-spherical shape of the 
the equation provides an approximate method of predicting referred to « O at interface 
diffusion coeflicient of solute 


mass transfer rates to growing drops 
mass of solute transferred during drop growth 


that mass transfer during the 


It has been pointed out [5 
its vibration time of drop growth 


formation ofa drop may be largely due to 


just after it breaks off. This is confirmed by the fact that volume of drop after time f 
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\. Guemer: Gleichgewichtsgase der Verbrennung 


und Vergasung. Wiarmetechnische Berechnung. 


Springer, Berlin, 1958 148 Seiten mit SS Abbildungen 


Buch hat der Vertasser nicht nur die bisher 


schon vorlievgenden Erkenntnisse aus zahlreichen ver 


IN diesem 


streuten Schrifttumsquellen in willkommener Weise zusam 
menyefasst, sondern er zeigt dartiber hinaus die Moglichkeit 
vemeinsamer Berechnungsverfahren fiir feste, flissige und 
vasformige Brennstoffe Als Grundlage dienen allvemein 
viltige Gesetze, nach denen Verbrennunys-und Vergasungs 
vorgiinge abhiufen ks werden eine gvrosse Zahl von 
Diagrammen gebracht, die zum Teil bisher noch nicht 
veroffentheht waren Sic yvestatten mit tragbarem 
\Nufwand auch schwierivere Probleme EKinblick zu 
vewinnen 

Der Stoff ist in cinen umfangreichen Lhauptabschnitt 
die Verbrennung und einen zweiten von geringerem 
Umfang tiber die Vergasung eingeteilt. An den ersten Teil 
schliesst sich ein Anhang an, der vorzugsweise die in den 
Text vielfach eingestreuten Beispiele ibersichtlich ordnet 
Auf dem eimfachsten Fall einer vollkkommenen Verbrennung 
mit Luftiberschuss ohne Dissoziation aufbauend, wenden 
sich die Betrachtungen der Verbrennung unter Luftmangel 
sowie der Verbrennung bei Luftiiberschuss mit ¢ O, und 
H.O-Dissoziation zu ound fiihren in den folgenden 
Abs hnitten auf die bei der Dissoziation im Spiel befind 
lichen Teilgase sowie auf die simultanen Gleichungssysteme 
zur Berechnung der Teildriicke mit den Gleichgewichts 
konstanten. Uber die Enthalpien und Entropien det 
beteiligten Gasmischungen yveben Ht Ix und IS 
Diagramme Aufschluss, die sich u.a. an die bekannten 
Begriffe von Mollier anlehnen Mit andern nambhaften 


Fachleuten, wie Bosnjakovie hat der Vertasser bei der 


Bearbeitung in Verbindung vestanden und so sein ges 
chlossenes Bild des wissenschaftlichen Standes der Dinge 
veschaffen Schliesslich geht er auf die technischen 
Arbeitsfiihigkeiten der auf di 
Umsetzung cines Wirmegefilles in AusstrOmyeschw indig 
keit mit Hilfe von Diisen und auf die Thermodynamik 
chemischer Prozesse ein und leitet so zu wichtigen 
\nwendungsfragen und dem zweiten Hauptabschnitt ber, 
der im wesentlichen zwischen Kohlenstoffvergasung und 
Brennstoffvergasung unterscheidet und ebenfalls reich 
haltig mit Zustandsdiagrammen ausgestattet ist. 
Ausklingend kommt die Brennstoffvergasung mit und 
ohne Beriicksichtigung einer Methanbildung zur Sprache. 
Fir die Ermittlung der bei den Vorgiingen mitwirkenden 
Teilgasmengen ist die Grenztemperatur massgebend. Das 
Krmittlungsvertahren und die Beziehunyen, denen dic 


Grenztemperatur unterliegt, werden angegeben. 


Hi. Picnier 
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Electrostatics in the Petroleum Industry. 
Edited by A. KuinkenBerc and J. L. VAN per 


Elsevier, 1958 191 pp. $7.00. 


Tuts most interesting book is concerned with the hazard 
of explosion and fire, arising during the handling of in 
flammable liquids, by reason of the generation and discharge 
of static electricity. Its particular purpose is to present the 
results of certain research and development work by the 
Shell group of companies against a background of general 
knowledge and experience in this field 

The main material is in fact presented three times over 
first, very briefly in an introductory chapter (4 pages) 
second, ina somewhat qualitative manner, but with duc 
regard to the sort of quantities that practical people need 
to know, in three succeeding chapters (32 pages) : and the 
third, with extensive theoretical background, in the space 
of five more chapters (84 pages). The first two presentations 
make up Part I of the book ; and the second occupies a 
considerable proportion of Part Il. The former will appeal 
more to the practical man concerned with storage and 
handling problems, and the latter to the research worker or 
student desirous of acquiring a deeper understanding of 
the subject than is perhaps necessary for day-to-day 
applications. Whether either would be quite happy about 
having to pay for that part of the book in which he is less 
interested is uncertain. 

Both interests are however combined in the later 
chapters of Part Il. One of these deals with practical 
methods of measuring the relevant electrical properties 
of liquids (charging tendency and conductivity) and the 
electrification of liquids in actual equipment (charging 
current and resultant field strength). Two other chapters, 
in still more readable vein, describe field experiments on 
the veneration and discharge of hazardous amounts of 
static electricity in storage tanks and fuelling equipment. 

The theoretical chapters of Part II are documented 
fairly fully. There is also, at the end of the volume, a short 
bibliography which claims to cover general surveys on static 
electricity. Some of the references included at this point 
are however of rather specialized significance, and they, 
having already been given at the end of the appropriate 
chapter, might well have been omitted here. On the other 
hand, surveys should perhaps have been included which 
extend the discussion to static hazards in the related fields 
of textiles, explosives and powdered materials generally, 
for these naturally receive little or no attention in the main 
text. 

The principal thesis of the book can be briefly stated. 
It is that the solution to the problem of static electricity 
in inflammable liquids is to be found along two lines which 


must be regarded as complementary. They are : 
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The complete bonding and earthing of plant and 
equipment. 


The improvement of the electrical conductivity of 
the liquid, where necessary by the incorporation of 


effective additives 


The latter have recently been the subject of extensive 


research and highly effective additive mixtures are an- 
nounced. It is argued, quite cogently, that by increasing 
the conductivity of a product in this way immediately after 
manufacture are safeguarded 


all subsequent processes 


provided that carthed equipment is employed. This being 
so, pumping speeds, for example, need no longer be re 
stricted 

Although the book is concerned almost exclusively with 


the petroleum industry so far as application is concerned, 


the principles are relevant to other industries, especially 


the chemical. Indeed, all who are concerned with handling 


low-flash-point, low-conductivity liquids on any appreci 
able 


has to say 


scale should consider the implications of what it 


J. H. 


Tierische Leime und Gelatine 
Vil 335 pp. DM 40-50 


kK. SAUER 


Verlag. Berlin, 1958 


Springer 


Fecknik in 


von 


Der rasche Fortschritt von Wissenchaft und 


den letzten 10 Jahren erfordert die Herausgalx 


Monographien, da es selbst dem Fachmann schwerfallt 


Neuerscheinungen in 
las 


dem Laien aber unmdéglich ist, alk 


den verschiedensten Zeitschriften zu lesen vorlie- 


gende Buch will einen Uberblick tiber den heutigen Stand 
der Technologie von Leim und Gelatine geben. 

Die Einteilung des Buches ist tibersichtlich geordnet. 
Nach einem kurzen geschichtlichen Abriss folgt ein Kapitel 
liber die Struktur des Kollagens und Glutins und den 
Vorgang der Umwandlung der beiden Substanzen ineinan- 
der. Weiterhin finden sich viele Angaben tiber die physi- 
kalischen Eigenschaften von Glutin, wie z.B. Viskositat, 
Oberflaichenspannung, Quellung und optische Eigenschaften 

Verwendungsarten des Glutins 


Fir die verschiedenen 


sind kurze Beschreibungen vorhanden clie 


Theorie des Klebevorgangs, Gelatine fiir photographische 
Zwecke 


Sehr ausfiihrlich mit vielen Skizzen von Maschinen und 


Spezi Ihe 


Speisegelatine und technische Gelatine 


Anlagen ist das Kapitel iiber die technische Herstellung 
von Glutin ausgestaltet, wobei noch besonders die Priif 
und Analvysenmethoden zu erwihnen wiren 

Neben IHlaut 
noch Fischleim und Kasein ausfiihrlich behandelt 

Als Nachteil wird beim Studium des Buches empfunden 


Knochenleimen wird 


und speziell 


dass an vielen Stellen wie z.B. bei der Beschreibung der 


Quellung, Gelbildung, Gerbung usw. nicht mehr an neuerer 
Literatur beriicksichtigt wurde 

Trotzdem kann gesagt werden, dass es dem Praktiker in 
Tabellen 


Schaubildern manchen Hinweis geben kann, jeder aber 


der Glutinherstellung mit vielen und 


der sich iber Gelatine und Leime unformieren will, erhalt 


einen umfassenden Uberlick 
Mou 
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